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Abstract 



The low-energy sector of Quantum Chromodynamics (QCD), which is domi- 
nated by a very strong interaction between the quarks and the gluons, is an 
excellent example of how an analytical and a numerical approach can work 
together and thereby mutually increase their predictive power. 

Chiral Perturbation Theory (xPT) offers a systematic method to predict 
analytical dependencies of different physical quantities in the low-energy region 
of QCD. It is, however, an effective theory which is based on the - generally 
accepted - assumption that chiral symmetry is spontaneously broken in massless 
QCD and relies on data from other sources to fix the low-energy constants in 
its Lagrangian. 

Lattice QCD, which is the only non-perturbative regularization of QCD 
known at the moment, is heavily based on numerical methods. It offers the 
possibility to calculate quantities in QCD from first principles. But it has shown 
to be a very delicate problem to incorporate chiral symmetry on the lattice. 
A fact which lead to many technical and fundamental problems. The recently 
rediscovered Ginsparg-Wilson relation is equivalent to chiral symmetry on the 
lattice and with this a long standing problem in Lattice QCD was solved. Hav- 
ing chiral fermions on the lattice leads to a much better control in calculations 
where chiral symmetry is important and allows to test the assumption of the 
spontaneous breakdown of chiral symmetry in QCD in a cleaner way than with 
traditional lattice fermions. 

Fixed Point (FP) Dirac operators are designed to solve another big problem 
in lattice simulations, the discretization errors introduced by putting a theory 
on the lattice. Because FP Dirac operators have been shown to satisfy the 
Ginsparg-Wilson relation, they combine two important virtues one can wish a 
lattice Dirac operator to have. 

In this work we construct a parametrization of a FP Dirac operator and ap- 
ply it in quenched Lattice QCD. The symmetry requirements for a lattice Dirac 
operator are discussed and an efficient way to make a practical construction 
of general lattice Dirac operators is provided. We use such a general lattice 
Dirac operator to approximately solve the Renormalization Group equation that 
defines the FP Dirac operator in an iterative procedure. We discuss the prop- 
erties of this parametrization and show that its breaking of chiral symmetry 
is much reduced compared to the most frequently used lattice Dirac operator, 
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Abstract 



the Wilson Dirac operator. Furthermore, we discuss the overlap construction 
with the parametrized FP Dirac operator. With this construction the remaining 
chiral symmetry breaking of the parametrization can be reduced systematically 
to a desired level. We show that the locality of the overlap construction with 
our Dirac operator is improved with respect to the common overlap construc- 
tion with the Wilson Dirac operator. Using the overlap construction with the 
parametrized FP Dirac operator we perform several test calculations, where a 
precise chiral formulation of the lattice Dirac operator is needed or desirable. 
Using the Atiyah-Singer index theorem, which is a consequence of chiral sym- 
metry, we calculate the quenched topological susceptibility. The results are 
consistent with other recent determinations, but they are not accurate enough 
to make a controlled continuum extrapolation. We also confirm recent results 
about the local chirality of near-zero modes of the Dirac operator. Finally, we 
determine the low-energy constant S of quenched Chiral Perturbation Theory, 
which in full QCD with 3 light flavors would correspond to —(■tpip), the order 
parameter of spontaneous symmetry breaking. 



Chapter 1 

Introduction 



Over the last 30 years Quantum Chromodynamics (QCD) has been very suc- 
cessful in explaining many phenomena of the strong interactions. Using stan- 
dard perturbation theory, which was developed in the framework of Quantum 
Electrodynamics (QED), processes involving high energies can be studied in a 
systematic manner. However, at energies below 1 GeV the strong coupling con- 
stant a s gets of 0(1) such that standard perturbative methods fail - a fact 
which has many consequences. Not only that the interesting physics of parti- 
cles like protons and neutrons, out of which nuclear matter in our surroundings 
consists, is governed by the non-perturbative properties of QCD, but also that 
confinement of quarks, a profound property of QCD, makes that in experiments 
always bound states of quarks and gluons - the so called hadrons - are ob- 
served. This implies that in order to compare the perturbative calculations with 
experiments one typically needs information about the low-energy behaviour of 
QCD, e.g. in form of hadronic matrix elements. 

Certain non-perturbative problems of QCD can be treated analytically using 
an effective theory for the Goldstone bosons degrees of freedom. This approach 
assumes that chiral symmetry is spontaneously broken This spontaneous 

breaking of chiral symmetry leads to the appearance of massless excitations - 
the Goldstone bosons - if the mass of the quarks is set to 0. In nature the 
masses of the lightest quarks are not 0, but a few MeV for the up and down 
quark and roughly 120 MeV for the strange quark. Hence, chiral symmetry is 
broken also explicitly by the quark mass terms. It however shows that these 
masses are still small enough compared to the intrinsic scale of QCD, Aqcd, 
that the effect of the spontaneous breakdown of chiral symmetry dominates the 
low-energy structure. Chiral Perturbation Theory (^PT), which is a systematic 
expansion in the quark masses and momenta of the light quarks, has been 
very successful in predicting many of the low-energy properties of QCD, such 
as e.g. quark mass ratios, scattering lengths and, in general, properties of the 
Goldstone bosons (7r,K,K,?7). However, being an effective theory the predictions 
of xPT rely on accurate determinations of low-energy constants, which have to 
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be extracted from various experimental data. For certain of these constants this 
shows to be rather difficult, if it is possible at all. Furthermore, x^T works well 
only within a limited range of energies, which moreover may vary from process 
to process under consideration. 

The lattice approach to QCD is a radically different way[] to have access 
to the low-energy structure of QCD, as it gives a completely non-perturbative 
definition of QCD. The idea, which has been proposed by Wilson in 1974 0, is 
to discretize the QCD action by replacing the continuum space-time by a discrete 
four-dimensional lattice with a finite lattice spacing a. It is the only definition 
of QCD beyond the perturbative level known at the moment and it provides, 
at least in principle, a tool to study QCD from first principles. Calculations 
in Lattice QCD involve the numerical evaluation of the path integral defining 
the discretized QCD action by stochastic sampling (Monte Carlo) techniques. 
The way from the point where one writes down the discretized QCD action 
to the point where one can extract accurate physical data which is relevant 
for QCD phenomenology has, however, shown to be more tedious and longer 
than expected in the early 1980's. In order to understand this we focus on 
two fundamental problems in lattice QCD: Chiral symmetry and discretization 
errors. 

In defining a disretized version of the gluonic as well as the fermionic part 
of the Lattice QCD action one has, in fact, an infinite choice of discretizations. 
The discretizations which are the most popular are the Wilson gauge action 
and the Wilson Dirac operator. These actions define the simplest way how the 
theory can be discretized with the right continuum limit a — » and essentially 
without destroying the subtle structure of the anomalies of QCD. The Wilson 
Dirac operator, however, in order to avoid the problem of species doubling — an 
effect we will discuss below — introduces a term that breaks chiral symmetry 
explicitly. This so-called Wilson term is an irrelevant operator and therefore 
chiral symmetry is restored in the continuum limit, but the costs are high. 
Due to the breaking of chiral symmetry the Wilson operator introduces errors 
of 0{a) into the lattice discretization and, furthermore, the quark mass gets 
additively renormalized, operators in different chiral representations get mixed 
and additional renormalization factors have to be calculated. All these problems 
are essentially of technical nature, but they make the extraction of physical data 
at least very tedious, thereby hiding the underlying chiral structure of QCD. For 
chiral gauge theories like the electroweak sector of the Standard Model, however, 
the explicit breaking of chiral symmetry has so far precluded all attempts to 
achieve a realistic lattice formulation. Another approach to put fermions on the 
lattice, which is used very often in simulations, are the so-called staggered or 
Kogut-Susskind fermions ||. They partially solve the problem of the species 
doubling by reducing the number of doublers from 16 to 4 and thereby also 



J QCD sum rules are another technique used to get non-perturbative results in QCD. 
As we do not use any sum rule results in this work, we refer the reader to the review 
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solve the problem of chiral symmetry breaking to a large extent by leaving a 
chiral U(l) subgroup invariant. However, they introduce a new headache: they 
break flavor symmetry, which makes the construction of hadron operators with 
the correct quantum numbers rather difficult. Even though staggered fermions 
have only 0(a 2 ) discretization errors, these errors can be large [g], making it 
necessary to go far to the continuum limit in order to get physical results that 
are not contaminated too heavily by the discretization^. 

In order to reduce lattice artifacts Symanzik has proposed the idea that one 
can add irrelevant higher order terms to the initial action such that the leading 
cut-off effects are cancelled [11, 12|. This idea is implemented in the Luscher- 



Weisz action 1 13 ] and the Sheikholeslami-Wohlert (clover) Dirac operator [1^ 
Similarly, other improvement programmes were pursued, most prominently the 
non-perturbative clover improvement, which eliminates all the O(a) effects 
This programme has shown to be very successful in the determination of hadron 
masses, which show a very nice scaling behaviour, whereas for other quantities 
the 0{a 2 ) cut-off effects can still be large. 

This brings us to the main topic of this work, the fixed point approach to 
QCD. Wilson's Renormalization Group (RG) approach offers a radical way to 
treat the problem of lattice artifacts fll6|-|l8|]. Using so-called quantum perfect 
actions one gets entirely rid of the lattice artifacts. Such actions are unfortu- 
nately very difficult to approximate and to really solve the defining RG equations 



is out of reach. Hasenfratz and Niedermayer have shown |19fl that for asymp- 
totically free field theories - like QCD - it is possible to define a so-called 
classically perfect or fixed point (FP) action that reproduces the properties of 
the corresponding classical theory without discretization errors. In fact, this 
corresponds to an on-shell Symanzik improvement at tree-level to all order in a. 
For asymptotically free theories, where the continuum limit corresponds to the 
limit where the coupling goes to 0, such a FP action is expected to be a close 
approximation to the quantum perfect action even at non-zero coupling. 

The FP Dirac operator -Dpp does not only have reduced cut-off effects, but 
as pointed out by Hasenfratz in 1997, it satisfies the Ginsparg-Wilson (GW) 

{©,751 = 2)752^®, (1.1) 



relation [|0| 



or equivalently 

{Sr\ 75 } = 752 J R, (1.2) 

where £> is the lattice Dirac operator and R is a local, hermitian operator which 
commutes with 75. The GW relation is equivalent to have chiral symmetry 
on the lattice, as Luscher has shown in 1998 |21J| . Hence, FP actions offer 

2 There are several promising attempts to improve the cut-off effects of staggered 
fermions and to reduce the breaking of flavor symmetry (see references in fic||). 
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a solution to two of the fundamental problems - chiral symmetry and cut-off 
effects - which have been making the extraction of physical results from lattice 
simulations in QCD difficult over the last 20 years. Let us mention at this point 
that even though this work is concerned mainly with the chiral aspects of the FP 
approach, i.e. its application to very light quarks, the FP QCD action has also 
the potential to be used in lattice simulations with heavy quarks - especially the 
charm quark - because these simulations, which are performed at momentum 
scales near the cut-off, are particularly subject to distortions by the discretized 
action. The application of FP actions to this topic, however, remains unexplored 
also in this work. 

There are other (approximate) solutions to the GW relation which are es- 
sentially focussed on the aspect of chiral symmetry, such as e.g. the approach of 
Gattringer, Hip and Lang to approximate a solution of the GW relation by a sys- 
tematic expansion in operators built out of paths up to a certain length [22-24]. 
The most prominent examples are, however, the 5 dimensional domain-wall 
(DW) fermions [25-27] and the related overlap fermions p^-30], which are 
both based on a proposition made by Kaplan in 1992. For more informations on 
these approaches we refer to the recent reviews in [f3~l]-33| for the DW fermions 
and to [ 33 1 for the overlap fermions. Both constructions, in contrast to the FP 
Dirac operator, are defined in such a way that it is more or less straight forward 
to implement them in simulations, even though they are numerically very ex- 
pensive. In fact, the costs of simulations with (approximately) chiral fermions 
in quenched QCD are roughly the same as full QCD simulations with Wilson 
fermions, i.e. that it is O(10) — 0(100) times more expensive than quenched 
simulations with Wilson fermions. 

Most of the simulations with (approximately) chiral fermions and also the 
largest ones have been performed with DW fermions. Long standing problems, 
where chiral symmetry plays an essential role, have been tackled. Examples are 
the weak interaction matrix elements, such as those relevant for kaon physics, 
i.e. the 5-parameter Bk, the AI =1/2 rule and the parameter of direct CP 
violation e'/e ]34|-36|. The results of these very intricate calculations, however, 
are clearly a deception, as the sign of e'/e measured on the lattice is opposite 
to the experimental value and for the magnitude of relevant parameter of the 
AI = 1/2 rule the two large lattice simulations (RBC, CP-PACS) differ roughly 
by a factor of 2 from each other. This raises the question where these problems 
come from. The answer is not definitively clear, but it has shown that it is rather 
difficult to have good control over the remaining breaking of chiral symmetry in 
the DW fermion approach, even though the linear extent of the fifth dimension 
N s can be used to extrapolate to the limit N s — > oo, where the DW fermions 
get exactly chiral. But there is evidence that the convergence is rather slow [32]. 
There are propositions how this problems can be solved [37, 33] and also the 



overlap fermion approach offers very good possibilities to control the remaining 
breaking of chiral symmetry. Future simulations will show whether the problems 
of the DW approach are really related to the remaining symmetry breaking or 



7 



whether other possible explanations as e.g. discussed in [^,35] are more relevant. 

The definition of chiral fermions on the lattice has yielded the hope that 
finally chiral gauge theories can be formulated on the lattice and, indeed, the 
last years have seen a lot of progress in this area [21, 39-|62] . There remain, 
however, several open questions, like the relative weight factors between the 
topological sectors, the fate of CP symmetry and a related question about the 
definition of Majorana fermions |9, 53, 54|]. Early hopes that a GW like approach 
might work for supersymmetry also did not realize either. 

Before we embark upon a more detailed discussion of chiral symmetry in 
QCD and chiral fermions on the lattice — in particular fixed point fermions — 
let us give an outline of this work. 

In Chapter 2 we discuss the structure of a general lattice Dirac operator 
that satisfies all the basic symmetry conditions, which are gauge symmetry, 
hermiticity condition, charge conjugation, hypercubic rotations and reflections. 
We give examples of terms that can occur in such a general operator, which are 
specific for our parametrization of the FP Dirac operator, and show how one 
can use them in a practical application. 

The details of the parametrization of the FP Dirac operator are discussed in 
Chapter 3. The parametrization turns out to be quite difficult and a lot of effort 
has been put into this part of this work, because the results in the simulations 
clearly depend on the quality of the parametrization. 

In Chapter 4 we give a summary of the properties of the parametrized FP 
Dirac operator. We mainly focus on the chiral properties, however, also scaling 
properties are discussed to a certain extent. The main source of data on scaling 
properties is the hadron spectroscopy study which is presented in the PhD thesis 
of Simon Hauswirth [65]. 

Chapter 5 explains Neuberger's overlap formula and its application to the 
overlap construction with the parametrized FP Dirac operator. Furthermore, 
properties of this specific overlap operator are discussed. In particular, we show 
that our overlap construction is more local than the usual construction with the 
Wilson Dirac operator. 

In Chapter 6 we present results for the quenched topological susceptibility 
and the local chirality of near-zero modes. 

The determination of the low-energy constant £ of quenched Chiral Per- 
turbation Theory, which in full QCD with 3 light flavors is equal to the order 
parameter of spontaneous chiral symmetry breaking (tpijj), is given in Chapter 
7. 

Conclusions and prospects are given in Chapter 8. 

Finally, this thesis covers only part of the work that has been done in col- 
laboration with Simon Hauswirth, Kieran Holland, Peter Hasenfratz and Ferenc 
Niedermayer in an ongoing project, where the parametrized FP Dirac operator 
is tested and applied to various calculations in quenched QCD. Parts of it have 
already been published in [|6^-|69| and much more information can be found 
in 1651. 
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1.1 Chiral Symmetry in QCD 

In this section we give a short overview of various aspects and consequences of 
the symmetries of QCD. We will, however, not discuss interesting properties of 
QCD, such as asymptotic freedom, confinement and various other topics and 



refer the reader to one of the many books about QCD e.g. [70, 7|]. We will 
keep the discussion in the framework of QCD, because the discussion of the 
same topics in quenched QCD is heavily loaded by various technicalities that 
rather hide the structure of chiral symmetry and its consequences. We will, 
however, mention the differences between the quenched approximation and real 
QCD there where it is needed and we will, in particular, point to the differences 
in those chapters where we measure quantities in quenched QCD. 

Symmetries are a very important concept in physics since their presence 
always simplifies analysis and in certain cases allows one to obtain exact or 
semi-exact results. In QCD the presence of light quarks are the reason for an 
approximate symmetry and this allows one to extract a lot of consequences 
concerning the dynamics of the theory. 

In Euclidean space, the Lagrangian of QCD with the gluon field strength 
tensor F^ v , the quark fields tpf, the Dirac operator !pand the gauge coupling g 
reads 

1 6 

Cqcd = j^KvKv + E^/^+ m /)^/ > (°) 

y /=i 

where mj are the quark masses. The up, down, and strange quarks are relatively 
light, with masses m u fx 3 MeV, m d fx 6 MeV, and m s fx 120 MeV|. There 
is a clear gap between m Uj d, which are especially small, and the masses of the 
charm, bottom and top quark with m c fx 1.25 GeV, m& ~ 4.2 GeV and m t = 
174.3±5.1 GeV [72|. It makes sense to consider the chiral limit when Nf = 2 or 



Nf = 3 quarks become massless and the other quark masses are sent to infinity. 
In this limit, the Lagrangian (|1.3|) is invariant under the transformations 



5 a ip = ieTa^ , 5 a 4> = -ieipT a , (1.4) 

and 

5 a ip = ie^Tatp , 5 a 4> = ie4>T a ^ , (1.5) 

where T a (a = 0, 1, . . . ,N* — 1) are the (hermitian) generators of the flavor 



S\J(Nf) group. The symmetry in eq. (1.4) is the vector symmetry and eq. (1.5) 



represents the axial symmetry. While the vector symmetry is still present even 
if the quarks are given a mass (of the same magnitude for all flavors), the 



3 The light quark masses are current quark masses in the MS scheme at /i = 2 GeV. 
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axial symmetry holds only in the massless theory. The corresponding Noether 
currents are 

jg = $T al ^, jf = $T al ^. (1.6) 
They are conserved upon applying the classical equations of motion. 

1.1.1 Singlet Axial Anomaly 

Let us discuss first the singlet axial symmetry (with T a = 1). An important fact 
is that this symmetry exists only in the classical case. The full quantum path 
integral is not invariant under the transformations 

Sip = ie^ifi, Sip = , (1-7) 

where the flavor index is now omitted. This explicit symmetry breaking due to 
quantum effects can be presented as an operator identity involving an anomalous 
divergence, 

where F£ v = e^uapF^p- There are many ways to derive and understand this 
relation. Historically, this so-called Adler-Bell-Jackiw (ABJ) anomaly was first 
derived by purely diagrammatic methods and showed up in the anomalous tri- 



angle graph [|73|. Using the index theorem of Atiyah and Singer |74|], which 
states 

rih-n R = Q , (1.9) 

where n^R is the number of the left-handed (right-handed) zero modes and Q 
is the topological charge of the gauge field configuration, 

Q = 32^/ d 4 xF- u F- u , (1.10) 

Fujikawa showed that in the case of the singlet axial symmetry the anomaly 
can also be understood in connection with the Jacobian of the change of the 
measure in the path integral under a global chiral transformation [J75[] . 

The axial anomaly is by far not only of theoretical interest because it implies 
that QCD has no singlet axial symmetry and therefore no associated Goldstone 
boson can be found in the particle spectrum. This is the explanation for the 
large mass of the rf particle, which is 958 MeV, whereas the mass of the tt° 
is 135 MeV. The absence of a massless flavor singlet particle also provides a 
linear relation between the topological susceptibility and the quark mass. In 
contrast to this the quenched topological susceptibility does not vanish in the 
chiral limit and is equal to the topological susceptibility of pure SU(3) gauge 
theory. Finally, through the coupling to the electroweak sector of the standard 
model the axial anomaly is responsible for the decay tt° — > 27. 



10 



Chapter 1. Introduction 



1.1.2 Spontaneous Breaking of Non-singlet Chiral Symme- 
try 

Consider now the whole set of symmetries from eqs. (|0|) and ( |l,5|) . It is 
convenient to introduce 

^L,R = 7^(1^75)^, ^L,R = -^(1±75) (1.11) 



and rewrite eqs. (|l.4j) and (|l.5| ) as 

-» VL^, ^R -» WR , (1-12) 

where Vl and Vr are two different U(Nf) matrices. The singlet axial transfor- 
mations with Vl = VS = e 1 ^ are anomalous as in the theory with a single quark 
flavor. Therefore, the true fermionic symmetry group of massless QCD is 

G = SU L (JV» x SU R (JV» x U v (l) • (1.13) 



There is ample experimental evidence that the symmetry in eq. (LIS) is actually 
spontaneously broken, which means that the vacuum state is not invariant under 
the action of the group Q. The symmetry Q is, however, not broken completely. 
The vacuum is still invariant under transformations with Vl = Vr, generated 
by the vector current. 

Thus, the pattern of breaking is 

SU L (JV» x SU R (iV» ->SUv(iV». (1.14) 
The vacuum expectation values 

^ = ml) (1-15) 

are the order parameters of the spontaneously broken axial symmetry. The 
matrix is referred to as the quark condensate matrix. 

The non-breaking of the vector symmetry implies that the matrix order 
parameter from eq. ( 1.15| ) can be cast in the form 



S a/3 = -^5 af3 (1.16) 

by group transformations from eq. (|1.12| ). This means that the general con- 
densate matrix is a unitary SU(iV/-) matrix multiplied by the real constant 
S. 

By Goldstone's theorem, spontaneous breaking of a global continuous sym- 
metry leads to the appearance of purely massless Goldstone bosons. Their 
number coincides with the number of broken generators, which is Nj — 1 in 
our case. As it is the axial symmetry which is broken, the Goldstone particles 
are pseudoscalars. They are the pions for Nf = 2 or the octet (ir,K,K,r}) for 



1.2. Chiral Symmetry on the Lattice 
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Nf = 3. It is a fundamental and important fact that spontaneous breaking of 
continuous symmetries not only creates massless Goldstone particles, but also 
fixes the interactions of the latter at low energies: a fact which is also called 
soft pion theorem. 

The Goldstone particles are massless whereas all other states in the physical 
spectrum have nonzero mass. Therefore, we have two distinct energy scales and 
one can write down an effective Lagrangian depending only on slow Goldstone 
fields with the fast degrees of freedom corresponding to all other particles being 
integrated out. The corresponding effective theory is called chiral perturbation 
theory. 

The Lagrangian of real QCD from eq. (|1.3[) is not invariant under the axial 
symmetry transformations just because quarks have nonzero masses. However, 
the symmetry in eq. ( 1.1$) is still very much relevant to QCD because some 
of the quarks happen to be very light. For Nf = 2, spontaneous breaking of 
an exact SUl(2) x SUr(2) symmetry would lead to the existence of 3 strictly 
massless pions. As the symmetry is not quite exact, the pions have a small 
mass. However, their mass goes to zero in the chiral limit m u d — * 0. This 
fact is encoded in the Gell-Mann-Oakes-Renner (GMOR) relation which can be 
derived from the chiral Ward identities 

fiml = 2(m u + m d )Z + 0(m 2 q ). (1.17) 

The constant f n appears also in the matrix element 

(0\A+\7T) p = if nP l 

of the axial-vector current = d'y^^u and determines the charged pion decay 
rate. Experimentally, f w 131 MeV. 

In contrast to f n it is not possible to give an accurate number for E from 
experimental data and therefore data from lattice calculations can help to im- 
prove the estimates for the scalar condensate. Closely related is the problem of 
the overall scale of the quark masses; x^T can provide accurate predictions for 
the quark mass ratios, but the exact scale can not be determined within this 
framework. 

This shows that the predictive power of x^T relies on measurements which 
fix its low-energy constants. In certain cases lattice calculations are the only 
viable source of data. This brings us to our next topic, which is the lattice 
Dirac operator and the problems to realize chiral symmetry on the lattice. In 
this short introduction we do not cover the basic ideas of the lattice approach to 
QCD and refer the reader to the books |76-78| or the more recent review |7S]. 



1.2 Chiral Symmetry on the Lattice 

This section deals with the general properties of lattice fermions. After recalling 
the fermion doubling problem and the Nielsen-Ninomiya No-Go theorem we shall 
discuss the recent progress made in formulating chiral symmetry on the lattice. 
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1.2.1 Fermion Doubling and the Nielsen-Ninomiya No-Go 
Theorem 

Suppose we want to describe massless free fermions on the lattice. The fermionic 
part of the QCD lattice action can be written as 

Af = Y J i>(?)D(x-y)ip{y), (1.18) 
x,y 

where D denotes the lattice Dirac operator. In particular, one would like to 
formulate the theory such that D satisfies the following conditions 

(a) D(x — y) is local, i.e. the absolute values of its couplings p{x — y) are 
bounded by an exponential function exp(— v\x — y\) with v > (see also 
Section |1.2.2| ). 

(b) D{p) = i^ Pu + 0(ap 2 ) 

(c) D(p) is invertible for p ^ 

(d) 75 D + L> 75 = 

Locality is required in order to ensure renormalizability and universality of the 
continuum limit; it ensures that a consistent field theory is obtained. Further- 
more, condition (c) ensures that no additional poles occur at non-zero momen- 
tum. If this is not satisfied, as is the case for the so called "naive" discretization 
of the Dirac operator, additional poles corresponding to spurious fermion states 
can appear: this is the famous fermion doubling problem. Finally, condition (d) 
implies that Af is chirally invariant. 



The main conclusion of the Nielsen-Ninomiya No-Go theorem [ 80 1 is that 
conditions (a)-(d) cannot be satisfied simultaneously. Since one is not willing 
to give up locality and condition (b), this implies that one is usually confronted 
with the choice of tolerating either doubler states or explicit chiral symmetry 
breaking. This is manifest in the two most widely used lattice fermion formu- 
lations: staggered fermions leave a chiral U(l) subgroup invariant, but only 
partially reduce the number of doubler species ||. Wilson fermions, on the 
other hand, remove the doublers entirely at the expense of breaking chiral sym- 
metry explicitly. This is easily seen from the expression for the free Wilson Dirac 
operator 

Dw = \l, (V„ + V*) - ^V* V„. (1.19) 

Using the definitions for the forward and backward lattice derivatives, V M and 
V*, one easily proves conditions (a)-(c), while it is obvious that (d) is not 
satisfied. 

But even though the Nielsen-Ninomiya No-Go theorem is correct, there is a 
way to have chiral symmetry on the lattice; the solution is to relax the condition 
(d) in a particular way, which we discuss in the following. 
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1.2.2 The Ginsparg- Wilson Relation and some Consequences 

As we already mentioned, formulations of chiral fermions on the lattice have 
been found that result from completely different constructions. However, they 
are all connected by the fact that they satisfy the GW relation in eq. Ql.lj) , 

Before we start the discussion about the GW relation, let us introduce some 
notations. We denote a Dirac operator that satisfies the GW relation, as it 
stands in eq. (|1 . 1[) , i.e. with a general operator R, by £> and, in particular, ©pp 
for the FP Dirac operator. In the special case when R = 1/2 we use V to 
denote the Dirac operator and the GW relation reads 

{V, l5 } = V l6 V, (1.20) 

or equivalently 

{P- 1 , 75 } = 75 . (1.21) 

From this equation one sees that the Dirac operator V is normal, i.e. it com- 
mutes with its hermitian conjugate and therefore the eigenstates of V are si- 
multaneously eigenstates of . The GW relation in combination with the 
75-hermiticity of the Dirac operator implies furthermore that V and 75 com- 
mute in the subspace of the real eigenmodes of V. Hence, real eigenmodes ip 
of V, in particular zero modes, have definite chirality, i.e. 75^ = ±ip. Another 
consequence is that the eigenvalue spectrum of V lies on a circle in the complex 
plane with radius 1 and the origin at (1,0). Note however that this is not the 
case for a general R. In this case the eigenvalue spectrum is bounded by 2 
circles, one with radius A~ in /2, the other with radius A~^ !C /2, where A m j n and 
Amax are the smallest and the largest eigenvalue of R, respectively^. In all cases 



the circles touch the imaginary axis at the origin, as shown in Figure [LI 
The Dirac operator -D can always be rescaled as follows 

{2R) 1 / 2 t){2R) 1 / 2 = V, (1.22) 

which is often very convenient to derive relations in a simpler way and therefore 
we will use V in our following discussion; even though all the results - with 
the obvious modifications - are also valid for &). The rescaled FP solution is 
denoted by £>fp- 

In order to see that the GW relation is indeed equivalent to chiral symmetry 
with a non-zero lattice spacing, we perform an infinitesimal change of variables 
ip — ► ip + ie5xp and tp — > tp + iedtp with a global flavor singlet transformation pi 



Sr ,, (l-ip)r (1.23) 

,)r r[ l-ip)-,-, 



4 The operator R is local and hermitian and its eigenvalues are real and bounded by a 
constant A > 0. 
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Figure 1.1: The eigenvalue spectrum of GW Dirac operators in the complex 
plane. In the case R oc 1 the spectrum lies exactly on a circle (left), whereas 
in the general case the spectrum is bounded by two circles and lies in the 
shaded area (right). 



From the GW relation eq. ( |1.20| ) it follows that the action Af = i/jVi/j is 
invariant under this transformation, i.e. 

5$ViP)=0. (1.24) 

The corresponding flavor non-singlet transformation is given by 

S a ip = T al5 (l-~V\T/;, (1.25) 



5 a ^ = ^[1- -V) l5 T a . 



Analogous to Fujikawa's observation discussed in Section 1.1 the fermionic 
integration measure is not invariant under the singlet transformation in eq. ( 1.23 ) 



= Tr( 75 V)[dil)di/>] 

= 2Nt{n h - nn)[d$dip] , 



(1.26) 



where the trace Tr denotes the trace over all indices, i.e. color, flavor, Dirac and 



spatial indices. To derive eq. ( |l.26D we make use of the lattice version of the 



Atiyah-Singer index theorem [42]. Hence, the measure breaks the flavor singlet 
chiral symmetry in a topologically non-trivial gauge field where the index(P) = 
nh — nR 7^ 0. Before we discuss the index theorem in more detail, let us mention 
that the non-singlet symmetry is not anomalous, since eq. ( 1.26| ) for this case 
contains Tr^PT) oc tr(T) = 0. 
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Having a Dirac operator satisfying the GW relation one can easily show that 
the following identity, the lattice index theorem, holds 

index(P) = ^TV(7 5 P). (1.27) 



Furthermore, one can derive 



\ti(j 5 V) (1.28) 



ir 



n the continuum limit [47, 4S[. This shows that the index of the lattice Dirac 



operator is indeed connected to the topological charge or winding number of 



the gauge fields (see also eqs. (|1.9|) and ( [LIOD ). Moreover, the index theorem 
provides a new definition for the topological charge density on the lattice, namely 

q(x) = ^tr( 75 V(x,x)), (1.29) 

where the trace is over Dirac and color indices. 

It has been shown that GW Dirac operators can not be ultra-local [81-84], 



i.e. their couplings p(x — y) can not vanish for all \x — y\ > r for any finite 
r. This fact is more a technical than a fundamental problem in simulations, 
because the couplings of the Dirac operator - at least for any acceptable solution 
of the GW relation, such as the FP Dirac operator - fall off exponentially. Such 
Dirac operators are still local in a physical sense, as the localization range of 
the couplings measured in physical units shrinks to in the continuum limit. 
However, the fact that no ultra-local solutions exist already gives a hint that 
it is not easy to find an approximation to a GW fermion that can be used in 
numerical simulations. But let us discuss now one of the solutions of the GW 
relation, the FP Dirac operator and the FP approach to QCD, more in detail. 



1.3 The Renormalization Group and Fixed Point 
Actions 

Field Theories are defined over a large range of scales. However, not all the 
scales are equally important for the underlying physics. Whereas the low energy 
excitations carry the information on the long range properties of the theory, the 
degrees of freedom associated with very high (and unphysical) scales do influence 
the physical predictions only indirectly through a complex cascade process. This 
makes the connection between the local form of the interaction described by the 
Lagrangian and the final physical predictions obscure and, moreover, introduces 
technical difficulties in treating the large number of degrees of freedom. That 
is the reason why one attempts to separate the low energy scales from the 
scales associated with very high energies by integrating them out in the path 
integral. The method accomplishing this, taking into account the effect on the 
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remaining variables exactly, is called the Renormalization Group transformation 

(RGT) nm-m. 

The repeated use of the RGT on an initial theory gives a sequence of theories. 
These theories can be described by a flow trajectory in the space of couplings. 
Fixed points (FP) of this transformation are theories that reproduce itself under 
the RGT. Since the correlation length £ of the theory scales by the scaling factor 
of the RGT, its value has to be or oo at the fixed point. Yang-Mills theory has 
a non-trivial fixed point (the Gaussian FP) with correlation length £ = oo whose 
exact location in coupling space c\, C2, . . . depends on the RGT used. There is 
one so-called relevant coupling whose strength increases as one is starting near 
this FP and performing RGTs. In QCD the quark masses are additional relevant 
couplings of the RGT. The flow along these relevant scaling fields whose end- 
point is the FP is called the renormalized trajectory (RT). Simulations performed 
using an action which is on the exact RT would reproduce continuum physics 
without any discretization errors. Let us finally mention that a detailed review 
on FP actions can be found in p||. 



1.3.1 Fixed Point Action of QCD 

We consider the QCD Lagrangian which consists of a SU(3) non-abelian gauge 
theory^and a fermionic part describing the quarks in interaction with the gauge 
fields in four dimensional Euclidean space defined on a periodic lattice. The 
partition function is defined through 

Z(fi) = J DUD^D^exp[-PAg(U)-Af(U,^,ip)], (1.30) 

where DU is the invariant group measure, DtpDtp the (anti)fermion integration 
measure and j3 the gauge coupling. The lattice actions A g {U) and Af(U,^,ip) 
are some lattice regularizations of the corresponding continuum gauge and 
fermion action, respectively. We can perform a real space renormalization group 
transformation (RGT), 

eM-W'A' g (V)+A' f (V,x,x))] = 
J DUD^Dip e W [-{(3(A g (U)+T g (U, V))+(A f (U, i>)+T f (U, X , X, i/>))}], 

(1.31) 

where V is the blocked link variable, x ar, d X the blocked fermion fields. Finally, 
T g (U,V) and Ty([7, x, x, VO are the blocking kernels defining the RGT, 

T g (U,V) = -|E (Retr(^(n B )Qt(n B )) - (1.32) 



5 The discussion is kept in the physically relevant case where the number of colors N c 
is 3. The discussion might, however, easily be translated to an arbitrary number of colors 
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and 

Tf(U, 



us n 

$ n ). (1.33) 



In eq. ( |1.32|) Q M (ne) is a 3 x 3 matrix representing some mean of products of 
link variables U^{n), connecting the sites 2rag and 2(ng + /}) on the fine lattice 
and Afjx is a normalization constant ensuring the invariance of the partition 
function. By optimizing the averaging function in and the parameter K g , it 
is possible to obtain an action on the coarse lattice, which has a short interaction 
range. Such an optimization has been done and we refer to [^|] for the explicit 
form of the RGT block transformation. 

The fermionic blocking kernel in eq. ( |l .33 ) is defined by the averaging func- 



tion Q,{U) nB)n , which connects the fermion fields on the fine lattice ip to the 
fermion fields on the coarse lattice \- It ' s a generic gauge invariant func- 
tion with hypercubic symmetry that on trivial gauge configurations satisfies the 
normalization condition Y^ n ^(^ n B — n) = 2 _ 2~~ for d space-time dimensions. 

On the critical surface at j3 — > oo eq. ( |1.3ip is dominated by the gauge part 
and it becomes a saddle point problem representing an implicit equation for the 
FP action, A FP , 

A FP (V) = min {A F g p (U) + T(U, V)} . (1.34) 

The normalization constant in the blocking kernel, J\f^, becomes in the limit 

(3 — > oo 

= ^ {Retr(^Q/,t)} . (1.35) 

Note that the FP eq. ( |1.34| ) defines the FP gauge action as well as the field on 
the fine configuration U[V] implicitly through a minimization condition. 

As the gauge part completely dominates the path integral for large values 
of the gauge coupling (3, the fermionic part does not have any influence on the 
solution of the FP gauge action. This is, however, completely different for the 
fermionic FP action as the gauge fields on the fine as well on the coarse lattice 
show up in the fermionic FP equation [p4| 



®£(Y) = — + n(u[v])^(u[v])n\u[v}) (l.se) 
/ 

or equivalently 

© FP (y) = Kf - K 2 f n(u[v])[® FP (u[v}) + K / nt(i7[y])n(c/[F])]- 1 nt( [ /[y]) . 

(1.37) 
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Combining eq. ( 1.36| ) with the GW relation from eq. (|1.2|) , we obtain the fol- 
lowing FP equation for the operator R appearing in the GW relationp] 



r fp (v) = — + n(u[v])RMu[v])rt(u[v]) 



(1.38) 



The FP eqs. ( fL3g) , ( |i~3p , ( |L37D and ( |TT38t ) , respectively, can be studied 
analytically up to quadratic order in the vector potentials []93|, |94|. However, 
for solving the FP equations on coarse configurations with large fluctuations - 
as they are needed for numerical simulations in QCD - one has to resort to 
numerical methods, and a sufficiently rich parametrization for the description 



of the solution is required. The numerical solution of the eq. ( |L34|) defining 
the FP gauge action is described in [95,^], whereas a substantial part of this 



work is devoted to the numerical solution of the fermionic FP eqs. ( 1.361 ) and 
1,37 ) 7 [. But before we give an account of the parametrization problem, we 



first discuss the properties of a general Dirac operator, because we need a Dirac 
operator with a rich structure to be able to capture some of the very important 
properties of the FP Dirac operator. 



Notice that this d efinit ion o f R is very specific for the FP approach. 
7 Even though eqs. ( 1.36) and (1.37) are equivalent their influence in the parametrization 
shows to be rather different, as we will discuss in Chapter^ 



Chapter 2 



General Lattice Dirac 
Operator 

In this chapter we discuss the steps to construct Dirac operators which have 
arbitrary fermion offsets, gauge paths, a general structure in Dirac space and 
satisfy the basic symmetries on the lattice, which are gauge symmetry, hermitic- 
ity condition, charge conjugation, hypercubic rotations and reflections. We give 
an extensive set of examples and provide an efficient factorization of the opera- 
tors occurring in the construction of a Dirac operator. Although the discussion 
about the lattice Dirac operators is kept very general, all the examples are spe- 
cific to the construction of the parametrized FP Dirac operator £> pa r, which we 
discuss in Chapter |3[ 

2.1 Introduction 

The construction of a general lattice Dirac operator that satisfies all the basic 
symmetries (gauge symmetry, 75 hermiticity, charge conjugation, hypercubic 
rotations and reflections) is a fundamental kinematic problem. We discuss this 
problem in d = 4 in a very general way. Notice that a similar discussion with 
different notation can be found in |22|,|24]]. For any fermion offset and for any of 
the 16 elements of the Clifford algebra we describe the steps to find combinations 
of gauge paths which satisfy all the basic symmetries. We construct explicitly 
paths for all the elements of the Clifford algebra in the offsets of the hypercube. 
We show also how to factorize the sum of paths, i.e. writing it as a product of 
sums, in such a way that the computational problem is manageable even if the 
number of paths is large. 

Most of the contents of this chapter are published in |T5]] as well as in 
hep-lat/0003013, where an easy-to-use Maple code dirac .maple is provided. 
This code can be used, if one wants to add additional terms to a parametrization 
of a Dirac operator. 

Let us add a few remarks at this point: 
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There are many ways to fix the parameters of the chosen Ansatz for the 
Dirac operator. As opposed to production runs this problem should be 
treated only once. It is useful to invest effort here, since the choice will 
influence strongly the quality of the results in simulations. 

There are compelling reasons to use the elements of the Clifford algebra 
beyond 1 and 7 M in the Dirac operator. For an operator © satisfying 
the GW relation in eq. (|l7l|), tr^i?©) is the topological charge density 
[21, 42 1 , i.e. the 75 part of £> is obviously important. Similarly, a term 
is already required by the leading Symanzik condition. 

If the parametrization is close to ®fp> the GW relation is approximately 
satisfied, which means the operator JvA under the square root in Neu- 
berger's overlap construction is close to 1. An expansion converges very 
fast in this case, as we will show in Chapter ||. In the present situation 
it is, however, unclear whether this speed-up of the convergence compen- 
sates for additional expenses of the improved operator [33,97,98|- But, 
there are many reasons to use improved Dirac operators and performance 
issues in the overlap construction are most likely not more important than 
good scaling properties of the resulting operator. 

The basic numerical operation in production runs is ® par ^, where £> par is 
the parameterized Dirac matrix and r\ is a vector. The matrix elements 
of ®par should be precalculated before the iteration starts. Using all the 
Clifford algebra elements and arbitrary gauge paths, the computational 
cost per offset of the operation © par 77 is a factor of ~ 4 higher than that 
of the Wilson action. Using all the points of the hypercube (81 offsets), 
the cost per iteration is increased by a factor of 4 x 81/8 ~ 40 relative 
to the Wilson Dirac operator. Note however, that this number may vary 
quite a bit depending on the underlying computer architecture, as one of 
the main issues in present lattice simulations is rather fast memory access 
and fast communication than fast floating point units p9|]. 



2.2 Symmetries of the Dirac Operator 

We define the basis of the Clifford algebra as T = 1, 7^, ia^ u , 75, 7^75. We 
use the notation S,V,T,P and A for the scalar, vector, tensor, pseudoscalar and 
axial-vector elements of the Clifford algebra, respectively. Notice that the tensor 
(T) and axial-vector (A) basis elements of the Clifford algebra are anti-hermitian 
as one can see from the explicit representation given in Appendix [A|. It will later 
be useful to list the basis elements of the Clifford algebra by a single index, as 
Tj, i = 1, . . . , 16. We choose the ordering 



Ti = 1 , 71 , 72 , 73 , 74 , io-12 , «cri3 , icru , i&23 , ^24 , i&34 , 75 , 7i 75 , 7275 , 7375 , 7475 • 

(2.1) 
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Under charge conjugation, the Dirac matrices 7^ transform as 

C 7 JC- 1 = - 7 „. ( 2 - 2 ) 

Having chosen the T and A basis elements of the Clifford algebra to be 
anti-hermitian, all the basis elements have the property that 

c 75 r* 7 5C- 1 = r . (2.3) 

We define a sign er by the relation 

7sr f 75 = err, (2.4) 

which gives es = ep = £a = 1 and ey = ex = — 1. 

Under reflection of the coordinate axis 77, the basis elements of the Clif- 
ford algebra are transformed as V = V^TV' 1 , where V v = 7^75- (When 
T is written in terms of Dirac matrices 7^ this amounts simply to reversing 
the sign of 7^.) The group of reflections has 2 4 = 16 elements. Under a 
permutation (1,2,3,4) — ► {pi,P2,P3,P4) of the coordinate axes the basis ele- 
ments of the Clifford algebra transform by replacing fj, — > p^ and accordingly 

75 — * e piP2P3P4^5- 

The matrix elements of the Dirac operator D we denote as 

D(n,n';U)£, (2.5) 

where n, a and a refer to the coordinate, color and Dirac indices, respectively. 
From now on, we suppress the color and Dirac indices. The gauge configuration 
U appearing in D is not necessarily the original one entering the gauge action 
- it could be a smeared configuration V. Using the Dirac operator D(V) 
should be interpreted as defining a new Dirac operator D' on the thermally 
generated configuration U: D'{U) = D(V), where D' has a large number of 
additional paths relative to D. Provided the smearing has appropriate symmetry 
properties^ all the following constructions in this chapter remain valid, without 
any extra modifications. 

The Dirac operator should satisfy the following symmetry requirements: 

Gauge symmetry 

Under the gauge transformation U^{n) — > U^(n) 9 = g(n) U ^(n) g(n + fry, 
where g(n) G SU(iV) we have 

D(n, n'\ U) -» D{n, n; U 9 ) = g(n) D(n, n'; U) g{n'y*. (2.6) 



X A11 conventional smearing schemes have appropriate symmetry properties, in particu- 
lar also the two smearing schemes used in the parametrization of the FP Dirac operator. 
A discussion of these smearing schemes can be found in Section 'iA. 
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Translation invar iance 

Translation symmetry requires that D(n,n + r) depends on n only through the 
n-dependence of the gauge fields. There is no explicit n-dependence beyond 
that. In particular, the coefficients in front of the different paths which enter 
D do not depend on n. 

Hermiticity 

D(n,n';U) = l5 D(n',n,U)^ 5 , (2.7) 
where f is hermitian conjugation in color and Dirac space. 

Charge conjugation 

D(n,n';U) =CD(n',n;U*) T C- 1 , (2.8) 

where T is the transpose operation in color and Dirac space. 
It will be useful to combine eqs. ( |2.7| ) and ( f2~g| ) to obtain 

D(n, n'; U) = C 75 D(n, n'; U*)* ^C' 1 . (2.9) 
Reflection of the coordinate axis 7] 

D(n,n';U) = V' 1 D(n,n';U V ^)V v , (2.10) 

where V v = 7^75 and h v = n v if v / rj, while h v = —n v . The reflected gauge 
field U Vr > is defined as 

<"(m) = U n (m-fj^, (2.11) 
U u v (m) = U v {fh) , v^i]- 

Permutation of the coordinate axes 

These are defined in a straightforward way, by permuting the Lorentz indices 
appearing in D. Note that rotations by 90° on a hypercubic lattice can be 
replaced by reflections and permutations of the coordinate axes. 

2.3 Constructing a Matrix D which satisfies the 
Symmetries 

To describe a general Dirac operator in compact notations, it is convenient to 
introduce the operator of the parallel transport for direction \x 

(%) =^(n)<WAX> ( 2 - 12 ) 

V / nn' 
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and analogously for the opposite direction 

(U-v) =U^n-p)U n ^ nl . (2.13) 

V / nn' 

Obviously (u^\ = U-JQ The Wilson Dirac operator at bare mass equal to 
zero reads: 

It is also useful to introduce the operator U(l) of the parallel transport along 
some path I = [Zi , I2, • • • , h] where Zj = ±1, . . . , ±4, by 

U{l) = U h U h ...U h . (2.15) 

In terms of gauge links this is 

(U(l)) =(u h (n)U h (n + h)...)s n+rhnl , (2.16) 

V / nn' V / 

where 77 = l\ + . . . + Z& is the offset corresponding to the path /. Note that 
U(l)i = U(l) where I = [-l k , ... , is the inverse path. In particular, one 
has 

(#([2,1,-2])) =C/ 2 (n)C/i(n + 2)C/ 2 (n + i)t5 n+in , (2.17) 

V / nn' ' 

for the corresponding staple. 

As another example, the Sheikholeslami-Wohlert (or clover) term |14] intro- 
duced to cancel the 0{a) artifacts is given up to a constant prefactor by 

icr^u (tr(\p,v,-fi,-v]) +U([i/,-fj,,-v,fi])+ 

U([-fi, -v, ^ v)) + U([-u, /x, v, -ft]) - h.c.) . (2.18) 
We consider a general form of the Dirac operator 

D = ^r i ^cf ) i7(0. (2.19) 

i I 

The Dirac indices are carried by Tj = 1, 7^, . . . , the coordinate and color indices 
by the operators U(l). The Dirac operator is determined by the set of paths I 

(i) 

the sum runs over, and the coefficients c, . 



In the case of D w (see eq. (|2.14|) ), for T = 1 one has I = [1], [-1], . . . , [-4] 



and I = [] (the empty path I = [] corresponds to U([]) = 1), while for T = 

2 Note that in terms of these operators the forward and backward covariant derivatives 
are: 9 M = ?7 M — 1, = 1 — U^, and c£9 M = U IJ , + U} l — 2. 
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7 M : I = [fj] and [—fi\. In the clover term, eq. (2.18) for F = ia^: I 



[/j,, v, — /j,, — u], . . . (altogether 6 x 4 = 24 plaquette products). As these well 

(i) 

known examples indicate, the coefficients e, for related paths differ only in 
relative signs, which are fixed by symmetry requirements. 

Our aim is to give for all F's and offsets r on the hypercube a set of paths 
and to determine the relative sign for paths related to each other by symmetry 
transformations. We give the general rules for arbitrary offsets and paths as 
well. 

Eqs. ( ^3| , |2.9| ) imply that the coefficients of in eq. ( 2.19| ) are real. Further, 



from hermiticity in this language it follows that the path I and the opposite 
path I (or equivalently, U(l) and U(l)') should enter in the combination 

r(^(0+er^(O f ) , (2-20) 

where the sign er is defined by 75^75 = erT, eq. (|2.4| ). 

The symmetry transformations formulated in terms of matrix elements in 
the previous section can be translated to the formalism used here. The reflec- 
tions and permutations act on operators U(l) in a straightforward way. Under 
a reflection of the axis r\ one has U(l) — ► U(l') where l[ = —k if \k\ = 77 and 
unchanged otherwise. Under a permutation (pi,P2,P3,P4) a component with 
li = ±/Lt is replaced by Zj = ±p^, as expected. The number of combined symme- 
try transformations is 16 x 24 = 384. We denote the action of a transformation 
a = 1,... ,384 by F V^ a \ U(l) -> U(l^). Acti ng on the expression in 
eq. (|2.20| ) by all 384 elements of the symmetry group and adding the resulting 
operators together, the sum will satisfy the required symmetry conditions for a 
Dirac operator. 

Let us introduce the notation 

d(T, l) = ^J2 r(Q) (^( /(Q) ) + er#(J (a) ) + ) • (2.21) 

a 

A general Dirac operator will be a linear combination of such terms, unless 
one chooses the coefficients to be gauge invariant functions of the gauge fields, 
which can take different values for different offsets generated by the reference 
path. We will be more specific about this possibility to extend the construction 
of a general Dirac operator using color singlet factors. The normalization factor 
M will be defined below. The total number of terms in eq. ( |2.21| ) is 2 x 384 = 
768. Typically, however, the number of different terms which survive after the 
summation is much smaller. It can happen that for a choice of starting V and 



I the sum in eq. (|2.21|) is zero. In this case the given path does not contribute 
to the Dirac structure F. 

To fix the convention for the overall sign we single out a definite term in the 



sum of eq. ( 2.21| ) and take its sign to be +1. Denote by Fq, Iq the corresponding 



quantities of this reference term, and by ro = r(Zrj) the offset of Iq. This term is 
specified by narrowing down the set {F^J^} to a single member as follows: 
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a) Given an offset r = r(l) = (n, r2, ^3^4) the reflections and permutations 
create all offsets (±r Pl , ±r P2 , ±r P3 , ±r P4 ) where (pi,J>2,P3,P4) is an 
arbitrary permutation. We choose for the reference offset ro the one from 
the set {r(l^)} which satisfies the relations roi > ro2 > ro3 > ro4 > 0. 

b) If several T matrices are generated to this offset then choose as To the 
one which comes first in the natural order, eq. (|2.1|). 



c) Consider all the paths {l^} having offset ro and associated with To, 
i.e. r(Z( Q )) = ro and T^ = Tq. To single out one path Iq from this set, 
we associate to a path I = [h,l2, • • • ,h] a decimal code d\d 2 ■ ■ ■ c4 with 
digits di = li if li > and di = 9 + k for ij < 0. The path with the 
smallest code will be the reference path Iq. (In other words we take the 
first in lexical order defined by the ordering 1, 2, 3, 4, —4, —3, —2, —1.) 

Of course, one can take To, Iq as the starting T and I, and we shall refer 
to the expression in eq. ( |2.21| ) as (1(Tq,Iq) to indicate that it is associated to a 



class rather than to a specific (T,l). 

We turn now to the normalization of d(To,Zo)- l n general, there will be 
K different paths in the set {l^ \ T^ = ro,r(/( Q )) = ro}, i.e. corresponding 
to the same offset ro and Dirac structure 1~V The normalization is fixed by 
requiring that the coefficient of the reference term TqU(Iq) is +1/K. 

Consider a simple example explicitly. Let ro = (1,0,0,0), Tq = io\2 
and Iq = [2,1,-2]. The starting term, in eq. ( pT20| ) is ia 12 (U ([2,1, -2}) - 



U([2, — 1, — 2])). Applying all the 16 different reflections gives 

8ia 12 (U([2, 1, -2]) - U([2, -1, -2]) - U{{-2, 1, 2]) + U([-2, -1, 2]) 



(2.22) 

Applying all the permutations on this expression results in: 

d(T , Jo) = ^ {l6m 12 (U{[2, 1, -2]) - U{[-2, 1, 2])) + 
16m 13 (U([3, 1, -3]) - U([-3, 1, 3])) + 

lGiau (U([A, 1, -4]) - U([-4, 1,4])) +...}, (2.23) 

where only the terms with the offset r = ro are written out explicitly. Their 

total number is 96. The whole generated set has 8 different offsets giving 

altogether 768 terms. Notice the form of the contribution in eq. ( [2. 23 ). There 



is a common factor (16 in this case) multiplying all the different operators. 
Only the tensor elements of the Clifford algebra enter, since we started with a 
tensor element. Beyond the common factor the path products have a coefficient 
±1. These features are general. The number of different paths with Tq = 

ia 12 and r = (1,0,0,0) is K=2: the paths [2,1,-2] and [-2,1,2]). The 
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normalization factor in this case is J\f = 16 x 2 = 32, so that one has d(To, Iq) = 
\ia l2 (#([2, 1,-2]) -#([-2,1, 2])) + .... 
A general Dirac operator is constructed as 

D = £ jf.Y,fi( l{a) ) r i a) (W^ + erW^) ■ (2-24) 

i=F ,l 1 a 

The coefficients fi(l^) which are the free, adjustable parameters of the Dirac 
operator are real constants or more generally gauge invariant, real functions of 
the gauge fields, respecting locality, and invariance under the symmetry trans- 
formations. The requirement that such a function is gauge invariant and real, 
restricts the choice to the real part of traces of closed loops of link products. 
The invariance under symmetry transformations can be assured by constructing 
these functions along the same lines as a scalar operator which satisfies the sym- 
metry requirements. The following example should make clear what is meant by 
such a color singlet function. We consider again the coupling generated by the 
staples and the tensor Clifford algebra elements i.e. vq = (1,0,0,0), Tq = iau 
and l = [2,1,-2]. A simple choice for such a function which satisfies all the 
requirements would e.g. be a polynomial of the trace of all the plaquettes gen- 
erated by the operator r = (0,0,0,0), T Q = 1 and l = [1,2,-1,-2]. But 
there are other possible choices, such as 

/(/«) = /(x(Z«)) 

with 

x(Z«) = l-i-tr[X(Z«)X(^))t], (2.25) 

where X is the gauge operator^ generated by tq = (1,0,0,0), To = 1 and 
lo = [2, 1, —2], N c the number of colors and iw, i = 1, . . . ,n ff sets denotes the 
different offsets generated by Iq. This choice gives the possibility that different 
offsets generated from the same To and lo can have different couplings in the 
sense that the value of the function f(l^) is not necessarily the same for every 
offset^. In the parametrization of ©pp ar| d R we always use a linear function 

f(x(lV))=c + c 1 x(lV) (2.26) 

for the couplings. 

3 The Dirac algebra structure is given only to indicate the transformation properties, 
i.e. the operator X is defined without the trivial Dirac structure. 

4 Neuberger's construction does not generate such terms when one starts with the Wil- 
son Dirac operator. 
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2.4 Tables for Offsets on the Hypercube 

Choosing offsets and paths to be included in the Dirac operator is a matter 
of intuition. It is also influenced by considerations on CPU time and memory 
requirements. In Tables 2.1-21] we give the reference paths Iq for offsets on the 



hypercube and general Dirac structure that are used in the current implementa- 
tion of the parametrized FP Dirac operator S par - The first 3 columns give Tq, 
Iq and the number of paths K as defined in the previous section. The 4th col- 
umn gives those Clifford basis elements which are generated in eq. ( p.21| ) to the 



offset rrj . In Table ^6| we give the gauge invariant, real functions of the gauge 



fields which are used as fluctuation polynomials in the current implementation 



Of ©par- 



It is also of interest how the given Dirac operator behaves for smooth gauge 
fields, i.e. to obtain the leading terms in the formal continuum limit. The 



expression in eq. (|2.21| ) generated by given Tq and Iq contributes in this limit 



to one of the expressions (of type S,V,T,P,A) below 

(c s + c s d 2 ); cvT/A,; c^a^F^ ; 

1 

Cp-^fivpaF^uFpa] CAiliTI^ iivpudvFpv . (2.27) 

Here is the covariant derivative in the continuum. The coefficients cx(Tq, Iq), 
which determine the continuum behavior of the Dirac operator, are presented 
in the last column of Tables pTl]-|2 .51. For Tq = 1 the first and second entries 



correspond to cs and cs, respectively. We list their meaning below. 
Introduce the notation 



Cx = J2f( T o,h)cx(T ,l ), (2.28) 

where f(T ,l ) is the formal continuum limit of fi(l^) defined in eq. ( |2.24| ). 
The bare mass itiq is given by 

C s = m . (2.29) 
The normalization condition on the D ~ Ipd^ term in D gives 

C v = 1 • (2.30) 
The 0{a) tree level Symanzik condition reads 

C s + C T = . (2.31) 



The coefficient Cp is interesting if the parametrization attempts to describe 
(approximately) a GW fermion. In this case it is related to the topological 
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r 


ref. path Iq 


K 


T's generated 


cx 


1 


D 

[1,2,-1,-2] 




1 

48 


1 


o o 

I— I 1—1 


7i 


[1,2,-1,-2] 


24 


y,...,7 4 





icrn 


[1,2,-1,-2] 


8 


MT12, • • • ,iCT 3 4 


-1 


75 


[1,2,-1,-2,3,4,-3, 


-4] 


384 


75 


-1/6 




[1,2,3,4,-1,-2,-3, 


-4] 


384 




-1/6 


7i75 


[1,2,-1,-2,3,4,-3, 


-4] 


192 


7i 75, •••,7475 







[2,1,-2,-1,3,4,-3, 


-4] 


192 








Table 2.1: Reference paths for different r 's for offset (0000). 



r 


ref. path Iq 


K 


T's generated 


cx 


i 


[1] 


1 


1 


8, 1 




[2, 1, -2] 


6 




8, 1 




[2,3,1,-3,-2] 


24 




8, 1 


7i 


[1] 


1 


7i 


2 




[2,1,-2] 


6 




2 


72 


[1,2,3,-2,-3] 


16 


72,73,74 





ia 12 


[2,1,-2] 


2 


i0"l2, ^13,^14 


4 


i&23 


[1,2,3,-2,-3] 


16 




-4 


75 


[2,1,-2,3,4,-3,-4] 


96 


75 


4/3 


7i 75 


[2,1,-2,3,4,-3,-4] 


96 


7i75 





7275 


[1,3,4,-3, -4] 


16 


7275,7375,7475 


-2 



Table 2.2: Reference paths for different r 's for offset (1000). 



r 


ref. path 


K 


T's generated 


cx 


i 


[1,2] 


2 


1 


24, 6 


7i 


[1,2] 


2 


7i,72 


12 


73 


[1,3,2,-3] 


8 


73,74 







[1,2] 


2 


2(712 


-2 




[1,3,2,-3] 


4 







icr 34 


[1,2,3,4, -3,-4] 


32 




-4 


75 


[1,2,3,4, -3,-4] 


32 


75 


-2 


7175 


[1,2,3,4, -3,-4] 


16 


7175,7275 


4 


7375 


[1,4,2,-4] 


8 


7375,7475 


-4 



Table 2.3: Reference paths for different IYs for offset (1100). 
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r 


ref. path Iq 


K 


T's generated 


cx 


1 


[1,2,3] 


6 


1 


32, 12 


7i 


[1,2,3] 


4 


71,72,73 


24 


74 


[1,2,4,3,-4] 


24 


74 







[1,2,3] 


4 


£C12> ^13,^23 


-8 


iau 


[1,4,2,-4,3] 


8 


ifJi4,icr 2 4,'i0'34 





75 


[1,4,2,-4,3] 


12 


75 


-8/3 


7175 


[1,4,2,-4,3] 


8 


7175,7275,7375 


8 


7475 


[1,2,3] 


6 


7475 


-4/3 



Table 2.4: Reference paths for different Tq's for offset (1110). 



r 


ref. path Iq 


K 


r's £ 


;enerated 


cx 


i 


[1,2,3,4] 


24 


i 


16, 8 


7i 


[1,2,3,4] 


12 


7 1 


...,7 4 


16 


i<?i2 


[1,2,3,4] 


8 




• • • ,«0"34 


-8 


75 


[1,2,3,4] 


24 


75 


-2/3 


7175 


[1,2,3,4] 


12 


7175 


• • • , 7475 


8/3 



Table 2.5: Reference paths for different Tq's for offset (1111). 



offset 


ref. path Iq 


K 


(0000) 


[1,2,-1,-2] 


48 


(1000) 


[2,1,-2] 


6 


(1100) 


[1,2] 


2 


(1110) 


[1,2,3] 


6 


(1111) 


[1,2,3,4] 


24 



Table 2.6: Reference paths for the fluctuation polynomials as defined in 
eqs. fl2.24j ) and ( |2.25| ). Note however, that the fluctuation polynomial for the 
offset (1000) in the current implementation of the operator R differs from the 
one given in this table. It is given by the trace of the sum of the plaquettes 
along a particular link U(l®) i.e. x(l®) = 1 - l/iV c Re tr [X(l®)U(l®)i], 
where _/V c is the number of colors and the operator X is given by the operator 
for the offset (1000) in this table. 
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charge density tr (75 £> i2) (see Section [T.2.2| ) , 

Y] f(To,lo)cp(To,l )Rfr ee (r ) = e 5 -^- 



(2.32) 



Of course, in the expressions above To should be of the corresponding type (S,T, 
. . . ). The operator i?f ree (ro), whose couplings of are specified in Appendix 0, 
is given by R(n, n + ro; U = 1) of the GW relation eq.( |l.l|) and €5 is defined 
through 

75 = £571727374 , (2.33) 

which means that €5 = —1 for our choice of the Clifford algebra elements 
(see Appendix^). The relation ( 2.32 ) can be derived easily from the work of 
Fujikawa in pTfl T In the parametrization of $Dpp we make use of the relations 
( |2.29| )-( |2~32" ) to constrain some of the parameters in the formal continuum limit 



(see Section 3.2.2) 



2.5 Factorizing the Paths 

As the Tables |2.1[ - ^ show the number of paths, in particularly for the offset 
ro = (0,0,0,0), is large. It is important to calculate the path products effi- 
ciently. An obvious method is to factorize the paths, i.e. to write the sum of a 
large number of paths as a product of sums over shorter paths. We shall try to 
factorize in such a way that these shorter paths are mainly plaquette, or staple 
products. 

Because this factorization is a very technical issue, we defer the detailed 
discussion of the factorization of all the paths involved in the construction of 
© par and R to Appendix |B|. We however explain the idea on a simple example. 
Consider the contribution generated by ro = (0,0,0,0), To = 75 and Iq = 
[1, 2, —1, —2, 3, 4, —3, —4]. It contains together with its hermitian conjugate 
768 different paths of length 8. The sum over all these terms can, however, be 
written as 

fiupa 

with 

p(as) _ p( ) _ k 

p( ) — p _P _p _L p 

where P^ v is the plaquette in the ^-direction, which is a path of length 4. 
Hence, the whole contribution can be written as the sum of 24 different prod- 
ucts of plaquette combinations. The reduction in computational expense in 
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comparison to the calculation of all the 768 length 8 paths is obvious. The 

fas) 

observation that the plaquette combinations P^ v ', as well as other combina- 
tions of plaquettes and staples, appear also in other path factorizations makes it 
favourable to precalculate them on the whole lattice, before one starts to build 
up the Dirac operator, thereby making the concept of factorizing the paths even 
more attractive. 

Using the factorization of paths the building of the matrix ® par with all 
its 41 different terms and linear fluctuation polynomials in Tables (|2.1| )-( |2~lf ) 
is as expensive as O(20) matrix vector multiplications with ® par on PC's and 
Alpha workstations. On supercomputers, like the Hitachi SR8000 in Munich, 
the cost of building the operator is no longer as favourable, since such machines 
can optimize the matrix vector multiplication much more efficiently and there- 
fore the cost of the building of © par is as expensive as 0(500) matrix vector 
multiplications. Such a cost becomes a noticeable fraction (10 — 20%) of the 
calculation of a full quark propagator. This shows that the factorization of the 
paths is crucial in order to have a Dirac operator whose build-up time is in a 
range where it is computationally still reasonable. 



Chapter 3 



Parametrization of the FP 
Dirac Operator 



The construction of a FP Dirac operator £>fp is a non-trivial task, since it is 
defined through a highly non-linear Renormalization Group equation. While the 
solution can be calculated analytically in the free case [^4|], this is clearly not 
possible for the interacting case. Therefore, one has to find a way to approximate 
the solution as well as possible, always keeping in mind that the resulting Dirac 
operator shouldn't be too expensive for practical use. In the following, we 
describe the details of our construction of an approximation to ©fp- We choose 
to approximate £>fp by a general parametrization of a hypercubic Dirac operator 
as described in Chapter ||. It is clear that such a hypercubic parametrization, 
which we denote by £> par , is quite a drastic reduction compared to ®fp, which 
has an infinite number of couplings, like any other Dirac operator satisfying the 
GW relation [j8l)-84|. Results from e.g. the 0(3) cr-model, however, support the 



hope that a compact parametrization which encodes the main features of the 
full FP operator can be found fllPCR . The use of a parametrization that contains 
the full Clifford algebra is a very natural choice for the parametrization of £>fp, 
since all the Clifford elements are generated through the RGT in eq. fll-37 ); a 



fact which was already noticed in [94, 101 1. Furthermore, important properties 



of Dirac operators satisfying the GW relation, like that the topological charge 
density is given by trf/ys-RSD), imply that the full Clifford algebra is necessary to 
achieve a parametrization that reproduces these properties accurately. 

The SU(3) gauge configurations we use for the parametrization and later 
on in all the production runs are generated exclusively with a parametrization 
of the FP gauge action. This gauge action has scaling properties which are 
much improved compared to the Wilson gauge action. A detailed analysis of 
the properties of this FP gauge action is given in [95, 102 1 and for the corre- 



sponding anisotropic FP gauge action see [gg]. Another important part of the 
parametrization is the use of a Renormalization Group inspired smearing tech- 



nique, which we describe in Section |3.4| . The smearing of the gauge fields is 
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a very helpful procedure in the definition of a Dirac operator, because it re- 
duces the unphysical short-range fluctuations inherent in the gauge fields quite 
significantly. This makes that the problems with additive mass renormalization 
and chiral symmetry breaking, which are common to most of the traditional 
fermion formulations, are already reduced substantially, without even changing 
the fermion action 1 103 1 . The effect of this RG smearing is illustrated in Figure 
3.6 , where the eigenvalue spectrum of the Wilson Dirac operator for a pair of 
unsmeared and smeared configurations is shown. The smearing clearly reduces 
the additive mass renormalization and the fluctuations in the small (real) eigen- 
modes, which are the reason for the so-called exceptional configurations. On 
such an exceptional configuration the quark propagator can not be calculated, 
because some of the real eigenmodes get so close to zero - notably at a non- 
zero bare quark mass - that the inversion algorithms break down. Finally, we 
also explain our parametrization of the operator Rfp, which is defined through 
eq. CI). 

In the following, we first explain the strategy of the parametrization and 
then discuss the ingredients of the parametrization, such as fitting procedure, 
the blocking kernel and the smearing in detail. The fundamental part of any 
parametrization of a Dirac operator is already explained in Chapter ||, namely 
the general structure of a Dirac operator on the lattice and how one can do 
a practical construction of such a general Dirac operator. The parametriza- 
tions used in the crucial steps of the parametrization procedure and the final 
parametrizations used in the production runs are specified in Appendix |C|. 



3.1 The Equations for the FP Dirac Operator 

In order to construct our approximation to £>fp. we solve the RG equations 
that define £>fp and its inverse as good as possible within the limitations of the 
chosen ansatz of the parametrization. This means that we try to get an approx- 
imate solution for the following two equations, which are actually equivalent as 
long as &f(U) has no zero mode: 

® c (y) = K fi-K 2 f n(u[v])[® f (u[v}) + Kf n t (tf[v])fi(E/'[v])]- 1 n t (E/'[v]) 

(3.1) 

and 

n-^v) = - + n(u[v})^ 1 (u[v])n\u[v]). (3.2) 

Kf 

The labels c and / stand for the coarse and fine lattice and indicate that the 
Dirac operators on both sides of the equation are not the same operators as 
long one is not in the fixed point, where the relation $D C = $Dy = ©pp holds. 
The gauge field on the fine lattice U[V] is determined from the gauge field V on 
the coarse lattice by a minimization procedure. The exact relation of the U[V] 
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and V fields is encoded in the FP equation for the gauge fields eq. ( 1.34| ). The 
fermionic blocking kernel Q,{U) and the parameter kj are discussed in detail in 



Section 3.3 



3.2 The Fitting Procedure 



In principle one would like to have an approximate solution of eqs. ( |3.1| ) and 



(3.2) which is valid for a large range of gauge couplings. This however shows 
to be a difficult thing to achieve, because the characteristic fluctuations in the 
gauge fields change very much with the gauge coupling, as one can see in Table 



3.1| , where the expectation value of the trace of the plaquette is shown. For this 
reason we adopt an iterative method for the whole parametrization procedure 
that starts in the weak coupling regime, i.e. at large values of j3 where the 
fluctuations of the gauge fields are small and finally ends at j3 as 3.0, which 
corresponds to a lattice spacing of a rj 0.15 fm. This is amongst the coarsest 
lattice spacings where typical lattice simulations with fermions are performed. 



p 


(v) 


(w) 


(u) 


100 


2.92 


2.989 


2.998 


10 


2.33 


2.90 


2.987 


5 


1.68 


2.76 


2.97 


3.4 


1.24 


2.62 


2.95 


3.0 


1.16 


2.56 


2.94 


2.7 


1.05 


2.49 


2.94 



Table 3.1: Average plaquette values (v) for the unsmeared coarse configu- 
rations V, (w) for the RG smeared coarse configurations (see Section |3.4j ) 
and (u) for the minimizing fine configurations U[V] at different values of the 
gauge coupling (5. 



Let us first have a closer look at one step in this iterative procedure which is 
also sketched in Figure 3.1. Assume that in the n th iteration, we have a Dirac 



(n) (n) 

operator S pa r which is a good approximation to ®fp- Actually, £> pa r has to 
be a good approximation only on minimized gauge configurations [/( n )[V( n )] 

which have a certain level of fluctuations, since we use £> P ai as S> y in the RG 
equations ( |Q| ) and (|3 . 2|) , respectively. After the RGT, the resulting Dirac 
operator on the coarse lattice ® c is approximated as good as possible within 
the chosen ansatz for the parametrization. This Dirac operator © P ai^ is the 
parametrization of the next level of the iterative procedure. The idea behind 
this step is that the fluctuations of the gauge fields on the minimized gauge 
configurations [/( ra )[F( n )] are much smaller than the ones on the correspond- 
ing coarse gauge configurations 7^ and thus ©pai^ is an approximation to 
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Figure 3.1: The iterative parametrization procedure used for the parametri- 
zation of ©fp- Starting with a hypercubic parametrization of the FP Dirac 
operator ©/ = ©par on a set of fine configurations [/("■) [y( n )] ; a Renormal- 
ization Group transformation (RGT) yields © c on the corresponding coarse 
configurations W n ). Subsequently, © c is approximated using again a hyper- 
cubic parametrization which results in ©par" 1 ^- This parametrization is the 
starting point for the next step of the iterative procedure. The crucial point 
is that the fluctuations in the gauge fields are increased substantially with 
each RGT. Hence, ©pa^ ^ describes ©fp on a larger range of gauge couplings 
than ©par- It is also possible to use the resulting parametrization ©par' 1 '' on 
the same fine configurations t r ( n+1 ) [\/( n+1 )] = {/("■) [y( n )] again in order to 
drive © Dar closer to the FP. 
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(n) 

£>fp on gauge fields with much larger fluctuations than £> par . At this point it 
should be noted, that it is much easier to find a good parametrization of $Vp 
on configurations with small fluctuations. Therefore, the step which leads from 
©par to ©plr^ 1 ' 1 is the crucial step in the whole parametrization procedure and 
it has to be done with great care. The next step is now to choose a coarse con- 
figuration y( n+1 ), such that the fluctuations of the corresponding minimized 
configuration Jj( n+1 ^ [y( ra+1 )] are roughly equal to the fluctuations of and 
then £> P ar fl ' ) can be used as £>/ in the RG equations (|3~l| ) and (|3^) . Starting 
from minimized configurations which are very far in the continuum limit and 
performing these steps several times, this procedure finally yields a parametriza- 
tion of -Dpp at intermediate to strong gauge coupling. This explains roughly the 
principle of our parametrization procedure. In the following, we indicate how 
the intermediate parametrization steps are done in detail. The Dirac operator 
used in the first step and the details of the subsequent steps of the iterative 
procedure will be discussed later in this chapter. 



3.2.1 The Details of the Fit 



Let us first focus on the solution of eq. (|3.l| ), because in this equation S) c and not 



its inverse enters and this is much simpler to use for parametrization purposes 



than eq. (3.2). But for computational reasons we can not afford to invert the 
full operator in eq. ( jO| ) for the lattice sizes we use in our parametrization, 
which are 10 4 for the fine i.e. 5 4 for the coarse lattice^ and therefore we use 
this operator equation acting on a set of normalized vectors. These vectors are 
chosen from the following two sets: 

• Random vectors^] with entries from a uniform random distribution in the 
interval [—1,1] for the real and imaginary part. These vectors give a 
large weight to the high-lying modes of the blocked operator, because the 
density of the eigenmodes of the Dirac operator is heavily peaked towards 
the upper end of the spectrum, as illustrated in Figure |37 



Low-lying eigenmodes of the operator Q=Q(U) &\(U)&f(U)& (U) de- 



termined with the Ritz functional method [104, 105(] . In clear contrast 
to the random vectors these vectors give a large weight to the low-lying 
mode contributions of the blocked operator. We choose the eigenmodes 
of O, because this operator is defined on the coarse lattice and therefore 



1 For weak gauge coupling, i.e. at f3 = 100 and (3 — 10 , 6 4 lattices for the fine and 3 4 
lattices for the coarse configurations are used, since the use of larger lattices does not lead 
to a substantial improvement in the parametrization. 

2 Instead of random vectors one may also take local vectors, i.e. vectors with one single 
non-zero entry. It shows that their effect on the parametrization is in fact the same as 
the one of the random vectors. The local vectors are, however, less efficient in carrying 
information about the bulk behaviour, which is the reason that we prefer the use of random 
vectors. 
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the determination of its eigenmodes is computationally not very expen- 
sive. This would not be the case, if we chose to determine the low-lying 
eigenmodes of ®/(?7) and then block them to the coarse lattice, even 
though this would be the more natural choice. 




Figure 3.2: Eigenvalue density p(|A| 2 ) of $D par ° n a 

4 4 lattice at /3 = 3.0. 



Our ansatz for ® par contains only linear couplings, i.e. it can be written in 



the following way 



par 



(3.3) 



=l,n 



where the q are the couplings and n op the number of terms, which is 82 in the 
present parametrization. In comparison to eqs. ( 2.24| )- ([2.26D we have expanded 



all the linear fluctuation polynomials in eq. (3J3) in order to make the linear 
structure obvious. 



3.2.2 Constraints 

In order to have a Dirac operator which respects properties that are highly de- 
sirable for a decent parametrization of $Dfp we add several constraints for some 
of the couplings q during the fit. The properties we enforce at various levels of 
our parametrization procedure are relations that provide correct normalizations 
in the formal continuum limit for the following quantities: 

C\. The bare quark mass is fixed to 0. 
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C2: The speed of light in the free energy-momentum relation is 
normalized to 1. 

C3: 0{a) tree level improvement is imposed. 
C4: The topological charge density tr (75/?©) is normalized to 
^P° F pv F pa> wnere in our case £ 5 = -1 due to 
the choice of 75 = — 7i7273740- 

C5: The free field limit is such that it coincides with the hypercubic 
parametrization of the free massless FP Dirac operator given 
in Appendix^. 

Note that these constraints, which are encoded by eqs. ( QjM gjgD, only affect 
the constants in the fluctuation polynomials from eq. fl2.26|), since the linear 



terms as defined in eq. ( 2,25 ) and Table 2Jj vanish in the formal continuum 
limit. 



3.2.3 Determination of the Couplings 

We fix the coefficients Cj of our parametrization by minimizing the following 
X 2 -f unction 



— E ll^cC^^-Spax^^ir + E^, (3.4) 



x 2 

i=1 .»conf 



where the sum runs over n con f different configurations with a total of n vec = 
Y\_i _ niec vectors v^'^ used in this parametrization step. The constraints 
Ci are enforced with weights A, > 0, which can be freely chosen. Typically, we 
use Xi = 10 4 such that the constraints are enforced to a high level. 

Due to the linear ansatz in eq. fl3.3|) and the fact that the constraints Ci are 
all linear in the coefficients c, this ^-minimization amounts simply to a matrix 
inversion and therefore poses no serious computational problems. We use the 
quasi minimal residual (QMR) matrix inverter to do the actual calculation [106]. 



3 In steps I-IV of the parametrization procedure (see Table 3.2) the constraint on the 
topological charge density was imposed with the wrong sign, because we initially over- 
looked that the results in ^| are obtained with a different definition of 75. There are, 
however, several reasons why we expect this error to be of minor consequence. The con- 
straint is only imposed on the constant terms of the fluctuation polynomials, i.e. that the 
linear terms can correct for the wrongly imposed constraint, and as the cross check of the 
index calculated with the parametrization IV and the corresponding overlap construction 
(see Chapter |^) shows, this is indeed the case, i.e. that the wr ong s ign in constraint C4 



is straightened out by the linear polynomials (see also Section 4.2.5). Furthermore, the 
pseudoscalar terms are of very small size compared to the scalar, vector and tensor terms 
and there fore have a limited influence on the parametrization procedure described in Sec- 



tion 3.2.1. Finally, the final as well as all the intermediate parametrizations describe a 
valid Dirac operator that satisfies all the symmetry conditions described in Chapter |^. 
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Unfortunately, it shows that in the region of stronger gauge coupling, i.e. (3 < 10, 
the described fitting procedure can no longer cope well with all the requirements, 
especially one observes an additive mass renormalization and, more importantly, 
a spread in the low-lying eigenmodes that is too large for the purposes for which 
we intend to use this Dirac operator. A solution to this problem can be found 
within the fixed point approach itself. Namely, one can give more weight to 
the important low-lying modes by including also the FP propagator relation 
from eq. (|3.2|) into the fit. We can include this additional information from the 
propagator by setting up a different x 2 -function, which is a combination of both 



RG eqs. (3.1) and Q3.2|) and the additional constraints C\ 



X 2 = — £ \\[®c(vW)-® p3X (vV)]v^\ 



n 



vec 



3 



=1 „w 



vec 



i=1 ' n conf 



a 

^VCC 



— E \\i^HV^)-^(V^)}w^ 2 + J2^xl, (3.5) 



where a > is chosen such that the two contributions in the x 2 -function are 
approximately equal at the beginning of the ^-minimization. The vectors w^ 1 '^ 



are the same type of random vectors as discussed in Section |3.2.1| . The gauge 
configurations that are used for the propagator fit have to be in the Q = 
topological sector, because (approximate) zero modes would dominate the \ 2 ~ 
function completely. As starting point for this minimization, which is now a 
highly non-linear problem through the use of the propagators, we use the co- 
efficients Cj obtained from the best fit in the corresponding linear x 2 -pr°blem 
(a = 0). The non-linear optimization is performed with a generic simulated 



annealing algorithm [107]. Because our ansatz for © par contains 82 parameters, 
the parameter space is rather large and each function evaluation with a new set 
of parameters involves the inversion of the Dirac operator on several sources, 
making this non-linear optimization computationally very expensive and there- 
fore an exhaustive search in the parameter space can not be afforded with our 
computing resources. This is however not a big problem, as the value of the 
X 2 -function typically decreases quite rapidly in the first few sweeps through the 
parameter space and afterwards only small improvements can be achieved. It 
shows that this non-linear minimization can improve the \ 2 f° r the propagator 
relation eq. ( ET3 ) quite substantially, without destroying the quality of the fit for 



the relation eq. (|3.l[) . An illustration of these facts is given in Figure |3.3|. The 
result of this non-linear minimization is a parametrized Dirac operator which 
has only a very small remaining additive mass renormalization and small fluc- 
tuations in the low-lying modes. Furthermore, it shows that the FP equations 
enforce the GW relation in terms of the propagators (|1.2|) to hold to a good 



extent even without adding this relation as an additional constraint during the 
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fit. One even finds that the x 2 f° r the propagator relation eq. 
for the GW relation are highly correlated []65|. 



and the x 2 



1.2 



1.1 - 




0.8 - 



0.6 



0.5 



Figure 3.3: Typical convergence history of the non-linear ^-minimization 
in the simulated annealing procedure. Even though the normalized total x 2 
(thick line) only decreases by 10%, the normalized x 2 f° r the propagator 
relation eq. (3.2) (lower line) decreases by 33%, while the normalized x 2 
for the relation eq. ( |3.1|) (upper line) increases only by 9%. Notice that 



the horizontal axis does not exactly correspond to the number of steps in 
simulated annealing procedure. The graph does only show those steps in 
which the total x 2 decreases. 



3.2.4 The Initial Dirac Operator and the Parametrization 
Steps 

At the first level of the iterative parametrization procedure shown in Figure 



3.1 we use a hypercubic truncation of the free massless £>fp, which was ob- 
tained analytically in [Q for the overlapping block transformation we use in the 
parametrization procedure. This hypercubic parametrization is explicitly given 
in Appendix^. This operator has properties that are already much improved 
compared to the Wilson Dirac operator, as we show in Chapter |j. A more 
detailed account on different hypercubic parametrizations of £>fp in the free 



case can be found in ]94fl. Starting with this hypercubic parametrization of 



the free massless £>fp, the complete parametrization procedure involves 5 main 



steps, which are shown in Table |3.2| . In step I this hypercubic parametrization 
of the free massless £>fp is used on minimized j3 = 100 configurations. The 
parametrization procedure yields $D par for thermal configurations at [3 = 100. 
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This operator is inserted again on minimized /3 = 100 configurations, until the 
X 2 of the fit does no longer decrease. In step II the same iteration is performed 
on the = 10 configurations, starting with the result of the = 100 parame- 
trization. Step III, which now involves the propagators in the fit, is an iteration 
on the j3 = 2.7, . . . ,3.4 configurations, starting with the result of (3 = 10 pa- 
rametrization. In step IV a new idea is used, namely a reparametrization of the 
Dirac operator on minimized configurations. This is the only step that deviates 
from the FP idea, because it involves an overlap expansion using £> par . We 
include this step, because characteristics of minimized f3 = 3.0 configurations 
are different from thermal configurations and therefore it is difficult to find a 
parametrization that has fluctuations of the low-lying eigenmodes which are 
small enough for our purposes. Having small fluctuations on the level of the 
minimized configurations is crucial, as the RGT magnifies these fluctuations 
by roughly a factor 2. More precisely, we use the same fitting procedure as 
described in Section 3.2.1 , but the vectors used in the fit are generated by an 
order 5 Legendre expansion of the overlap Dirac operator with $D par [|. Since 
£> par is already close to satisfy the GW relation, the overlap expansion is only 
a small correction to £> pa r, making the remaining breaking of chiral symmetry 
even smaller. The effect of the reparametrization is shown in Figure 3.4. Finally, 
in step V one last RGT leads to the final parametrization on j3 = 2.7, . . . ,3.4 
configurations. 



3.3 The Fermionic Blocking Kernel Tf 



In order to perform the Renormalization Group transformation in eqs. (|3.l| ) and 
( |3.2D numerically, the averaging function Q(U), which defines how the fields tp 
from the fine lattice are blocked to the fields x on tne coarse lattice, has to be 



specified. We choose the overlapping block transformation discussed in [94], 
which has the property that in the free case all the sites on the fine lattice 
contribute equally to the average on the coarse lattice. There are essentially 
two ways how one can introduce a block transformation in the presence of gauge 
fields. One can fix the gauge and then one can work with the same definition 
of the block transformation as in the free case, i.e. that one can do a simple 
averaging of the fields on the fine lattice. We find, however, a gauge invariant 
procedure more attractive. For the overlapping block transformation gauge 
invariance can be achieved by introducing a product of gauge links connecting 
the points ng on the coarse with the points n on the fine lattice, i.e. parallel 
transporting the fine fermion fields before averaging them. Using the definitions 



4 The Legendre expansion of the overlap Dirac operator with 5) par is described in detail 
in Chapter ra. 
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Figure 3.4: The effect of the overlap reparametrization on minimized [3 = 3.0 
configurations on the low-lying eigenvalues is shown (notice the horizontal 
scale) . The initial parametrization (filled circles) has a larger additive mass 
renormalization and the fluctuations in the low-lying eigenmodes are appar- 
ently larger than for the reparametrized operator (open circles). 
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step 


p 


size 


^■conf 


n v 


^conf 




Ci 


c 2 


c 3 


c 4 


c 5 


<— > 


I 


100.0 


3 4 


10 


5 






X 


X 


X 


X 


X 


X 


II 


10.0 


3 4 


10 


10 






X 


X 


X 


X 


X 


X 


III 


2.7 


5 4 


2 


5 


_ 


_ 


X 


X 


X 


X 


X 


X 




2.8 


5 4 


2 


5 




















2.9 


5 4 


2 


5 


l 


3 
















3.0 


5 4 


2 


5 


3 


3 
















3.1 


5 4 


2 


5 




















3.2 


5 4 


2 


5 




















3.3 


5 4 


2 


5 




















3.4 


5 4 


2 


5 


2 


3 














IV 


min. 3.0 


6 4 


10 


5 






X 


X 


X 


X 


X 




V 


2.7 


5 4 


2 


5 






X 


X 












2.8 


5 4 


2 


5 




















2.9 


5 4 


2 


5 


1 


3 
















3.0 


5 4 


2 


5 


3 


3 
















3.1 


5 4 


2 


5 




















3.2 


5 4 


2 


5 




















3.3 


5 4 


2 


5 




















3.4 


5 4 


2 


5 


2 


3 















Table 3.2: The 5 main steps of the parametrization procedure of ©par- The 
number of configurations for the n v random vectors and n v eigenvectors is 
given by n con f , while the number of configurations for the m v propagators is 
given by m con f. C\ - C$ indicate which of the constraints discussed in Section 
3.2.2 are imposed during the fit. Finally, the <-» symbol indicates, whether 
the parametrization procedure was iterated on the same configurations. As 
one can see, only steps III and V involve a non-linear minimization of the 
X 2 -function given in eq. (13? 
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from eq. ( 2.1 6| ) and defining 

as well as the completely symmetric combinations of the shortest paths to a 
given offset 

y(sym) = ^ + 

Vil Vm) = + y^-v + V-w + V-^-v 

^Up ^ = ^AWP + V/XjUj—p + Vp,,—V,p + VlX,—V,—p 

~\~ *—p,,h>,p V—fi,v,—p ~i~ V—(j,,—v,p ~i~ V— p,,— v,— p 

VvZiT = ^1,2,3,4 + ^1,2,3,-4 + ^1,2,-3,4 + ^2,-3,-4 

+ ^1,-2,3,4 + ^1,-2,3,-4 + ^1,-2,-3,4 + ^1,-2,-3,-4 
+ ^-1,2,3,4 + ^-1,2,3,-4 + ^-1,2,-3,4 + ^-1,2,-3,-4 
+ ^-1,-2,3,4 + ^-1,-2,3,-4 + ^-1,-2,-3,4 + ^-1,-2,-3,-4 

the overlapping block transformation in the presence of gauge fields can be 
defined as follows 

0(77) = b f Lj{U) 

with 



(sym) 
1234 

f-t p<v p<u<p 



where the coefficients q are given by 

ch = 2~ d ~\ (3.6) 

with the spacetime dimension d = 4. The parameter bf is determined by the 
engineering dimension of the field ip and is therefore fixed to bf = 2^ d ~ 1 '' 2 |94| |. 
With all these definitions the fermionic blocking kernel now reads 

Tf(U, V>, ^) = «/ Yl ( Xn B - ^n^H U )n,n B ) 

us n 

x(^-^fi(P)» B) «4). (3-7) 

n 

This leaves Kf as the only free parameter of the block transformation. It can be 
used to optimize the range of the blocked Dirac operator £> c . Using a parame- 
trization of £>fp that has couplings on the hypercube only, this optimization of 
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the interaction range becomes a very crucial point in finding a good approxima- 
tion. In order to find out for which value of kj the couplings of blocked Dirac 
operator have the fastest fall off, however, is computationally a very expensive 
thing to do in the interacting case. Therefore, we choose kj such that the 
couplings of the corresponding $Vp i n the free case have the fastest fall off, 
which is at Kf 3.45. Finally, we do not use any smearing in the definition of 
the blocking kernel even though, in principle, everything in this section is valid 
for smeared configurations. 



3.4 The Smearing 

During the parametrization of £>fp we use two different types of smearing 
schemes for the gauge fields, namely a smearing that is similar to the APE 
smearing [10S] and a smearing scheme that is closely related to the RG map- 



ping discussed in Jl09|. The APE smearing is used practically throughout the 
whole parametrization procedure, only in step V when we make the final pa- 
rametrization of ©fp on configurations with gauge couplings (3 ~ 2.7, ... ,3.4 
we use the more sophisticated RG inspired smearing. Also in the parametri- 
zation of i?Fp the RG inspired smearing is used only in the last step of the 
parametrization, whereas in the other steps the APE smearing is used. 

Modified APE Smearing 

The modified APE smearing^] is defined through the following transformation of 
the gauge links Up(n) into W^n) 

W M (n) = P SU (3) [UM + dQ M (n) + c 2 Q M (n)C/t( n + £)Q M (n)} , (3.8) 

where Psu(3) denotes the projection to SU(3) and where Q^{n) is the average 
of all the staples minus the corresponding link 



U\(n - X)U^(n - X)U x (n - A + £)] - U^{n) . (3.9) 



6 



This transformation of the original gauge links U^(n) into W^(n) has the prop- 
erty that it reproduces the trivial gauge configuration for any choice of the 
parameters c\ and C2- The resulting gauge configuration W can be used as 
starting point for further smearing steps. 

In the parametrization of £>pp we use a 2-level APE smearing. We do not 
use more smearing steps in order not to endanger the locality of the smearing. 
In steps I and II we use c\ = 0.23 and c 2 = —0.25, whereas in steps III and 



r> In the original APE smearing the third term in eq. (3^) is not present 
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IV we use c\ = —0.10 and C2 = —0.70. Using these values for the coefficients 
the fluctuations measured by the expectation value of the plaquette trace are 
reduced substantially. 



RG Smearing 

Performing a minimization and subsequently a blocking step one defines a map- 
ping of a coarse configuration V onto another coarse configuration V"[£/[V]]; a 
procedure which is called RG cycle in | 110 |. The result of the RG cycle is a 
configuration that has much smaller UV fluctuations than the original config- 
uration. It however still encodes the same long-range properties as the initial 
configuration due to the properties of the RG. The basic idea of the RG smear- 
ing is to combine the effect of a minimization and blocking step into one local 
transformation of the coarse gauge links as shown in Figure 3.5 . 

RG Smearing 




Figure 3.5: Schematic picture of the idea of the RG smearing procedure. 
The RG cycle, which consists of a minimization step from the configuration 
V to the configuration U[V] and a subsequent blocking step from the con- 
figuration U[V] to the configuration V"[Z7[V]], is approximated by the RG 
smearing, which is a local transformation of the coarse field V. 



In order to approximate such a RG cycle we use a pure gauge operator, 
which contains in addition to the simple staples also diagonal staples. These 
are built by first going in the planar or spatial diagonal directions orthogonal to 
the direction of the link /i, followed by a step in direction p, and finally returning 
along the corresponding diagonal to n + fx. To be be specific, let us first build 
the matrices 

W {0 \n,n + 0) = 1 (3.10) 
W {1) (n,n + v) = U u {n) (3.11) 

W W(n,n + v + p)= V m {S){\[U v {n)U p {n + u) + U p (n)U u (n + p)]} 

(3.12) 

W (3) (n,n + v + p + A) = P S U(3){ ^[U v (n)U p (n + v)U x (n + p + v) + perms.] 

(3.13) 
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Here u, p and A go over all directions (positive and negative) different from \i 
and from each other. In eqs. ( [3.10 )-( |3~l"3| ) the sum is taken over all shortest 
paths leading to the endpoint n' of the corresponding diagonal. Out of the 
lfW(n, n') we construct the generalized staples y( a )(ra, n + p) given by 



V (1 \n,n + fi) 
V {2 \n,n + {i) 
V {3) (n,n + fi) 



(3.14) 

n> 

(3.15) 

I W (3) (n> n')U„(n')W (3 \n' + fi,n + fi)- U^n) . 

n' 

(3.16) 

With these definitions we can now define the fuzzy link by 

3 

(n, n + fj) 

m=l 

3 

+ d ™ y{m) ( n ' n + ® U K n + ti v(m) in,n + p)}. (3.17) 



m=l 



The definition of the V^ a \n,n + /}) in eqs. (|3.14|) - (|3.16 ) ensures that for a 
trivial gauge configuration V^(n,n + jj) = and therefore W^(n) reduces to 
U^n). 

The calculation of W M (n) out of U^(n) defines one level of smearing. Using 
the resulting configuration as the starting point for another smearing step this 
procedure can be repeated as often as wanted. We denote the result of an 
fc-level smearing by W^ k \n). In our parametrization we choose the level of 
smearings to be 2, because this is clearly superior to one single smearing, but 
nearly as good as 3 levels of smearing as shown in Table 



k 


x 2 


1 


0.0133 


2 


0.0052 


3 


0.0050 



Table 3.3: The dependence of the minimal x 2 , as defined in eq. ( |3.18| ), on 
the number of smearing levels k. 



In order to parametrize the effect of a RG cycle we make a full nondinear 
optimization of the parameters c m and d m - We even allow the value of these 
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parameters to be different at each level of smearing. The x 2 -function which has 
to be minimized is given by 



X 



"conf "vol 4 



N 



(3.18) 



=1 n=l fi=l 

2 is the matrix norm and the normalization is set 



where ||M||' _„ ; 
to N = 4:N^n con {n vo \ yielding the square deviation per matrix element. The 
non-linear minimization is done with a simplex algorithm [111, 1 1 2f| . In the fit 
10 j3 = 3.0 and 5 (3 = 3.4 configurations are used. The resulting optimal 



coefficients for the 2-level RG smearing are given in Table [3.4 . 

The effect of the RG smearing on the eigenvalue spectrum of D\v is illus- 



trated in Figure 3.6. One can see that the spectrum calculated on the smeared 
configurations has clearly a smaller additive mass renormalization and that the 
spread of the low-lying eigenvalues is smaller than on the unsmeared configura- 
tion. 



k 


Cl 


C2 


C3 


di 


d 2 


d 3 


1 


0.109 


0.169 


0.098 


0.002 


0.018 


-0.014 


2 


0.309 


-0.007 


-0.068 


-0.082 


-0.121 


-0.044 



Table 3.4: The results for the coefficients c% and di in eq. (3.17) for the 2 
different levels of the RG smearing. 



3.5 The Parametrization of the Rpp Operator 

A parametrization of the operator Rpp, which shows up in the GW relation 
eq. ( |LT| ) for © = ©fp. is needed for many reasons, the simplest being just to 
check the accuracy to which £> par satisfies this very important relation. We will 
also see that a parametrization of Rpp is indispensable, because it enters the 
definition of the topological charge density eq. ( 1.29j ) and is used in the definition 



of the overlap operator with © par as we see in Chapter^; furthermore, it also 
used in the definition of the mass term of $D par used in the hadron spectroscopy 
measurements 



Using the RG equation ( p.2[ ) and the GW relation eq. (1.2) the following 
RG equation can easily be derived 

Rc(V) = — + n(U[V})R f (U[V])rt(U[V}) , (3.19) 
K f 

where the subscripts c and / stand again for coarse and fine lattice and Kf 
and Q(U) together define the blocking kernel, which is described in Section 



3|. The RG equation eq. (|l.3q) for Rfp has by far a simpler structure than 
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Figure 3.6: Eigenvalue spectrum of the Wilson Dirac operator on a FP gauge 
configuration with (top) and without RG smearing (bottom). 



the corresponding equation for ©fp. because it is a linear equation for Rfp- 
Furthermore, Rpp is an operator whose couplings vanish exactly outside the 
hypercube for the given block transformation Q(U), i.e. it is an ultra-local 
operator. Note that this is in total contrast to the RG equation for the FP 
propagator which has exactly the same form of the operator equation, 

but the FP propagator itself is an inherently non-local operator and therefore not 
suited for parametrization. That Rpp is a hypercubic operator can be seen from 



the RG equation (1.38) and the fact that Rpp is a hypercubic operator in the 
free case. A RGT involves the combination U(U)Rf(U)^V(U), which is strictly 
zero outside the hypercube, if Rf(U) is strictly zero outside the hypercube, 
i.e. such a RGT never takes R c outside the hypercube. 

3.5.1 The Details of the Fit 

The principle of the parametrization procedure for Rpp is the same as for ©fp> 



i.e. that one can simply take Figure [3J] and replace £> by R in order to have an 



idea of the basic steps of the parametrization. However, the facts that i?FP is 
a hypercubic operator that is trivial in Dirac space and that the corresponding 
RG equation ( 3.1S| ), through which it is defined, is of a much simpler structure 



than the corresponding RG equations for £>fp ( |3.1| ) and (|3^), make that the 
fitting procedure can be set up in a simpler way than the one for the Dirac 
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operator. In contrast to the Dirac operator fit we can use eq. ( |3.19 ) as it 



stands, i.e. that it is computationally feasible to perform the matrix multipli- 
cation Q(U)Rf(U)Stf(U). This allows to make use of the linear structure of 
eq. ( |3.19|) to fit the 5 different types of offsets of the hypercube (denoted by 
a) independently, i.e. that we can set up the following x 2 -functions for the fit 

4») = j}- E H a HV^) - Ri a 2(V^)\\ 2 + X X 2 C(a) , (3.20) 

( a ) i=l,n conf 

where ||M|| 2 = |Mjj| 2 is the matrix norm. The normalization is provided by 
■A/fa) = n conf n vo\ n ( a ) > where ni a \ is the number of offsets of the corresponding 
offset type. The only constraint Cr a \ on the fit is that in the formal continuum 
limit the coefficients of the free Rfp are reproduced. We choose a linear ansatz 
for the parametrization of Rfp, i.e. 

J Rpar= Yl ^00' ( 3 - 21 ) 



where the number of operators n op is 12 in the present parametrization in the 
production run code. The constraints are again linear such that the minimization 
of the x 2 -functions in eq. fl3.20 ) amounts simply to matrix inversion. Since this 



minimization can be done very fast, it is even possible to optimize the non- 



linear parameters q of the APE smearing (see Section |3.4[) . This non-linear 



optimization is performed with the simplex algorithm [ 111 , |112[| 



The 3 steps we perform in the iterative parametrization procedure, which 
are also described in Table [3~5| , are the following. In step I the free hypercubic 
^fpQ is used on the minimized = 5000 and = 100 configurations. The 
fit on the coarse lattice yields a parametrization that is used on minimized 
= 5000, = 100 and = 10 configurations in step II. This leads to the 
parametrization that is used in step III on the minimized = 3 configurations, 
where the last blocking step and fit yields the hypercubic i? par used in production 
runs. In all steps 3 4 and 6 4 configurations are used on the coarse and fine 
lattices, respectively. The constraint that fixes the coefficients in the formal 
continuum limit to the free Rpp is also kept during the whole parametrization. 
In Table 3J3 the minimal x 2 -values for the different offsets at the 3 levels of the 
parametrization procedure are given. There one can see that, like for the Dirac 
operator, the step which is by far the most difficult is the parametrization on 
the = 3 configurations. 



The free hypercubic R pal actually corresponds to the exact free Rfp, because the 
corresponding RGT equation can be solved exactly. 
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step 


P 


^conf 


X (0) 


V 2 


*(2) 


4) 


X( 4 ) 


I 


5000.0 


5 


0.000097 


0.023826 


0.003066 


0.000159 


0.000006 




100.0 


20 












II 


5000.0 


5 


0.000338 


0.091302 


0.025730 


0.000809 


0.000010 




100.0 


20 














10.0 


20 












III 


3.0 


100 


0.103887 


1.193480 


0.251903 


0.014465 


0.000786 



Table 3.5: The steps in the parametrization procedure of Rpp. The values 
X? -j given in units of 10~ 6 are the minimum of the x 2 -f unc ti° n i n e Q- ( |3.20| ) 
for the different offsets (0000), ... , (1111). In all steps 3 4 and 6 4 configura- 
tions are used on the coarse and fine lattices, respectively. The constraint 
that fixes the coefficients in the formal continuum limit to the free -Rfp is 
kept during the whole parametrization. 



Chapter 4 

Properties of the 
Parametrized FP Dirac 
Operator 



In this chapter we try to characterize our parametrizations of ©pp and -Rfp- 
From our experience in the fitting procedure, we expect the parametrization to 
perform better on gauge configurations with smaller fluctuations, i.e. at larger 
values of the gauge coupling (5. We compare properties of ® par in different 
situations. First of all we show, how our hypercubic parametrization behaves 
in the free case, where also comparisons to exact analytical results are possible. 
It is an important test case, because we do not expect a parametrization in 
the interacting case to perform better than in the free case. Furthermore, we 
discuss properties related to chiral symmetry of the parametrization used on 
minimized configurations and the production run parametrization. In addition, 
we compare the production run parametrization to the blocked Dirac operator, 
defined through the RGT in eq. (|3.l| ). Where it is possible, a comparison to 
the Wilson Dirac operator is also given. Finally, we summarize some results of 
a larger study of the hadron mass spectrum, which is presented in more detail 
in [ 55 1 . This study gives first hints for the quality of the parametrized Dirac 
operator with respect to scaling. 



4.1 Free Case 

The properties of the hypercubic approximation of £>fp in the free case are cer- 
tainly of interest, because this approximation is the starting point for the whole 
parametrization described in Chapter ^. Moreover, it shows the limitations of 
our hypercubic parametrization $D par , because it is very unlikely that a parame- 
trization will perform any better in the interacting than in the free case. The 



analytical results used in this section can be found in |94j] 
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4.1.1 Energy-Momentum Spectrum 

In the free case, one of the few physical quantities that is of interest is the 
energy-momentum spectrum. It shows how much the continuum spectrum is 
distorted by the discretization of the Dirac operator. However, the quality how 
different lattice Dirac operators approximate the continuum result varies quite 
a lot. There are simple discretizations, like the Wilson Dirac operator, which 
deviate from the continuum result already at quite small momenta, whereas ©fp 
coincides with the continuum result for all momenta. The question is, how well 



does the hypercubic parametrization perform in this respect. In Figure 4.1 the 
energy-momentum spectrum of 25 par is compared to the Wilson Dirac operator 
and to ©fp- Even though the spectrum of the truncated FP operator deviates 
from the continuum at high momenta, it is clearly closer to the continuum 
result than the spectrum of the Wilson Dirac operator over the whole range 
of momenta. Note that the high-lying branch influences the quality of the 
approximation only very little. 



3 - 



LU 




i i A » » 



- Continuum / Exact FP 

• Truncated FP 

* Wilson 



Figure 4.1: Energy-momentum spectrum of different lattice Dirac operators. 
The Wilson operator clearly deviates earlier from the continuum spectrum, 
which is also the exact FP spectrum, than the hypercubic parametrization. 
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4.2 Chiral Properties 

How much a Dirac operator deviates from exact chiral symmetry can be mea- 
sured by different quantities, such as the residual additive quark mass renor- 
malization, the occurrence of small A^A eigenvalues, the breaking of the GW 
relation on a random vector and related to this, the non-normality of the Dirac 
operator. While providing a quantitative measure for the breaking of chiral sym- 
metry, all these quantities do not show the effect of the breaking of the GW 
relation as impressively as the eigenvalue spectrum of the Dirac operator. 

We almost exclusively refer to £> par = (2i? par ) 1 / 2 ®p ar (2i? par ) 1/2 in this 
section, because it connects the two operators $D par and i? par that we have 
parametrized and it can be compared more easily to all the other (nearly) chiral 
formulations of lattice Dirac operators. We set up the following notations for 
the different types of Dirac operators used in the comparisons below/Q: 

I The hypercubic parametrization of the free FP Dirac operator together 
with the exact parametrization of the free Rfp- 

II The overlap reparametrization used on minimized f3 = 2.7 — 3.4 config- 
urations and the corresponding parametrization of Rfp- For the Dirac 
operator this corresponds to the parametrization obtained after step IV in 
the parametrization procedure and for i?pp it corresponds to the operator 
obtained after step II in the parametrization procedure (see Chapter ^). 

III The blocked Dirac operator that is obtained from the Dirac operator II 
by the RGT given in eqs. (|3.1|) and ( 3.19| ). 



IV The parametrizations 5) par and i? par used in all the production runs. 

Furthermore, we denote the Wilson Dirac operator on unsmeared and RG 

s 

w 



smeared configurations by D\y and D^, respectively 



4.2.1 Eigenvalue Spectrum 

The eigenvalues of a Dirac operator satisfying the GW relation with R = 1/2 lie 
on a circle in the complex plane with radius 1 and center at (1,0). Deviations 
from the GW relation with R = 1/2 can therefore be visualized by the eigenvalue 
spectrum: the more the eigenvalues scatter away from the circle, the more the 
GW relation is broken. Even though this is not a quantitative statement as the 
results we present later on in this chapter, it still represents the same facts. 
All the eigenvalue spectra shown in this section — apart from the free spectra, 
which can be obtained easily in Fourier space — have been calculated with the 
implicitly restarted Arnoldi method [ 113, 1 14 1 . 



In Figure we compare the parametrization in the free case with the 



production run parametrization at (3 = 3.0. The eigenvalues in the free case, 

1 The parametrizations for the ©par and _R par operators are given in Appendix ^ 



56 



Chapter 4- Properties of the Parametrized FP Dirac Operator 



which are highly degenerate, are aligned much better and lie closer to the GW 
circle than the eigenvalues of the production run parametrization, in particular 
in the region of the cut-off, where most of the eigenvalues are located. Still, 
the eigenvalues of the interacting T> vm lie fairly close to the GW circle, even 
in comparison with several other chirally improved Dirac operators [24, 32 , 97, 

In Figure \2 we show the physical branch, i.e. the small eigenvalue part, of 
eigenvalue spectrum of our parametrization II for a minimized j3 = 3.0 configu- 
ration on a 10 4 lattice. The eigenvalues lie very close to the circle. Performing 
a RGT, the eigenvalue spectrum of the resulting Dirac operator III on the cor- 
responding 5 4 lattice is roughly of the same quality, as the spectrum on the 
minimized configuration. This shows that the Dirac operator, which we approx- 
imate in the last step of the parametrization procedure, has very good chiral 
properties. In comparison we show the eigenvalue spectrum of the final para- 
metrization IV on the same 5 4 configuration. In this spectrum the alignment of 
eigenvalues to the circle is clearly inferior to the one of the blocked operator, 
showing the difficulty of parametrizing ©pp on rough QCD configurations. But, 
one has to keep in mind that the final parametrization only has a very small ad- 
ditive mass renormalization and compared to the Wilson Dirac operator, which 
is shown in the lower part of this figure, the fluctuations of the small (real) 
modes are reduced considerably. 

The eigenvalue spectrum of £> par does not lie on a circle (see Section 1.2.2 ) 
and therefore can not as easily be used to get an impression of the amount of the 
breaking of chiral symmetry as the eigenvalue spectrum of D par . But, because 
© par is the fundamental operator in the parametrization procedure, we show its 
eigenvalue spectrum on 8 4 lattice in the free case, where the eigenvalues are 
highly degenerate, and on a 4 4 lattice in the interacting case at (3 = 3.0 in 
Figure 4.4. 



4.2.2 Breaking of the Ginsparg- Wilson Relation 

The FP Dirac operator has an exact chiral symmetry on the lattice. The para- 
metrization 35p ar , however, satisfies the GW relation no longer exactly. In order 
to quantify the degree to which it breaks the GW relation, we use the following 
measure, which is also used in the literature [119, 120| ]p1 



A GW = (HCPpar +^ Par " Pp ar PparH|} , (4.1) 

where v is a normalized random vector. The brackets denote the expectation 



value and | \v\ | is the usual vector norm. Table [4JJ compares Aqw °f t ne Dirac 
operators l-IV. Let us give a few comments to this table: Aqw is smaller by 
one order of magnitude in the free case (operator I) than for the production run 



2 Note that in some publications also the square norm 
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Figure 4.2: The eigenvalue spectrum of the truncated £>fp in the free case 
on a 8 4 lattice (top) and the eigenvalue spectrum of V par at (3 = 3.0 on a 4 4 
lattice (bottom). The eigenvalues in the free case are highly degenerate and 
therefore it seems that there are less eigenvalues on the larger lattice than 
on the small lattice, which is however not true. The figures show clearly, 
that the GW relation is broken to a lesser extent in the free than in the 
interacting case; although the eigenvalues in the interacting case are still 
close to the circle. 
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0.5 



Figure 4.3: Top row: In these figures the physical branch of the eigenvalue 
spectra on a minimized 10 4 configuration and the corresponding coarse 5 
configuration at (3 = 3.0 are shown. The eigenvalues on the minimized 
configuration (left), which are calculated with operator II, are clearly aligned 
to the GW circle, which is also the case for the blocked operator III (middle). 
The alignment of the eigenvalues of the final parametrization IV (right) is 
clearly worse than the one of the blocked operator III. However, essential 
features like the small additive mass renormalization are conserved. Bottom 
row: The eigenvalue spectrum of the Wilson Dirac operator on the same 
5 configuration at (3 = 3.0 is shown without (left) and with RG smearing 
(right). The additive mass renormalization, which is roughly 1.2 and 0.6, 
respectively, is removed for this comparison. The fluctuations of the small 
(real) modes are clearly larger than for T> vax ] however, one has to be aware 
that it is only meaningful to compare the fluctuations close to the origin, 
because the spectrum of a Dirac operator satisfying the GW relation with 



a general R does not necessarily lie on a circle (see Section 1.2.2 ). 
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Figure 4.4: Eigenvalue spectrum of the truncated $Dfp i n the free case on 
a 8 4 lattice (top) and the eigenvalue spectrum the production run © par at 
(3 = 3.0 on a 4 4 lattice (bottom). This figure shows clearly that an operator 
which satisfies (approximately) the GW relation with R not proportional to 
1 does not lie on a circle (see Section 1.2.2| ). Notice that the eigenvalues in 
the free case are highly degenerate. 
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Type 





5 4 


8 4 


10 4 


I 


free 


0.0293(1) 


0.0294(1) 


0.0294(1) 


II 


min. 3.4 


_ 


_ 


0.0525(1) 




min. 3.2 


_ 


_ 


0.0530(1) 




min. 3.0 






0.0553(1) 




min. 2.7 






0.0604(1) 


III 


3.4 


0.0298(1) 








3.0 


0.0308(1) 








2.7 


0.0339(3) 






IV 


3.4 


0.1782(2) 


0.1779(2) 


0.1783(1) 




3.2 


0.1885(2) 


0.1938(3) 


0.1947(1) 




3.0 


0.2017(2) 


0.2176(4) 


0.2172(1) 




2.7 


0.2819(2) 


0.2816(3) 


0.2811(1) 



Table 4.1: The breaking of the GW relation Agw 5 as defined in eq. (|4.1| ), is 
shown for the different types of Dirac operators defined at the beginning of 
this section. The statistical error is given in the brackets. 

parametrization (IV), whereas the parametrization on the minimized configura- 
tions (II) performs only slightly worse than the free parametrization on trivial 
configurations (I). The blocked operator (III) is better than the parametrization 
on the minimized configurations (II) and almost reaches the same quality as the 
free parametrization (I). The increase in the breaking of the GW relation with 
decreasing (3 (i.e. increasing lattice spacing) is as expected. The fluctuations of 
Aqw f° r different vectors and configurations are very small and finally, Aqw 
does not show any dependence on the lattice size, as expected. 

4.2.3 Breaking of Normality 

A Dirac operator satisfying the GW relation is normal with respect to R, which 
means 

© t i?S = ®ii® t and VV^ = V^V, (4.2) 

respectively. Hence, V and are commuting, which implies that the eigenvec- 
tors of V are also eigenvectors of . A further consequence is that [2?, 75] = 
in the subspace of the real eigenmodes and therefore the real eigenmodes of a 
75-hermitian V have definite chirality. Thus, the magnitude of the breaking of 
normality indicates to what extent e.g. the zero modes are chiral. We measure 
the non-normality of a Dirac operator D by 

A n = <||(D£)+-DtD) v ||) ) (4.3) 

where v is a normalized random vector. Table shows A n for the Dirac 
operators ll-IV. The free case is not indicated, since free Dirac operators are 
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Type 


P 


5 4 


8 4 


10 4 


II 


min. 3.4 


_ 


_ 


0.0196(1) 




min. 3.2 


_ 


_ 


0.0199(1) 




min. 3.0 






0.0210(1) 




min. 2.7 






0.0235(1) 


III 


3.4 


0.0223(2) 








3.0 


0.0251(2) 








2.7 


0.0312(5) 






IV 


3.4 


0.0946(2) 


0.0950(1) 


0.0952(1) 




3.2 


0.1027(4) 


0.1043(2) 


0.1050(1) 




3.0 


0.1163(6) 


0.1197(2) 


0.1196(1) 




2.7 


0.1663(10) 


0.1662(3) 


0.1680(2) 



Table 4.2: The breaking of the normality A n , as defined in eq. (|4.3|), is 
shown for the different types of Dirac operators defined at the beginning of 
this section. 



normal. The blocked operator (III) performs slightly worse than the parametri- 
zation on the minimized configurations (II), the difference is however very small. 
The fluctuations of A n for different vectors and configurations are again very 
small. Finally, A n does not show any dependence on the lattice size, while the 
dependence on the lattice spacing is again as expected, i.e. the larger the lattice 
spacing, the larger the breaking of normality. 

4.2.4 Eigenvalues of A^A 

The GW relation eq. ( |1.20|) is equivalent to 

A* A = 1, A = l — T>, (4.4) 



i.e. that A is a unitary operator. Hence, the breaking of the GW relation is 
signaled by the occurrence of A^A eigenvalues deviating from 1. In Figure 4.E 
we show the density of the A^A eigenvalues for £> par at = 3.0 on a 4 4 lattice. 
The eigenvalues are distributed between 0.5 and 1.5 with a well-defined peak 
at 1. For larger physical volumes the distribution typically will have longer tails, 
which eventually leads to the occurrence of very small A* A eigenvalues. 

The ratio of the smallest A m i n to the largest eigenvalue A max of A' A is 
related to the asymptotic convergence of approximations to the overlap con- 
struction [ 33 1 . Hence, this ratio can be used to indicate how expensive the 
overlap construction is approximately for a certain Dirac operator and gauge 
action. In Figure 4J5, we plot the 50 smallest eigenvalues of A' A divided by the 
largest eigenvalue of A* A for 5 different gauge configurations at gauge coupling 
(3 = 3.0 and /3 = 3.5, respectively. We see that the A^A eigenvalues are much 



62 



Chapter 4- Properties of the Parametrized FP Dirac Operator 





3.0 


3.5 


D w 
u w 


16130(212) 
6176(212) 
682(18) 


17174(154) 
8726(122) 
944(50) 



Table 4.3: The average number of Dirac operator times vector multiplica- 
tions to find the 50 lowest A^A eigenvalues for D par and the Wilson Dirac 
operator on smeared (Dyy) and unsmeared (Dw) gauge configurations. This 
table shows that calculations using D par do not necessarily have to be as 
expensive as one would naively think, because the large difference in the 
number of matrix vector multiplications between P par and Dw almost com- 
pensates for the overhead of T> paT in one matrix vector multiplication. 



closer to 1 for 2? par than for Dw at both gauge couplings. At this point it 
should be mentioned that relatively to Dw the calculation of the A* A eigen- 
values, which is an important ingredient in the numerical implementation of the 
overlap Dirac operator (see Chapter ||), is by far less expensive for D par than 
one would naively think. Because P par is substantially closer to be a normal op- 
erator and because the A^A eigenvalue spectrum is less dense at the lower edge 
than for the Wilson operator, the Arnoldi and the Ritz functional method find 
the eigenvalues with considerably less matrix vector multiplications. The differ- 
ence is particularly large between Dw on unsmeared configurations and P par , 
as shown in Table A3, where the number of operator times vector products is 
given for the Arnoldi method^. At the upper edge of the spectrum the difference 
is much smaller and amounts to roughly 200 matrix vector multiplications for 
P par and 800 for Dw, independent of smearing. 

Let us finally remark that one might try to use the Chebyshev acceleration 
technique described in [ 12 1 1 to map the low-lying part of the A* A eigenvalue 
spectrum onto a larger interval, which in certain cases speeds up the calculation 
of the eigenvalues and eigenvectors |122|. It would be very interesting to know 
howD par and Dw, respectively, can profit from this technique for the calculation 
of the low-lying A* A eigenvalues. 



4.2.5 Index Theorem 
For GW fermions the Atiyah-Singer index theorem from eq. (1.27) holds. Re- 



calling the previous discussion about the distribution of the iU eigenvalues we 
can not expect the index theorem to hold for 2? par because it can be rewritten 



3 A comparison of T> paT and Dw with the same conclusion for the Ritz functional is 
given in 
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Figure 4.5: Density p{\) of the A^A eigenvalues of V pai on a 4 4 configuration 
at (3 = 3.0. On larger volumes the probability of small A^A eigenvalues 
increases, which makes that a tail on the left hand side of the peak develops. 
At the upper edge of the spectrum, however, the structure of the distribution 
does not change significantly at larger volumes. 
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Figure 4.6: The ratio of the 50 smallest A^A eigenvalues to the largest A^A 
eigenvalue on 5 different 12 4 configurations at (3 = 3.0 (top) and (3 = 3.5 
(bottom) is shown for T> par , D w and D\y. As expected this ratio is clearly 
closer to 1 for 2? par than for Dw- Notice that the legend is valid for both 
figures. 
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as follows 

index(P) = ^Tr(7 5 P) = ~Tr(^ 5 (l-V)) = -lTr( l5 A), (4.5) 

where we use TV 75 = 0. Hence, small eigenvalues of 75A and, in particular, 
an unsymmetric distribution of 75A eigenvalues^ tend to destroy the correspon- 
dence given by the index theorem. We will discuss our definition of the index of 
D par , which does not satisfy the GW relation exactly, and therefore has no exact 
zero modes with definite chirality in Section |6.1.1| . In Figure 4J we see that the 



relation is indeed distorted heavily for the production run parametrization. The 
definition of the index using the trace leads to a result which on average is too 



arge by a factor of 3.16. In Figure 4.S we compare the index of £> par with the 



index of the corresponding overlap operator constructed from T> paT (see Chap- 
ter |5|), which has a very precise chiral symmetry. We see that the index of the 
Dirac operator defined through the low-lying real eigenmodes and their chirality 
is a much more stable quantity than the index defined through 1/2 Tr(75 2? par ), 
because it coincides with the index of the overlap Dirac operator in most of 
the cases. We will use this observation in Chapter ||, where we calculate the 
topological susceptibility. 

In contrast to the production run parametrization the trace definition of 
the index works well on the minimized configurations when the corresponding 
Dirac operator is used (definition II), despite the fact that the wrong sign in 
the normalization condition for the topological charge density was used in the 



parametrization (see Section [3.2.21) . This means that the linear polynomials in 
the couplings really could straighten out the error in the normalization condition 
and therefore we can also conclude that the problems with the trace definition 
of the index, we see in the production run parametrization, are due to other 
problems. Maybe an even richer parametrization for the pseudoscalar part is 
needed to give an accurate description of 1/2 Tr(7s P par ) on configurations 
used in simulations. 



4.2.6 Residual Additive Mass Renormalization and Fluctu- 
ations of the Low-Lying Modes 

Another very important property of a Dirac operator satisfying the GW relation 
is the absence of an additive mass renormalization and related to this, the 
protection against exceptional configurations at non-zero quark mass. This 
means that the pions get massless if the bare quark mass is zero, making that 
no fine tuning of parameters is needed to reach the chiral limit [39|. Even more 
important is the fact that there are no fluctuations in the low-lying modes of 
the Dirac operator that make the simulation of (quenched) QCD in (or at least 



4 The distribution of the 75^ eigenvalues is symmetric for a GW Dirac operator, apart 
from the eigenmodes of T> at and 2, which lead to an asymmetry in the distribution that 
is needed for the index theorem as one can easily show. 
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Figure 4.7: The correlation between the index of 2? par and 1/2 Tr(7s D par ) 
is shown for 50 8 4 configurations at (3 = 3.0. One can see that a certain cor- 
relation between the two quantities remains, but the index theorem clearly 
does not hold for the parametrized operator. The definition with the trace 
gives a result which on average is too large by a factor of 3.16 as indicated 
by the dashed line. The solid line shows, where the points actually should 
lie if the index theorem would be satisfied. 
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near) the region of the physical up and down quark mass impossible. For ® par 
there is a residual additive mass renormalization and there are fluctuations in 
the small eigenmodes that do not allow to simulate at an arbitrarily small quark 
mass, but the situation is definitely improved compared to the Wilson Dirac 
operator, where the fluctuations in the low-lying modes make the simulations 
difficult at a n/p ratio of 0.5 to 0.4, depending on the lattice spacing and 
the volume [ |123[ |124-| ] . The clover improved Wilson Dirac operator is even 
worse in this respect, i.e. that the low-lying eigenvalues fluctuate even more. 
Exceptional configurations were already found at a ir/p ratio of 0.54 at a lattice 
spacing of 0.1 fm | |125| |. Finally, the situation does not seem to improve as 
much as expected in the case of full QCD simulations |126|. In our hadron 
spectroscopy study, which is described in detail in p5| , we were able to simulate 
down to a n/p ratio of 0.28, which corresponds to a quark mass of roughly 10 
MeV in the MS scheme at p = 2GeV; notably, without facing the problem of 
exceptional configurations in an ensemble of 200 16 3 x 32 configurations at a 
lattice spacing a = 0.155 fm (f3 = 3.0). From this study we can extract the 
residual additive quark mass renormalization aA m and an upper limit for the 
minimal bare quark mass am m j n that is needed to prevent from the appearance 



of exceptional configurations. The results are given in Table 4.4 and show that 
the additive mass renormalization and the fluctuations of the low-lying modes 
are small. The result on the largest physical volume at f3 = 3.0 is the most 
impressive as the bare mass can be chosen to be very small. 





aA m 




3.0 


-0.0006(4) 


0.013 


3.4 


-0.0180(4) 


0.029 


3.7 


-0.0194(4) 


0.0235 



Table 4.4: The additive mass renormalization and the smallest bare mass 
used in the hadron spectroscopy simulations in J6q ]. The values at f3 = 3.0 
are from a larger physical volume and therefore it is possible to go to even 
smaller quark masses. 



4.3 Scaling Properties of the Parametrized FP Dirac 
Operator 

The FP Dirac operator is expected to have small scaling violations. This is 
a crucial property that allows one to simulate at larger lattice spacings with 
only very small distortions of the continuum physics. We do not cover the 
scaling properties of £> par as detailed as the chiral properties; the reason for 
this is not the significance of scaling properties, but simply the fact that it 
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is computationally much more expensive to get information on scaling. This 
explains why there is not much data on scaling properties, yet. 

Presently, our only source of data on the scaling behaviour of © par is the 
hadron mass spectroscopy measurement described in [55|. The simulations are 
performed at 3 different lattice spacings. The extent of the lattice is chosen such 
that all 3 measurements are done in the same physical volume with a spatial 
extent of L « 1.22 fm, in order to allow for a continuum extrapolation. This 
volume is rather small, which makes that the masses of the particles (especially 
of the baryons) are distorted by finite volume effects [127|. But because the 
finite volume distortions are physical effects, the data is nevertheless suited for 
a scaling analysis, even though uncertainties in the scale determination lead to 
larger uncertainties in the scaling analysis than for larger volumes, where the 
dependence of the particle masses on the volume is negligible. The data at j3 = 
3.0,3.4 and 3.7 for the vector meson, the decuplet and octet baryons is shown 
in Figure |4~9|. It indicates a remaining 0{a 2 ) scaling violation. However, a more 
careful analysis is needed, because there are several sources of systematic errors. 
Most significantly is probably the uncertainty in the scale, as the interpolating 
formula the Sommer scale ro in [^5, 96] covers the range = 2.361, . . . , 3.4 and 
therefore an extrapolation for the value at (3 = 3.7 has been used. Furthermore, 
this scale determination has a systematic error which is not accounted for in 
Figure 4J. The data is not accurate enough to completely exclude 0(a) scaling 
violations - which are definitely not present in the hadron spectrum for a GW 
Dirac operator [52] - but, if there are remaining O(a) scaling violations, then 
they are very small. This is actually supported by the fact that also the energy- 
momentum dispersion relation for the pseudoscalar and the vector meson do 
not show any visible cut-off effect, as shown in Figure 4.1C. 



4.4 Discussion 

The parametrized FP Dirac operator, which we obtained from the fitting pro- 
cedure in Chapter [3], is clearly improved in many respects in comparison to the 
Wilson Dirac operator. It has a relatively small breaking of chiral symmetry 
that allows to calculate the hadron mass spectrum down to quark masses of 
roughly 10 MeV in the MS scheme at /x = 2GeV. Even though the physical 
point can not yet be reached with the present parametrization, it is a remark- 
able progress in comparison to Wilson fermions. Furthermore, it is possible to 
go to such small quark masses that there is a clear signal for quenched chiral 
logarithms, even in the data of the pseudoscalar spectroscopy with degenerate 



quark masses [55|. The comparison of the parametrization with the blocked 
operator obtained from the RGT shows the limitations of the present parame- 
trization, because the blocked operator has smaller breaking of chiral symmetry. 
Furthermore, the present parametrization does not satisfy the index theorem. 
But one should keep in mind that the topological charge density is a quan- 
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Figure 4.8: The correlation between the index of 2? par and the corresponding 
overlap construction with 2? par is shown for 50 8 4 configurations at (5 = 3.0. 
One can see that the value of the index defined by 2? par coincides on most 
of the configurations with the overlap definition. 
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Figure 4.9: Scaling of the hadron masses (vector meson, octet and decuplet 
baryon) calculated with $D par at 3 different lattice spacings. The data sug- 
gests that the remaining scaling violations are mostly 0(a 2 ) effects and are 
rather small. 
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Figure 4.10: Energy-momentum dispersion relation E(p) for pseudoscalar 
(P) and vector mesons (V) on the 16 3 x 32 lattice at = 3.0. The solid lines 
show the continuum value with E(0) given by the measured mass. 



tity which is of small order and in the present parametrization no particular 
effort was undertaken to make the fitting procedure more sensitive to the pseu- 
doscalar contribution and therefore some improvement might be achieved in 
future parametrizations. The data for the scaling of the hadron masses sug- 
gests that the present parametrization leads to remaining cut-off effects which 
are mainly 0(a 2 ) and as shown in [65] the cut-off effects are on the order of 
various other 0(a) improved Dirac operators. As discussed in Section A3, the 
sources of systematic errors are too large to make a clear statement on the size 
of the 0(a 2 ) scaling violations, apart from the fact that they are rather small. 
The data, however, indicates that order 0(a) cut-off effects are very small, as 
one can see no sign of their presence. These statements about the small scaling 
violations are supported by the fact that the energy-momentum dispersion rela- 
tion for the pseudoscalar and the vector meson do not show any visible cut-off 
effect. 



Chapter 5 

Overlap Construction with 
the Parametrized FP Dirac 
Operator 



One of the very important properties of ©fp is its exact chiral symmetry on 
the lattice [128]. Unfortunately, this exact symmetry gets lost through the 
parametrization. As expected, the violation of chiral symmetry is much smaller 
for £> par than for Wilson fermions. In addition, Neuberger's overlap construction 
[40,41] gives a recipe, how a Dirac operator that respects chiral symmetry on the 
lattice can be constructed from any Dirac operator that has no doublers in the 
continuum limit and is local [129]. In numerical simulations approximations to 
Neuberger's construction can be used to obtain Dirac operators whose violations 
of chiral symmetry are, in principle, below an arbitrarily chosen limit. Very 
precise chiral symmetry is required for some measurements, such as the chiral 
condensate 1 130 1 . However, the decision to what extent chiral symmetry is 
needed for a particular calculation has to be taken from case to case, since 
(approximately) chiral actions are computationally very expensive [98, 120 1 . 

In this chapter we first discuss Neuberger's overlap formula and then, in 
particular, the properties of the overlap construction with £> par and i? par . The 
actual construction with R ^ 1/2 is numerically somewhat more involved than 
the standard case R = 1/2 and we give the technical details in Appendix]^. Be- 
low, we explain the details of the Legendre expansion of the operator (A^A)^ 1 / 2 
in the overlap formula, because in our simulations with the overlap we exclu- 
sively use this type of approximation. We present results on the properties of 
the overlap construction with ® par and -R par and make some comparisons to 
the standard overlap construction with the Wilson Dirac operator. Finally, we 
conclude with a remark on scaling properties of the overlap construction with 

^par ■ 
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5.1 General Overlap Construction 

The overlap formula in terms of a Dirac operator D and the local gauge operator 
R from the GW relation eq. ( |1.1[ ) in the massless case is given by 

V = i _ A(AtA)-V2 = 1 - 75 e(75A) (5.1) 

with 

A = l + s— (2R) 1 / 2 D(2R) 1 / 2 , (5.2) 

where e{^A) is the matrix sign function of the operator 75AQ. The real param- 
eter s is important for the overlap construction with the Wilson Dirac operator, 
where it can be used to optimize e.g. the localization range or the convergence 
rate of an approximation [131]. In connection with the overlap with V pav the 



parameter s does not play a role, because there is no need to shift the A^A 



eigenvalues, since they are already centered around 1 as shown in Figure 4.5. 
It is easy to see that V satisfies the GW relation. More generally, it possible to 
define 

© = (2i?)- 1 / 2 (l - A(A^A)- 1 / 2 )(2Ry 1 / 2 = (2i?)- 1/2 (l " 7 5 e(75^))(2i?)" 1/2 , 

(5.3) 



where the operator A is again as defined in eq. (|5.2] ), such that satisfies 
the GW relation with a general R. The inverse square root of the hermitian 
operator A' A or the matrix sign function of 75^ are the objects that cause all 
the trouble in the overlap formula. They have to be approximated in a certain 
way, because a full decomposition of these operators into their eigenmodes is 
computationally not feasible, even for modest lattice sizes. In the last 5 years 
many different approximation schemes have been worked out and applied to 
this problem. One of the main directions of approximation schemes deals with 



polynomial expansions, like Legendre [131[, Gegenbauer [ 132 1 or Chebyshev 



polynomials [ |12C| ], while the other main direction is concerned with rational 
approximations, like the Remes method [ 133 1 , the polar decomposition [|134fl or 



the Zolotarev polynomials |98, 135 1 . Furthermore, there are also direct Krylow 



subspace methods used to approximate the matrix sign function [136, 137 1 and 
some attempts to formulate algorithms in a 5-dimensional space in order to 
avoid nested Krylow methods, which are used in the rational approximations. 
Such 5-dimensional approximations have a close relation to the domain-wall 



fermion approach [|138| , |139|1 



It is useful to introduce the mass term of a Dirac operator satisfying the 
GW relation in a particular way, namely 

V{m) = (l - y)^(0) + m (5.4) 



lr The matrix sign function of 75^4 and the square root of A' A are well defined, because 
75^4 and A^A are hermitian operators and therefore the matrix functions can be defined 
in terms of the eigenvalues and eigenfunctions of the corresponding hermitian operator. 
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and 

8>( m ) = (l - ^ ©(0) + m(2 J R)- 1 , (5.5) 
respectively. With this definition the scalar density S 1 is defined as follows 

S = ^(l-^Z>(0))V (5.6) 

and 

5 = Vi((2ii)" 1 -i©(0))^ (5.7) 

for the general case. Defined this way the scalar density transforms under chiral 
transformations like the corresponding continuum density ||140(| . This implies 
that this density, which is a dim=3 operator, is 0(a) improved, since it does 
not mix with any dim=4 operator. A detailed discussion of these densities and 
the construction of chiral covariant currents, which are also 0(a) improved, is 
given in Hl40f| . 

Our choice of the Legendre expansion might surprise a little bit when con- 
sulting the literature about the different approximations, because the Legendre 
expansion does not figure among the most efficient schemes [pq] . However, 
the discussion of these approximations is almost entirely concentrated on the 
overlap construction with Wilson fermions, which is completely different due to 
the fact that the operator A^A deviates much more from 1 than in the case of 
"Dpar (see Figure and therefore high order polynomials have to be used in 
the approximations. On the contrary, P par is already quite close to satisfy the 
GW relation and therefore low order polynomial approximations can be used to 
approximate e.g. the operator (A^ A)~ 1 / 2 . For low order approximations, how- 
ever, the difference between the different polynomial approximations is much 
less pronounced than at high orders. In fact, our experiments with the Cheby- 
shev expansion didn't show a visible advantage for the Chebyshev expansion at 
polynomial orders of < 4, but clearly more precise studies are needed to clarify 
this issue. 

5.1.1 Legendre Expansion 

Following the original work [ |131| ], we give a short description of Legendre ex- 
pansion of the operator (A^ A)^ 1 / 2 . 

In order to ensure convergence of the Legendre expansion we assume that 
the bounds 

u < (if>\A^ A\ifi) < v (5.8) 
hold for some strictly positive constants u < v and arbitrary normalized states 
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The Legendre polynomials Pk(z) may be defined through the generating 
functional 

oo 

(l-2tz + t 2 )- 1 / 2 = Y,t k Pk(z). (5.9) 

k=0 

Usually z is taken to be a number, but ( |5.9| ) remains meaningful if we substitute 

v + u - 2A^A 



v — u 



(5.10) 



The expansion ( |5.9| ) is convergent for all t satisfying \t\ < 1. 
Introducing the parameter 9 through 

v ~\~ u 

cosh0 = , 0>O, (5.11) 

v — u 

and t = e~ 6 (^) takes the form 



{rfAy l l 2 = KY^t k P k (z), K = J-^-, (5.12) 



k=0 



SI 



nee 1 — 2tz + t 2 is proportional to A^A for this choice of t. In practical 



applications the sum in (5.12) will be truncated after a finite number of terms 
and this then defines the Legendre expansion of the overlap operator to a certain 
order N, which we denote by P' w ' or JD^, respectively. 

5.1.2 Exact Projection of A^A Eigenvalues 

For large volumes and small values of j3 the eigenvalues of A* A tend to scat- 
ter further away from 1, even for P par . This eventually leads to big problems 
in the numerical expansions used for the inverse square root or the sign func- 
tion, because the parameter c governing the asymptotic convergence of these 



expansions is |33|] 



\\V^ -V\\ = e~ cN , c=J2m, (5.13) 

V ^max 

where A m i n is the smallest eigenvalue of A' A and A max the largest. Because the 
distribution of the eigenmodes of A^A is typically not very dense at the lower 
edge of the spectrum, one can use exact projection methods to decrease the 
condition number c substantially. The same is not true for the upper edge of 
the spectrum of A' A, because the modes lie much denser there. Let us specify 
what we mean by the exact projection of the n eigenmodes of with 
the eigenvalues \ of smallest modulus: 

n 

(aU)-v 2 = ^2j^mm+viA VP [(A^Ar 1/2 ]vi, (5.14) 
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75 



for the square root and 



e( l5 A) = J>(A<)I^>M +n App[e(75A)]^ 



(5.15) 



i=l 



for the matrix sign function, where the projection to the subspace orthogonal 
to the lowest n eigen modes is defined by 



(5.16) 



i=l 



and where App stands for one of the possible expansions which may be used 
in practice. The calculation of the A' A eigenvalues and eigenmodes can be 
performed e.g. with the Ritz functional [104, 105] or the implicitly restarted 



Arnoldi method [113, 114]. In our overlap studies we mainly use the Arnoldi 
method to calculate the eigenmodes, because it showed to be efficient in the 
search of rather large numbers (< 200) of small A^A eigenvectors. 



5.2 Properties of the Overlap with FP Kernel 

In the following discussion of the properties of the overlap construction with 
X? par as kernel we focus on two topics: the locality of the resulting operator and 
the dependence of the remaining breaking of chiral symmetry on the precision 
of the approximation. 

5.2.1 Locality 

A lattice action must be local for the continuum results to be universal, i.e. in- 
dependent of the details of the lattice action. Locality means that fields at 
large separations r = \y — x\ 3> 1 have an exponentially small coupling p(r), 
i.e. p(r) ~ exp(— vr) with v = 0(1). Optimizing the locality is essential so 
that e.g. the exponential fall-off of correlation functions can be separated from 
direct couplings in the lattice action even on coarse lattices. The locality of a 
lattice Dirac operator can be measured by 

f p (r) = m&x{\\Dv\\, \\y - x\\ p = r}, (5.17) 

where v is a vector with point source at x and \\v\\ p = (%2i=i n ' s tne 

vector norm. We use the usual vector norm (p=2) and the taxi driver norm 
(p=l) to discuss the data. These two norms are also the most common choices 
in the literature [116, 1 3 1 1 . In Figure 5.1, we present the locality measured 



by the expectation value of the normalized function /2(^)//2(0) of the overlap 
construction with 2? par and the Wilson Dirac operator on smeared {Dyy) and 
unsmeared configurations (Dw) at (3 = 3.0 and (3 = 3.5 on 12 4 lattices. We 
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3.5 


3.0 


p 


1 


2 


1 


2 


D w 


0.48(1) 


1.14(5) 


0.42(1) 


0.94(3) 


u w 


0.58(1) 


1.32(5) 


0.55(1) 


1.22(4) 




0.81(2) 


1.60(2) 


0.82(1) 


1.60(2) 



Table 5.1: Value of the exponent v as defined in eq. ( 5.18| ) for the overlap 
construction from eq. (5.1) with different Dirac operators and for the taxi 
driver norm (p=l) as well as the usual vector norm (p=2), which gives the 
locality in terms of physical distances. 



use a Legendre expansion of order 250 for the Wilson overlap and order 25 for 
the overlap with D par and project out the 50 lowest A* A eigenmodes to make 
sure that the deviations from the exact overlap Dirac operator are smaller than 
the couplings at r/a = 12. The locality for the Wilson overlap is optimized by 
adjusting the parameter s in eq. (|Q|). The overlap construction with £> par is 
clearly more local than the Wilson overlap, independent of smearing. While the 
exponential fall-off of the T> pax overlap stays the same on both gauge couplings, 
this is not the case with the Wilson overlap, whose locality gets slightly worse 
with stronger coupling. In Table [Tl] we list the exponents v obtained from an 
exponential fit to the expectation value of f p (r) 

(f p (r)) oc e — ' a , (5.18) 



where r/a > 13 for p = 1 and r/a > 7 for p = 2. The results in Table 5T show 
that with the overlap with £> par one can investigate masses up to ma « 1.5 
in spectroscopy calculations. Finally, the overlap construction from eq. (|5.3|) 



with ® par and i? par has the same exponential fall-off as the construction from 



eq. (5T) 



5.2.2 Chiral Properties 

In this section we consider the breaking of chiral symmetry defined through 
various quantities. We use again the breaking of the GW relation Aqw(N) and 
the breaking of normality A n (N) on normalized random vectors, as defined in 



eqs. (4.1) and (|4.3j) ; however, with the additional dependence on the order N 
of the Legendre expansion. But first we try to characterize the deviations of the 
eigenvalues from the GW circle. 

The eigenvalue spectrum is suited to show deviations from the GW circle. 
However, when the GW breaking Aqw drops below ~ 10~ 2 , it gets difficult to 
judge the quality of the approximation by eye, which makes it necessary to use 
a different measure to visualize the deviations from the circle. Hence, we use 
the fact that if the eigenvalues A of V lie exactly on the GW circle, then the 
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Figure 5.1: Locality of the overlap Dirac operator using Dyy, Dyy and 2? pa r 
on 12 4 configurations at ft = 3.5 (top) and f3 = 3.0 (bottom). The over- 
lap construction with 2? par is clearly more local than the Wilson overlap, 
independent of smearing. While the exponential fall-off of the 2? par overlap 
stays the same at both gauge couplings, this is not the case for the Wilson 
overlap, whose locality gets worse with stronger coupling. 
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stereographic projection 



A 



1 - * 

1 2 



(5.19) 



is purely imaginary. In Figure 5.2, we see that Re(A) falls off exponentially as 
we increase the polynomial order, showing that the eigenvalues A get closer and 
closer to the GW circle. The most rapid decrease is at the finest lattice spacing 
a as 0.10 fm {(3 = 3.4), which is what one expects. 

In Figure [T3] the breaking of the GW relation and the breaking of the 
normality is shown on 10 4 lattices at f3 = 3.0, (3 = 3.2 and (3 = 3.4. First of 
all we see that the breaking Agw(-W) and A n (N) fall off exponentially as we 
increase the order N of the Legendre polynomial approximation of (A^ A)^ 1 / 2 . 
Treating the the same number of A* A eigenvalues exactly, the fall-off is clearly 
the steepest for the smallest lattice spacing ((3 = 3.4), which is what one 
expects. When the number of projected A* A eigenvalues, however, is decreased 
such that the ratio A m i n /A max is approximately the same on all lattice spacings 
(this actually corresponds to projecting only the 20 lowest eigenvalues of A' A at 
(3 = 3.2 and the 5 lowest at (3 = 3.4) the approximation is of the same quality to 
a very high degree. In Table |5]2| we collect the exponents fi describing the fall- 
off of Aqw(A0 with increasing order of the Legendre expansion and compare it 



to the predicted asymptotic value c from eq. ( 5.13 ). Notice that the exponent fi 
is clearly larger than the asymptotic value in all cases, which indicates that the 
regime of asymptotic convergence is by far not reached, making the convergence 
in the region that is interesting for simulations much faster than in the worst 
case. But, all the same the value c seems to be in clear correlation with the 
measured fall-off \i. The results are based on rather low statistics and in order 
to make the statement more quantitative clearly more statistics is needed. 



p 




c 


3.4 


1.66 


0.67 


3.2 


1.51 


0.63 


3.0 


1.26 


0.55 



Table 5.2: Value of the exponent /i giving the exponential fall-off of the 
breaking of the GW relation in terms of Aqw(-W) is compared to the asymp- 
totic fall-off given by the condition number c, defined in eq. (5.13). The ex- 
ponential fall-off of Agw(A0 is clearly steeper than the asymptotic fall-off, 
but c is obviously correlated to fi. The data is obtained from 10 4 lattices. 
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Figure 5.2: Deviation of the P^' eigenvalues from the GW circle defined 
by the stereographic projection in eq. ( |5.19 ). 



5.3 A Final Remark 

The overlap construction with P par can be looked upon as a tool to cure the 
breaking of chiral symmetry that remains in the parametrized FP Dirac op- 
erator. Clearly, one hopes that the good scaling properties of the FP Dirac 
operator, which to a certain extent are still present after the parametrization 
process, are not distorted too much by the overlap construction. In the best 
case the overlap cures some small remaining O(a) cut-off effects in hadron spec- 
troscopy measurements, without creating larger 0(a 2 ) scaling violations. The 
hadron spectroscopy measurements performed in [E5l, which do not only contain 



simulations with £> par , but also include simulations with an order 3 Legendre 
polynomial expansion of the overlap Dirac operator from eq. (|5.3|) with $D par 
and -Rp ar at a lattice spacing of 0.155 fm (f3 = 3.0), may be interpreted this 
way. It is, however, very risky to make such a conclusion with so little data. The 
results on the dependence of the vector meson mass with respect to the pseudo- 
scalar mass show that there are larger discrepancies between the measurements 
with $D par and the corresponding overlap construction than one might expect. 
Furthermore, the energy-momentum relation for the overlap construction indi- 
cates cut-off effects. Even though the differences cancel in the Edinburgh plot, 
the data indicates different cut-off effects for £> par and the corresponding over- 
lap construction. To conclude this section we compare the physical branch of 
the eigenvalue spectrum of the blocked Dirac operator obtained from the RGT 
equations (|3~l| ) and ( |3.19D (Dirac operator II in the nomenclature of Chapter 0) 
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with the overlap construction with D par in Figure |5.4| . It shows that even if the 
overlap construction shifts the eigenvalues to the circle, there remains a visible 
difference between the two spectra, which indicates that the cut-off effects for 
the two Dirac operators might, indeed, be different. The question how large 
the cut-off effects are for the different Dirac operators discussed in this section 
has surely to be answered in future studies. 
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Figure 5.3: Breaking of Ginsparg- Wilson relation Agw(N) (top) and the 
breaking of normality A D (N) (bottom) on 10 4 lattices. Both quantities have 
almost exactly the same exponential fall-off, even though A Q (N) remains 
slightly smaller for all the orders of the Legendre approximation shown on 
this figure. The figure also shows that the exponential fall-off is rather 
governed by the condition number c, defined in eq. ( |5.13D , because the ex- 
ponent is roughly the same for the 3 different lattice spacings, where the 
lower boundary of the projected A* A eigenvalues is chosen such that the 
condition number c is approximately the same. 
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Figure 5.4: The physical branch of the eigenvalue spectra of the blocked 
Dirac operator (left) and the overlap Dirac operator with P pa r as kernel 
(right) are shown. One observes that the eigenvalues of the blocked oper- 
ator are shifted away slightly more from the origin than the eigenvalues of 
the overlap Dirac operator. This indicates that, even though the overlap 
construction cures the remaining breaking of chiral symmetry in 2? pa r> the 
overlap Dirac operator with T> p3ir might show different cut-off effects than 
the (approximate) FP Dirac operator, defined through the block transfor- 
mation in eqs. ( |3.1[ ) and ( |3.19j ). The lines on the figure are shown to make 
the different locations of the eigenvalues more evident. 



Chapter 6 

Topological Susceptibility 
and Local Chirality of 
Near-Zero Modes 



The zero modes as well as the low-lying modes of the Dirac operator have 
a particularly important influence on several important phenomena of QCD, 
because their contribution to the quark propagator gets very large in the chi- 
ral limit. Topological excitations like instantons might be responsible for the 
spontaneous breaking of chiral symmetry by building up a non-zero density of 
eigenvalues at zero in the infinite volume limit; a fact which is encoded in the 
Banks-Casher relation |141|. Instantons also provide an explanation for the large 
mass of the rf particle [ 142 1 . Related explanations for the mass of the rf predict 
a connection between the quenched topological susceptibility and the mass of 
the if jl43ip4ll . 

The question has recently been raised if it is possible to show that the 
near-zero modes are dominated by instantons or whether instantons are not 
important to explain the non-perturbative properties of the QCD vacuum. 

Answers to questions related to gauge field topology and its influence on 
the low-lying modes of the fermions living on such a gauge background can be 
explored within lattice QCD. Lattice Dirac operators which satisfy the GW rela- 
tion provide an excellent framework to study topological excitations of the QCD 
vacuum. They allow a definition of the topological charge and the topological 
charge density on the lattice which is unambiguous in the continuum limit (see 
Section 1.2.2| ). Hence, GW fermions can be used to study the topological sus- 
ceptibility and properties of the low-lying modes of the Dirac operator on the 
lattice in a very clean way. 

In the following we present some facts about the topological susceptibility. 
Then we will collect the data for the topological susceptibility which we obtained 
from different measurements with approximations to the overlap construction 
with V par as defined in Chapter |j| Finally, we present our results on the local 
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chirality of near-zero modes. 



6.1 Quenched Topological Susceptibility 



A mechanism that gives rise to the large mass of the rf was given in |]143| , |144 | 
and involves large vacuum fluctuations associated with confinement. In the 
leading order of the large N c expansion it allows to relate the mass of the ?/ to 
the strength of the topological fluctuations through the following relation 

xt = m,^Aii (fU) 

At V 2N f y J 

with 

m o = m %< + m rj ~~ , (6.2) 

where Q is the topological charge, \t the topological susceptibility and where 
(.) denotes the expectation value. Using the physical values for m v i, m^, ttik 
and fa one obtains \t ~ (180 MeV) 4 and lattice measurements typically yield 



Xt ~ (200 MeV) 4 |145| . Note that on the left hand side of eq. (fy) the 
quenched topological susceptibility \t enters as this equation has been derived 
in the large N c limit; for practical applications, however, one assumes that the 
physical value N c = 3 is already close to N c = oo. This assumption which is 
used quite frequently, clearly needs further testing, even though some results 
point into the direction that this approximation is reasonable, e.g. ||14rj[| . 

The topological susceptibility in quenched QCD differs substantially from 
the topological susceptibility in full QCD which vanishes in the chiral limit. 
The reason for this can be understood quite easily. After integrating out the 
fermion fields in the QCD partition function one is left with the determinant 
of the Dirac operator with mass m. As this mass goes to zero in the chiral 
limit the determinant vanishes, if the Dirac operator has an eigenmode at m, 
i.e. a zero mode in the chiral limit, and therefore configurations with a non-zero 
topological charge will be suppressed. In contrast to this the determinant in 
quenched QCD is set to 1 for purely technical reasons and hence there is no 
suppression of the zero modes in the quenched case. 

6.1.1 Lattice Determination of the Topological Charge 

The determination of the topological charge on the lattice has a long history and 
several different techniques have been developed in this time. The techniques 
vary from determinations from purely gluonic quantities like lattice discretiza- 



tions of the F^ U F^ U operator | 147 | to determinations from fermionic observables 



as the method by Smit and Vink [148], which in a certain way tries to make 



use of remnants of the Atiyah-Singer index theorem |74j]. While the traditional 
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gluonic methods run into trouble because the discretizations of F^ V F^ V typi- 
cally are very sensitive to ultraviolet fluctuations and require either cooling or 
smearing, the fermionic methods have the problem that the Atiyah-Singer index 
theorem does not hold for non-chiral lattice fermions. A solution is clearly pro- 
vided by a chiral formulation of lattice fermions as they satisfy the Atiyah-Singer 
index theorem also at finite cut-off and therefore a definition of the index of a 
gauge configuration can be given, even though the result may vary from Dirac 
operator to Dirac operator [42]. In the continuum limit, however, the results 
will ultimately agree. The FP Dirac operator offers, in principle, an even better 
solution, because the result for the fermionic charge agrees with the FP gauge 
charge at finite cut-off p2|| . While the discretizations of F^ V F^ V have to deal 
with large ultraviolet fluctuations even at small lattice spacing, the fermionic 
determination using chiral fermions and the chirality of the zero modes, i.e. the 
index theorem, does not have any problems with ultraviolet fluctuations because 
the zero modes are completely unaffected by the fluctuations at the cut-off. In 
the following we explain the method we use to determine the topological charge 
in our test studies. 



Counting Zero Modes and their Chirality 

The simplest way to make use of the lattice index theorem to determine the 
topological charge is to calculate the zero modes and to determine their chirality. 
In real simulations this actually means that one has to calculate the chirality 
of the small real modes, because in numerical simulations one cannot have an 
exact GW fermion and therefore the "zero modes" will scatter a little bit and 
become real modes \ip) with eigenvalues with a small modulus. Furthermore, 
their chirality, i.e. (ip\j5\ip}, will not be exactly ±1, but e.g. something like 
±0.7 to ±1 for the small real modes of V pav . In rare cases the real modes of 
D par can have an eigenvalue on the order of 0.3 and then also (V'lTslV') can 
be on the order ±0.2 to ±0.4. In our definition of the index of T> p&v we will 
sum up the chiralities of the small real modes with eigenvalues smaller than 
approximately 0.5 by simply counting the signs, independent of the magnitude 



°f ("0IT5 I V 7 ) ■ As we showed in Figure 14^ the real modes are quite stable, 
because their number and chirality does change only rarely when one compares 
the results from P par and the corresponding overlap construction with V pav . 
This already gives a hint that the level to which chiral symmetry is required to 
determine the number of small modes and their chirality is not very high. Using 
a moderate approximation of the overlap construction the small real eigenvalues 
will approach the origin and their chirality gets clearly closer to ±1, making the 
determination of the index by counting the modes less ambiguous than with 
P par . In order to calculate the low-lying eigenmodes we use again the implicitly 
restarted Arnoldi solver [pl3| , pT4|| . 
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6.1.2 Determination of the Quenched Topological Suscepti- 
bility 

In this determination of the quenched topological susceptibility, where we use 
the method described above to extract the topological charge of a gauge con- 
figuration, we collect data from various measurements that we have done with 
different approximations to the overlap with P par (see Chapter ||) and therefore 
the accuracy in the approximations varies substantially. In order to assure that 
for the counting of small real eigenmodes and their chirality there are essen- 
tially no differences between approximations with lower accuracy and such with 
higher accuracy we performed several tests, where we compare the outcome 
of the mode counting on samples of up to 100 configurations. The tests all 
lead to the same conclusion that no difference in the counting of the modes 
and their chirality was found between the different approximations we used. It 
shows that in the volumes we use it is actually enough to treat the 75 — 200 
lowest A^A eigenmodes - the eigenmodes which are responsible for the major 
part of the breaking of chirality - exactly and leave the rest of the "D par as it is, 
i.e. no Legendre expansion is used for the approximation for the square root in 
eq. ( 5.14| ). This shows that the degree to which chirality is needed to measure 



the topological charge with this definition is actually not very high. Obviously, 
the number of A^A eigenmodes which have to be treated exactly in order to 
provide enough accurate chiral symmetry depends on the gauge coupling and 
the volume. 

We determine the quenched topological susceptibility at several values of 
the gauge coupling (3, which lie in a range between 3.2 < j3 < 2.3, and on 



different lattice sizes. We give the details in Table p.l| . The scale is determined 
from the interpolating formula in |]96|| . The results are also shown in Figure 6.1. 
Notice that these measurements are pushed to the extreme: the resolution lies 
in between a = 0.13 fm and a = 0.36 fm, i.e. very coarse lattices are included 
and some of the measurements are on 4 4 and 6 4 lattices. Further, there does not 
exist, not even approximately, a unique way to fix r$ (and so a) on such coarse 



lattices. We believe, however, that Figure 6.1 contains interesting information 



6.1.3 Discussion 



The results shown in Figure 16^ give clearly the impression that the topological 
susceptibility measured with the overlap construction with D par has substantial 
cut-off effects. But one has to keep several things in mind. The decrease 
towards coarse lattice spacings is expected as the lattice gets too coarse to 
resolve topological objects, i.e. instantons fall through the lattice, and the results 
for the very coarse lattice spacings are only given for completeness. Obviously 
more important is the decrease towards small lattice spacings, which really 
indicates a relevant cut-off effect. It can be understood in different ways. The 
measurement at the smallest lattice spacing, which with a = 0.13 fm is still not 
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Figure 6.1: The quenched topological susceptibility xt measured with the 
overlap construction with V par . The open squares are used to denote in a 
clearer way that the corresponding symbol is shifted horizontally in order to 
display measurements done at the same scale, but with either different accu- 
racy of the overlap construction or on different volumes. The decrease of \t 
towards coarse lattices (large a) is an artifact due to the fact that instantons 
fall through the lattice, i.e. they can no longer be resolved. The decrease 
towards small lattice spacings is much more important as it indicates clearly 
a cut-off effect. 
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V/a A 


ro/a 


La [fm] 


^-conf 


(Q 2 ) 


4 

X* r 


N 


2.300 


4 4 


1.386(5) 


1.443(5) 


1000 


3.53(16) 


0.0508(23) 


5 


2.361 


4 4 


1.500(5) 


1.333(4) 


1000 


2.67(13) 


0.0528(25) 


2 


2.361 


6 4 


1.500(5) 


2.000(7) 


954 


13.28(62) 


0.0518(23) 


2 


2.400 


4 4 


1.576(6) 


1.269(5) 


2000 


2.51(8) 


0.0606(20) 


4 


2.500 


4 4 


1.787(6) 


1.119(4) 


1000 


1.87(9) 


0.0743(37) 


4 


2.680 


6 4 


2.237(7) 


1.341(4) 


1000 


4.24(19) 


0.0820(37) 





2.680 


9 4 


2.237(7) 


2.012(6) 


100 


24.46(3.37) 


0.0933(128) 


2 


2.927 


8 4 


2.969(14) 


1.347(6) 


990 


4.14(20) 


0.0787(40) 





3.000 


8 4 


3.197(16) 


1.251(6) 


600 


2.55(16) 


0.0650(41) 


4 


3.000 


8 4 


3.197(16) 


1.251(6) 


340 


2.54(24) 


0.0648(61) 





3.000 


10 4 


3.197(16) 


1.564(8) 


192 


6.18(65) 


0.0646(68) 


4 


3.200 


10 4 


3.947(26) 


1.267(8) 


200 


2.53(27) 


0.0614(76) 


2 



Table 6.1: In this table the different measurements of the quenched topo- 
logical susceptibility are collected. It shows the different gauge couplings 
j3, volumes V, the scale ro/a, the length of the box in fm (the error of the 
matching of ro to a physical scale is not included), the number of configura- 
tions n con f , the measured expectation value of the square of the topological 
charge (Q 2 ), the topological susceptibility measured in units of ro and fi- 
nally N gives the order of the Legendre polynomial approximation used for 
the overlap construction with 2? pa r kernel. Where no Legendre expansion is 
used, 75 to 200 A'A eigenmodes are treated exactly. 



very fine, is not very accurate and the picture might change with more statistics. 
Furthermore, the determination of the scale enters with (ro/a) 4 which makes 
that systematic errors become quite relevant. Clearly, the finite volume effects 
are not under control in this compilation of data and we can not exclude that 
the points at j3 = 3.2 and j3 = 3.0 are affected more by finite volume effects 
than the points at (3 = 2.927 and j3 = 2.68, which are measured on a larger 
volume. However, the finite volume effects are not expected to be large when 
the size of the box is 1.2 fm or more. Finally, we give a comparison to other 
recent measurements [ 149 , 150 1 of the topological susceptibility. Our results do 
not show any obvious discrepancy with the other data as the lattice spacing 
gets small, however, neither our data nor the data of the other groups are as 
accurate that a more precise statement can be made at this point. 



6.2 Local Chirality of Near-Zero Modes 

Exact zero modes of the Dirac operator tell us, via the index theorem, about 
the topological charge Q of the background gauge configuration. However, the 
exact zero modes alone cannot break chiral symmetry spontaneously. According 
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Figure 6.2: The quenched topological susceptibility obtained from different 
recent determinations. In this figure only the physically interesting part 
towards smaller lattice spacings is shown. We compare our results (FP), 
where we have combined the data from different determinations at the same 
lattice spacing in this figure, to the very recent determination in |150 |, which 
uses a truncated solution of the GW relation and counts the chirality of the 
small real modes (CHIR). The other two measurements are done by the CP- 
PACS collaboration and use traditional cooling techniques with the Iwasaki 
and the Liischer-Weisz action [14£], 
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to the Banks-Casher relation [141] (ipip) = —np(0) ^ 0, the Dirac operator 



must build up a finite density of near-zero modes, which does not vanish as 

V -» oo. 

One mechanism which explains the formation of near-zero modes involves 
instantons. Consider a gauge configuration containing one instanton and one 
anti-instanton. If the instanton and anti-instanton are separated by a large 
distance, the Dirac operator has a pair of complex eigenvalues lying close to 
0. The farther the instanton and anti-instanton are separated from each other, 
the closer the complex eigenvalue pair moves to the origin. If the instanton and 
anti-instanton are brought closer together, the complex eigenvalue pair moves 
away from the origin and disappears into the bulk of the eigenvalue spectrum. 
If the gauge configurations contain many instantons and anti-instantons, this 
could produce a non-zero density of near-zero modes, giving p(0) ^ in the 
infinite volume limit. 

The question has recently been raised if it is possible to show that the near- 
zero modes are dominated by instantons. From instanton physics, it is expected 
that the modes are highly localized where the instantons and anti-instantons sit. 
If this is so, then in these regions the modes should be close to chiral i.e. mostly 
either left- or right-handed, depending on whether it is sitting on an instanton 
or anti-instanton. In [151], the authors defined a measure of local chirality at 
lattice site x by 



tan 



7T 



4 



V'lVlOe) 



t(i +*(*)) = j ;vt ■ ( 6 - 3 ) 



An exact zero mode is purely either left- or right-handed, giving X(x) = ±1 
at all lattice sites x. If a near-zero mode is localized around instanton-anti- 
instanton lumps, then X(x) should be close to ±1 for the sites x where the 
probability density rp'tp(x) is largest. 

In the original paper |151[ , many near-zero modes of the Wilson operator 
D]y were analyzed for many gauge configurations and the finding was that in 
the regions where the modes are localized, the distribution for X{x) is peaked 
around and the modes do not display local chirality. This led to the conclusion 
that the near-zero modes are not dominated by instantons. Since then, several 
other groups have found the opposite conclusion ||117| , p.52 - 154 |, using Dirac 



operators with much better chiral symmetry than Dy/ (or even just an alterna- 
tive definition of a complete basis for the non-normal operator D\y). They find 
the distribution of X{x) is double-peaked with maxima at large positive and 
negative values of X, indicating that the modes are locally chiral. 

6.2.1 Results 

In our study we analyze the 10 smallest near-zero modes of the overlap operator 
with ©p ar as kernel for 60 different 10 4 gauge configurations at /3 = 3.2. We 
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use a Legendre polynomial of order 2 to approximate the operator (A^A) -1 / 2 
in overlap formula, with the 10 smallest A^A being projected out and treated 
exactly (see Chapter |] for details). The eigenvalues A and eigenvectors ip are 



found using the implicitly restarted Arnoldi method [112,114|. In Figure p^3| , 
we plot the distribution P(X) of the measure of local chirality X at the lattice 
sites where the density ip>ip(x) of a mode is largest. The three distributions 
correspond to taking, for each mode, 1%, 5% and 10% of all lattice sites that 
have the largest density t/)*ijj(x). We do not include exact zero modes, for which 
X(x) = ±1 at all lattice sites. We see a clear double-peaked distribution, whose 
maxima are farther from zero if we only include the sites where the modes 
are most localized. The maxima are not at X = ±1 as the modes are not 
exactly chiral. We find the same conclusion as [117, 152- 154| : where the near- 
zero modes are most localized, they are also very chiral. Recent work [ 155 1 
has shown that, using near-zero modes of the Dirac operator, the gauge field 
strength appears to contain lumps which are mostly either self-dual or anti-self- 
dual. This behavior is consistent with the picture of instanton-dominance of the 
near-zero modes, however it is no conclusive evidence that instantons are the 



driving mechanism for chiral symmetry breaking [156, 157 1 
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Figure 6.3: The distribution P(X) for the top 1%,5%,10% lattice sites x 
with the largest ip'ip(x). 



Chapter 7 

Determination of the 
Low-Energy Constant S 



It is the common expectation that for QCD with a number Nf > 2 of mass- 
less quark flavors, the chiral symmetry is spontaneously broken by a non-zero 
expectation value for the chiral condensate (^0). Chiral perturbation theory 
(xPT), which is based on this assumption, is an excellent description of many 
low-energy QCD phenomena Q. However, it is only possible via lattice QCD 
to test from first principles if the symmetry is spontaneously broken. 

The leading order effective theory of xPT contains the low-energy constants 
f w and S. In full QCD the Gell-Mann-Oakes-Renner relation (GMOR) 



r2 2 



(7.1) 



becomes exact in the m — > chiral limit and the chiral condensate (ipip) = 
(uu) = (dd) = . . . , defined at zero quark mass, is equal to -E. The low energy 
constant £ depends on the number of massless flavors Nf. 

In quenched QCD, the relation (^ip) = — S and eq. (7.1) receive corrections 
even in the chiral limit [ 15£ — 1160| ] . Actually, (tpi^) and m^/m are not even 
defined in the chiral limit due to diverging quenched chiral logarithms. In the 
case of rn^/m the logarithms, which are characterized by the parameter 5, are 
connected to the mass of the flavor singlet particle 7/Q, which is an additional 



Goldstone boson in QxPT [15£ , 159 [. The mass of the rj is an additional, third 
leading order low-energy constant mo of QxPT, which is not present in the 
normal xPT Lagrangian. The expectation that m^/m diverges in QxPT in 
the chiral limit is confirmed by different measurements [161-163], including our 



hadron spectroscopy measurements with $D par |65]. Similarly, the divergence of 
{ipij}) in the infinite volume limit seems to be confirmed in numerical studies 



of the Banks-Casher relation ||164| . On the other hand, it is possible to study 
and determine the low-energy constant E in the quenched theory. Under the 

1 The mass of the rj' in the quenched QCD Lagrangian does not correspond to the mass 
of the physical rj' , because the physical rj is mixed with the rj. 
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assumption that £(-/V/) is a smooth function of Nf and £(iV/ = 0) is close to 
E(iVy = 3) (which is the standard assumption when quenched results are used 
to estimate full QCD quantities) we get an estimate for the chiral condensate 
-{i>ip){Nf = 3) = E(Nf = 3) « E(Nf = 0), even though it is hard to give a 
systematic error for this estimate. Typically, the quenching error are said to be 
of order 10-20% @. 

On the lattice the calculation of the low-energy constant £ of (quenched) 
chiral perturbation theory requires good control over the chiral properties of the 
theory and therefore it has notoriously been difficult to measure it directly. Still, 
early attempts with staggered fermions have been tried already in 1985 [165] and 
the calculation can also be performed with other traditional discretizations, like 
0(a) improved Wilson fermions using Ward identities and the GMOR relation 
eq. (|7.l| ) |166| , 167]. However, the calculation is more transparent and cleaner 
in the framework of GW fermions due to many reasons. With Wilson fermions 
one has to deal with a power divergent subtraction that has to be determined 
numerically |166(| . Furthermore, the occurrence of exceptional configurations 
makes it difficult to go to quark masses which are small enough to provide an 
accurate estimate for £ in the chiral limit. In the case of GW fermions chiral 
symmetry is already realized on the lattice [21]. This makes it possible to define 
scalar and pseudoscalar densities which transform under chiral transformations 
exactly like the corresponding continuum densities []52]]. Using such a scalar 
density the power divergent subtraction, which has to be performed on the bare 
condensate, is exactly known and can be done without problems in numerical 
simulations |39]. A further advantage of GW fermions is that one does not 
encounter any exceptional configurations in the small quark mass regime. These 
facts altogether allow a much better control over the calculations. 

There are different techniques to calculate £ with GW fermions. One pos- 
sibility makes use of the GMOR relation and data from pseudoscalar spec- 
troscopy | 168 |. Another possibility is to use the distribution of the lowest-lying 
eigenmodes of the massless Dirac operator, because these distribution can be 



described in the framework of random matrix theory (RMT) |169j and provide a 
completely different way to extract £. Actually, the distributions of the low-lying 
eigenmodes contain a lot more of information than the value of £ as explained 
in [ |170fl , but in leading order they are essentially described by £. A technique 
which is closely related to the determination of £ from the eigenvalue distri- 
butions is the so called finite-volume scaling technique, which has been worked 
out in [16C], and has been used for the first time in [130]. In our calculation of 
£ we use a slightly modified version of this technique. 

All the 3 different methods provide a bare value of £ that subsequently has 
to be renormalized in order to give the phenomenologically relevant predictions 
£ and £j^g(2GeV), respectively. Again, there are different possibilities how 
the scalar renormalization factor Z$ can be calculated. It is known that lattice 
perturbation theory converges only slowly in many cases and is therefore very 
difficult to control in the range of gauge couplings where typical simulations 
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are performed [171]. The way to circumvent this problem is the use a non- 



perturbative renormalization scheme. There are essentially two different non- 
perturbative renormalization schemes used in the present day lattice simulations. 



One scheme is the RI/MOM technique []172| , which uses the fact that QCD 
can be treated perturbatively at large momentum scales, in order to match 
amputated quark Green functions to a perturbative scheme (e.g. MS), while 
the other scheme, the Schrodinger functional technique, is based on the fact 



that QCD in a very small volume can be treated perturbatively |]171][ . The 
size L of the box acts as a reference scale and can e.g. be connected to the 
normalization mass in the MS scheme using perturbative renormalization group 



arguments. We will follow the renormalization scheme proposed in [173|, which 
is based - even though not directly - on the Schrodinger functional technique. 

But, before we start to explain our determination of S, let us mention that 
more details of the calculation, in particular about the stochastic evaluation of 
the trace of the subtracted quark propagator, are given in Appendix 

7.1 Determination of the Bare E 

One possibility to determine £ is to study the chiral condensate in a fixed 
topological sector with charge Q in a finite volume V at finite quark mass 
M [|160| . The volume and the quark mass are chosen such that the finite size 
effects are dominated by the pions with zero momentum. In this situation the 
partition function simplifies considerably and can even be evaluated analytically 
in certain cases [ 160| , 174 1. Using (Q)xPT, or random matrix theory (RMT), the 



fermion condensate at finite volume and quark mass has been calculated in the 
continuum, both for full and quenched QCD. The quenched QCD condensate 
{'ipip)M,v,Q can be written in terms of the derivative of the logarithm of the 
partition function Zq in fixed topology 

d 

-(iH>)m,V,Q = QM lnZ Q> ( 7 ' 2 ) 

with 

Z Q = f dC/ Sdet(rJ Q )exp[ySReStr[MU ]], (7.3) 
JG1(1|1) 

where Sdet and Str are the supersymmetric generalizations of the determinant 
and the trace, respectively. The integration is over the graded Lie group Gl(l|l). 
For details of the supersymmetric formulation of QxPT see [158- |160| , |170| , |175 , 



176]. The analytic expression for this partition function is known [160] and 



therefore the condensate in fixed topology at mass M is given by 

-{H>)m,v,q = MVY?[I u {z)K u {z) + I lQl+1 (z)K lQhl (z)] + & (7.4) 
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where Iq and Kq are modified Bessel functions, z = MT.V and the low-energy 
constant £ is the quantity we wish to measure. By measuring (tpip) m,v,Q in 
different topological sectors, at different masses and volumes, the continuum 
prediction of the M and V dependence can be used to extract S. 

On the lattice the subtracted condensate in a fixed topological sector in 
a volume V and with bare mass m can be calculated by measuring the trace 



V 



Tr' 



2>- A (m) -R 



(7.5) 



with 



®(m) 



m 



©(0) +m(2R) 



-i 



(7.6) 



and where the primed trace is without the divergent contribution of the zero 
modes and (.)q is the expectation value on the ensemble of gauge fields with 
topological charge Q. Because we want to make use of the renormalization 
scheme proposed by Hernandez et al. in [173], which is based on pseudoscalar 
spectroscopy, we have to modify the strategy of the measurements performed 
n [116, 13Cf] slightly. The reason for this is that in our hadron spectroscopy 



study [5g] a mass definition m is used which differs from the usual definition 
M by a factor 2R (see e.g. eq. (ffTq 1 )). Instead of the trace in eq. dTT 



(^(m)2R)- 1 - -S)(0)S> -1 (m) 



we use 



(7.7) 



The scalar density S is the same as introduced in Chapter || and its properties 
are discussed in more details in [14C], i.e. 



S = i;((2R)- 1 --®(0))i>. 



(7.* 



The partition function Zq is given by 



Z Q = / DUqDtPD^ exp[-(^ ©(0)V> + mS) - A g ] 



9} ■ 



(7.9) 



with the gauge action A g . Furthermore, the path integral is restricted to gauge 
fields with topological charge Q. For masses m smaller than the modulus of the 
smallest eigenvalue A m i n of the massless Dirac operator D(0), i.e m 2 < |A m i n | 2 , 
the trace in eq. (|7.7|) reduces to 



m,V,( 



v 



(f){m)2R)' 



(7.10) 



Q 



7.1. Determination of the Bare £ 



97 



up to terms of 0{m 2 ), as one can show easily. Because the numerical results 
show to be less sensitive to a residual breaking of chiral symmetry with the defi- 
nition from eq. (7.10) than with eq. (|7.7[), we use the definition from eq. (|7.10| ) 
for our measurements, which clearly agrees with the definition in eq. (|7.7[) in 
the chiral limit. The correlation function (S) m y : Q differs from (^V')mVq by a 
factor (2R) . Together with the scalar renormalization factor Zg, which for 
chiral actions is given by Z$ = 1/Z m , however, the additional factor drops out 
in the physical prediction, because the additional factors from the mass defini- 
tion and (S) m y ! Q cancel. Let us be a bit more specific about what is meant 
with this: In our case the eigenvalue spectrum of the operator 2R is roughly 
bounded by 2/k and 4/k with k = 3.45 for our parametrization and therefore 
its eigenvalues are close to 1 and no exceptionally small or large eigenvalues can 
occur, which in a simulation might spoil the fact that the two additional factors 
of (2i?) _1 and 2R in (S) m y,Q and m, respectively, cancel to a high degree. 



7.1.1 Details of the Simulation 



We measure the trace in eq. ( 7.1C ) stochastically using random Z{2) vectors 
[177]. In order to measure the condensate at very small quark mass, the remnant 
explicit chiral symmetry breaking must be very small, so we use the overlap 
operator with the kernel £> par , as discussed in Chapter ||. Due to the \Q\ 
chiral zero modes of the Dirac operator, the quenched condensate contains a 
term \Q\/MV which diverges as the mass becomes small. The \Q\ zero modes 
have the same chirality and their contribution to the condensate is removed by 



measuring the trace TV in the chiral sector opposite to the zero modes [13C , 178 ], 
i.e. if the \Q\ modes have chirality +, the Z(2) vectors used to measure the 
trace are chosen to have chirality — . To determine (S) m y ! q, the stochastic 
trace is doubled to include both chiral sectors. 

We measure the condensate in volumes 8 4 and 10 4 at (3 = 3.2. We use 
10 random Z(2) vectors to measure the trace for each configuration and a 
BiCGstab multi-mass solver to invert the Dirac operator at all 12 masses in the 
range ma = 0.0001,... ,0.2048 simultaneously. In the overlap operator, we 
approximate (A* 1 A)' 1 / 2 with Legendre polynomials of order 7 and 10 for volumes 
8 4 and 10 4 , respectively. We project out the 10 smallest A eigenvalues, which 
are treated exactly. This altogether gives sufficiently precise chiral symmetry: 
increasing the polynomial order further, the relative change in the estimate of 



(("5) m ,v,Q0 3 )/(™) is < O(10~ 4 ). This fact is shown in Figure 7T, where the 
relative difference between the (nearly) exact result with a Legendre expansion 
of order 12 and lower order expansions with order N at ma = 10~ 4 for the 
quantity 

A(N) _ (m (N) (rn)2R)-i ^|(<p( 12 ) \m)2R)^ - \\<j>) 

m& 12 Xm)2R)-^ - \\<i>) ' 
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L 


order 


\Q\ 


-^conf 


8 


7 


l 


154 






2 


41 


10 


10 


1 


53 






2 


43 



Table 7.1: Polynomial order of (^A) l / 2 approximation and statistics for 
the measurement of S. 

where is a Z{2) random vector, is shown for several 8 4 configurations. Our 
statistics are given in Table |7~I|. In Figure |7~2| , we show the distribution of the 
estimates of ({S) m y t Qa?) / (ma) obtained from the different Z(2) random vec- 
tors in \Q\ = 1 topological sector on the 8 4 and 10 4 configurations, respectively. 
One observes that the distributions have quite a long tail towards large values, 
especially on the 10 4 lattice. 



7.1.2 Results 



Let us first explain how the bare value of £ can be extracted from eq. (|7.4|) 



and the measurements of the trace in eq. ( 7.10 ). The bare quark condensate at 
finite quark mass contains a ~ m/a 2 cut-off effect. As the quark mass m — > 0, 
one can rewrite eq. (\IA\) such that 



(^)M b V,Q 2 {S) m ,V,Q d 

— j^ = ~ a ^ l — = m + ^ ( } 

where c\ is an unknown coefficient, which has to be fitted, and a denotes the 
ratio between the different mass definitions: a = m/M . As the coefficient c\ 
comes from ultraviolet fluctuations, it is natural to assume that it is independent 
of the topological charge Q. The contribution of the \Q\ zero modes is removed 
by taking the trace as described above, and there is no artifact l/(ma 3 ) in 



eq. (7.12) due to the fact that the condensate is defined as the expectation value 



of a scalar operator that transforms covariantly under chiral transformations 



and has no mixing with the unit operator |39|]. From Figure [7.3| , we see that 
— ((S) m y t Qa 3 )/(ma) reaches a plateau at very small quark mass. In Figure |7~4| , 
we plot the value of — ((S) m y : Qa 3 )/(ma) at ma = 10~ 4 versus (Va~ A )/2\Q\, 
which we fit to the form of eq. ( |7.12j ). From the slope, we extract the bare 
low-energy constant as Sa 3 = 4.42(36) x 10~ 3 with \ 2 /df = 3.55. We can 
convert this result using the Sommer parameter (ro/a = 3.943(60) at /3 = 3.2 
from the interpolating formula in (H), giving Erg = 0.271(22)(12), where the 
first error is statistical and the second the uncertainty in the scale. 
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Figure 7.1: In this figure the dependence of the relative difference A(iV), 
defined in eq. ( 7.11| ), between the (nearly) exact overlap and lower order 
Legendre approximations with order N is shown for several different 8 
configurations at f3 = 3.2. 
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Figure 7.2: Normalized distribution p{s) of the estimate of s = 
{(S}m,v,Qa 3 ) I '{ma) at ma = 10~ 4 for the different Z{2) random vectors and 
configurations used in the measurements is shown for the Q = ±1 sector on 
8 4 and 10 configurations. 
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Figure 7.3: The estimated ({S) m y^QO?) / {ma) for different volumes, masses 
and topological charges. 
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Figure 7.4: ((S} m y t Qa 3 ) / (ma) for the smallest value of the mass ma = 10~ 4 
plotted versus (Va~ 4 )/(2\Q\). The fit used to determine the slope is shown 
by the line. 
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7.2 Renormalization of E 



In order to turn this bare result into the renormalized low-energy constant E 
we need the scalar renormalization factor Zg. In the present test study we ob- 
tained Zs combining the continuum extrapolated renormalization group invari- 



ant (RGI) quark mass of the ALPHA collaboration [17S] with our spectroscopy 
data. This method has been suggested recently by Hernandez et al. [ 173 1 . Even 
though the method is based on the Schrodinger functional, the special boundary 
conditions of Schrodinger functional make a direct application to GW fermions 
difficult, if not impossible, and therefore improved Wilson fermions were used 
in an intermediate step. In the following we give a short description of the 



calculations done in [ 173 1 and apply the result to our data. 



7.2.1 Short Introduction to the Technique 

The idea of this technique is to connect the bare quark mass m in the given 
discretization to the renormalization group invariant (RGI) quark mass M 



M = Z M (g )m(g Q ) 



(7.13) 



and to use the fact that for lattice regularizations that preserve chiral symmetry 
the renormalization factor of the scalar density is the inverse of that for the 
quark mass (see, e.g. [ |180| ) 

Zp{go) = Z s (g ) 



In order to fix Zm(9o) one can rewrite the ratio M/m(go) as 

M = M m w (g' ) 
m(g ) mwig'o) m (fo) 



^ {r m w )(g' ) 
(r Q m)(g ) 



(7.14) 

(7.15) 
(7.16) 



where g' ^ go in general, and the hadronic radius r$ is used to set the scale. The 
factor Zjjjr has been computed in |181| for a large range of couplings. The ratio 
M/m(go) is then obtained by determining (rnm) and (romw) at a reference 
value of x re f of a physical observable, say x re f = (romp) 2 , where mp is the 
mass of the pseudoscalar meson. In order to get rid of discretization errors the 
combination Z^(g )(romw)(9o) was calculated in the continuum limit, i.e. 

U M = Hm Z^(g' )(romw)(9o) , (7-17) 



leading to the following results ]173 | 

'0.181(6), 

U M = 



x Te f = 1.5736, 
0.349(9), rc ref = 3.0, 
w 0.580(12), x Te{ = 5.0. 



(7.18) 
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By combining ( 7.13j ), ( [7.1 5| ) and (7.17) one obtains 

1 



Zh[{go) = U M 



(r m) 



(r m P ) 2 =x rcf 



(7.19) 



where now all reference to the bare coupling g' Q and the use of the intermediate 
Wilson fermions has disappeared. Apart from cut-off effects, Zm{9o) is expected 
to be independent of the reference point x re f . 



7.2.2 Results 

Using the results of the hadron spectroscopy with the overlap with © par at 
P = 3.2 |^5|, we are able to calculate the scalar renormalization factor Zs- 
The pseudoscalar spectroscopy data is obtained from 32 9 3 x 24 lattices. Be- 
cause of the small statistics and the rather small volume the determination of 
the renormalization constants at f3 = 3.2 has a rather large error. For our 
determination we combine data from the pseudoscalar - scalar (PS) and the 
axial-vector (A) correlator measurements of the pseudoscalar mass. The reason 
for this is that the PS correlator gets contributions from the scalar at high bare 
quark masses, while the A correlator is sensitive to the topological finite size 
effects caused by the zero modes at small bare quark masses [182]. Using this 
combined data, we fit a second order polynomial through the data points in the 
region m = 0.01, . . . , 0.17 as shown in Figure |7~5|. With these fits we determine 



the bare quark at which the pseudoscalar takes on the three reference values 
^ref = (r m P ) 2 = 1.5736,3.0,5.0. We obtain 



(r m)| 



(r m P ) 2 =x Ici 



0.141(24), x re f = 1.5736, 
0.271(36), x rcf = 3.0, 
[0.447(24), x re{ = 5.0, 



(7.20) 



where we use the hadronic scale ro/a = 3.943(60) from the interpolating for- 
mula in [^|] to set the scale. In the determination of the r^m values a small 
additive mass renormalization oA m = 0.005(4) is taken into account. This 
residual additive mass renormalization is due to the fact that in the determina- 
tion of the pseudoscalar spectrum a Legendre expansion of order 3 was used. 
This leaves place for a breaking of the GW relation which is large enough to 
give rise to a residual mass on the order of 10~ 3 . 

Using (|7.19|) , we can combine the results for the bare quark masses in 
( |7.20| ) with the universal factor Um from ( 7.18 ) and obtain the mass and scalar 
renormalization factors. The results are given in the Table |7~^ . 

The dependence on the value of x re f is weak and well covered by the sta- 
tistical uncertainty. Following the arguments in [172] we choose the value at 
x re f = 3.0 as our best estimate of the renormalization constant, and thus we ob- 
tain the following final result for the renormalization factor Zm{9o) = 1.29(12) 
and Zs(go) = 0.869(32), respectively. These renormalization factors can be 
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0.15 



a m 



Figure 7.5: The (mpa) 2 versus ma data used to determine the scalar renor- 
malization factor at f3 = 3.2. The solid line is a fit using a second order 
polynomial, while the dot-dashed lines indicate the three reference values 
for the pseudoscalar mass x re f = (romp) 2 = 1.5736,3.0,5.0 used to extract 
Zm- 



■^ref 


Zm(9o) 


z s (go) 


1.5736 
3.0 
5.0 


1.28(25) 
1.29(12) 
1.30(7) 


0.78(15) 
0.78(7) 
0.77(4) 



Table 7.2: Renormalization factors at 0=3.2 for the different reference values 
of the pseudoscalar mass. 
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converted into the MS scheme at the reference scale \i = 2 GeV using the 
results from 111791 



0.72076, / u = 2GeV, (7.21) 



^Ms(^) = Z rn{gO,V) 

M ~ Z M (g ) 



with an error of 1.5% due to the uncertainty in the quenched value of Ajjg. 
This leads to Zm(9o,^) = 0.93(9) and Zs(go,^) = 1.08(10), respectively. 

We can combine the scalar renormalization factor with the bare value of E, 
which leads to the final prediction for the physical value of E and Ej^g(2GeV). 
We get 

rgS = 0.210(17)(9)(20), (7.22) 



r^S M g(2GeV) = 0.291(24)(12)(28) , (7.23) 



with statistical, scale and renormalization errors, respectively. Combining all 
errors, but the error associated with the scale ambiguity in quenched QCD in 
quadrature and using ro = 0.5 fm, we find at (3 = 3.2 



E IB (2GeV) = (262 ± 12 MeV) 3 X ( °~ 5 'j^ ) • (7.25) 

These results are still subject to discretization errors and clearly it would be 
desirable to have data at smaller lattice spacings in order to make a continuum 
extrapolation. 

Using the GMOR relation from eq. ( |7.1| ), our result for E and the slope B 
in the vs. m plot, which we determine to aB = 2.90(12), we can make a 
consistency check of our measurements. We get af n = 78(4) x 10~ 3 , which 
corresponds to 

(121.5 ±5.9 MeV)x(g^l), (7.26) 

i.e. we get sufficiently close to the value from other determinations of in 
quenched QCD, e.g. [ |183| |, which are quite close to the experimental value f w = 
131 MeV, to be confident that no big systematic errors are induced through the 
low statistics or the small volumes we used in this study. 



7.3 Discussion 



The result of this determination of the low-energy constant E in quenched QCD 
is in very good agreement with other recent determinations with chiral fermions 
[ 115 , 168j 173H , as shown in Figure 7^. One should, however, mention that 
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Figure 7.6: Comparison of different recent determinations of rQ$]j^g(2GeV). 



these studies - including ours - are based on poor statistics and the simulations 
were done in small physical volumes, which means that all these studies have a 
certain test character. Furthermore, none of these studies contains a continuum 
extrapolation. The small volume is a problem because the RMT prediction of 
the distribution of the low-lying eigenmodes is only valid up to the Thouless 



energy [184, 185], which means that only a small fraction of the eigenmodes 
are described properly. This clearly leads to a systematic error. The problem of 
small statistics is enhanced through the fact the trace in eq. ( |T . 1 0| ) gets large 
contributions from the smallest eigenvalues of the Dirac operator. Hence, it is 
probable that the results are distorted by the very rare but sizeable contributions 



from very small eigenmodes and therefore we think that the statement in [185] 



that the distributions of the lowest-lying eigenvalues of the Dirac operator would 
provide an approach to the determination of E, which has smaller systematic 
errors, is most likely true. 



Chapter 8 

Conclusions and Prospects 



The parametrized FP Dirac operator 

In this work we have presented the construction and parametrization of a FP 
Dirac operator using a general hypercubic ansatz for the parametrization. We 
have shown how general Dirac operators can be constructed in an efficient 
way such that it is feasible to use them in numerical simulations. We have 
discussed the details of the parametrization procedure and have shown where 
the main difficulties lie. The properties of the resulting parametrization of 
the FP Dirac operator $D par have been investigated. We have found that the 
parametrization has clearly improved chiral behaviour with respect to the most 
frequently used lattice Dirac operator, the Wilson Dirac operator D\y. This 
fact is most evidently seen in the eigenvalue distribution of the operator A, 
which for a Dirac operator satisfying the Ginsparg-Wilson relation is equal to 
1. It deviates much less from 1 for $D par than for Dy/- Furthermore, we have 
found that the Krylow space methods used to calculate the eigenvalues and 
eigenvectors of the operator A' A, which are important for an efficient numerical 
implementation of the overlap construction, use substantially less operator times 
vector multiplications in the case of © par than for Dyy. We have found that 
the difference in the number of operator vector products can get so large that 
the determination of the low-lying eigenmodes of A^A actually is only very 
little more expensive for © par than for D\y, even though one matrix vector 
multiplication is O(40) times more expensive than for Dyy- We have shown 
that the additive mass renormalization and the fluctuations of the low-lying 
eigenmodes of the Dirac operator are small, such that hadron spectroscopy 
simulations at remarkably small quark masses can be performed. Because of 
the possibility to go to very small quark masses clear signs of quenched chiral 
logarithms have been found in our hadron spectroscopy, which is described in 
detail in the thesis of Simon Hauswirth [65]. In the same study we have found 
that the scaling properties of $D par are, indeed, clearly improved to Dy/ and 
0(a) scaling violations can not not be seen, even though more precise studies 
are needed to give a more quantitative statement. Finally, we have found that 
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the definition of the topological charge through Tr(75i? par £> par ) is problematic 
with the present parametrization as it deviates substantially from the index of 
©p ar , i.e. it does not satisfy the index theorem. 

Overlap construction with $D par 

We have shown that the overlap construction with £> par is clearly more local 
in the sense that its couplings fall off faster than the common construction 
with Dy/ . In contrast to the Wilson overlap there is no need to fine tune a 
mass parameter to get optimal convergence or optimal locality in the overlap 
expansion with £> par . We have found that polynomial approximations to the 
overlap construction with $D par converge very fast, i.e. that chiral symmetry 
up to machine precision can be reached with Legendre polynomial expansions 
of 0(20) for the lattice sizes and gauge couplings used in this work. It has 
shown that such a precision, however, is barely needed in most simulations. We 
have found that low order expansions can be used in many cases to get enough 
accuracy for the chiral properties of the Dirac operator. 

Physical Results 

We have used low order polynomial approximations to the overlap construction 
with © par to calculate the quenched topological susceptibility at various lattice 
spacings and volumes. The results point in the same direction as other recent 
results. However, our data is not yet accurate enough to allow for a controlled 
continuum extrapolation and shows cut-off effects. We also confirm the findings 
of various other groups about the local chirality of the near-zero modes of 
a Dirac operator, which might result from instanton dominance in the QCD 
vacuum. Finally, we have determined the low-energy constant £ of quenched 
chiral perturbation theory using a finite-volume scaling technique in a test study. 
Our result for the renormalized £ is 

^£^(2 GeV) = 0.291(24)(12)(28) at (3 = 3.2 , 

with statistical, scale and renormalization errors, respectively. This result does 
not include a continuum extrapolation and therefore is still subject to discretiza- 
tion errors. Converting to physical units by using the scale ro = 0.5 fm, which 
is somewhat problematic due to the scale ambiguity in quenched QCD, we get 

^ms(2 GeV) = (262 ± 12 MeV) 3 x (^f^) ' at = 3.2 , 

where all errors are combined in quadrature. This result is in good agreement 
with other recent studies of £ using chiral fermions. It should be mentioned 
that all these studies, including ours, are based on rather small statistics. Nev- 
ertheless, we have shown that with $D par we can use low order polynomial ap- 
proximations to the overlap construction to perform calculations where chiral 
symmetry is needed to high precision. 
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Prospects 

Let us start with the most obvious weakness of the present parametrization: 
The topological charge definition through Tr(75-R par £> par ). In the present pa- 
rametrization no special attention has been given to the pseudoscalar part of 
the parametrization and it is presently under study whether this weakness can 
be fixed with additional care in the parametrization process. If this is possible it 
would provide a numerically very cheap definition of the topological charge and 
moreover of the topological charge density. The overlap construction with $D par 
is suited to study the structure of topological excitations of the QCD vacuum, 
because there are indications that this overlap construction can resolve the un- 
derlying structure of the gauge fields better than the Wilson overlap. Using 
approximations to the overlap operator with £> par , or simply © par itself, further 
problems where chiral symmetry plays an essential role might be investigated, 
because the list of such problems is long. Let us mention just a few: De- 
termination of further low-energy constants of (quenched) chiral perturbation 
theory, determination of quark masses of the light quarks, the resolution of the 
Kaplan-Manohar ambiguity in the first non-leading order Lagrangian of chiral 
perturbation theory and the calculation of the pion scattering length. Using the 
available data from the present study on hadron spectroscopy some of these 
problems can already be tackled. The data contains information about the pion 
decay constant, the light quark masses and further information about the pa- 
rameter 5 describing the quenched chiral logarithms. In order to extract some of 
this data one might use e.g. the RI/MOM scheme to calculate renormalization 
factors accurately. In connection with this renormalization scheme, which is 
sensitive to momentum cut-off effects, the good energy-momentum dispersion 
relation of © par might be a big virtue. Furthermore, first applications of con- 
served chiral currents, which we proposed in [69, 14C], might be interesting and 
useful. 
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Appendix A 

Representation of the 
Clifford Algebra 



In this appendix we give the explicit representation of the 16 elements of the ba- 
sis of the Clifford algebra in the Weyl representation, as they are used throughout 
this work. In the current implementation of the programs for the parametrized 
FP Dirac operator the same convetions are used. 
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Pseudo-scalar 
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Appendix B 

The List of Factorized 
Contributions 



In this appendix we provide a factorization of the paths used in the program 
code of the parametrized FP Dirac operator. The correspondig reference paths 



and their properties are given in the Tables 2.1-2.5. 

B.l Definitions 

We introduce the operator of the parallel transport for direction \i 

(%) = U^n)8 n+ix>n , , (B.l) 

V / ran' 

and analogously for the opposite direction 

([/_„) = U^U - ttfSn-^n' ■ (B.2) 
V / nn' 

It is also useful to introduce the operator U(l) of the parallel transport along 
some path I = [hjfa, ■ ■ ■ , h] where li = ±1, . . . , ±4, by 

U{l) = U h U h ...iJ h - (B.3) 
In terms of gauge links this is 

lU{l)) =[U h (n)U h (n + h)...) 

"ra+r;,ra' > 

(B.4) 

V / ran' V / 

We define the plaquette products (i.e. the operators of parallel transport 
along a plaquette) in the following way: 

P h , h = U([h,l 2 ,-k,-h}), (B.5) 

where = ±1, . . . , ±4. Their hermitian conjugate is given by (-P^,^)^ = Pi 2 ,h ■ 
Reflections and permutations act on them in an obvious way. 
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We also define the staple products as 

S h , h = U([h,h,-h}). (B.6) 

We have (S h ^ = S h ,- h . 

To describe the shortest path to an offset we introduce the notation 

V h = U([h}), 

V hM = U{[hM}): 

V h M,h = Hh,hM), (B.7) 
v h,h,hM = U([h,h, h, k]) ■ 

Introduce the following linear combinations transforming in a simple way 
under reflections: 



(B.8) 



The signs in the superscript denote the parity for reflections of the /x, v axes, 
respectively. Hermitian conjugation acts as interchanging both upper and lower 
indices, permutations simply by jjl — > p^, v — > p v . 

It is also useful to denote combinations which are symmetric/antisymmetric 
with respect to interchanging the axes: 



p(++) 

fMU 


— P 


+ ^M, 


-u + P- 


-\x,v 


+ P- 


"M, 


— 1/ 


p(+-) 

J fMU 


= Pfi,u 


-P* 


-u + P- 


-fM,U 


- p. 


"M> 


— i/ 


p(-+) 


— P 


+ P», 


-u-P- 


-\x,v 


- P. 


"M> 


— i/ 


p(~ ) 

1 fMU 


— p 

— Mi 1 ' 




-u-P- 


-fM,U 


+ P- 


"M> 


— i/ 



^ S / m) = ^ +) + ^ +) = Pi% +) + h-c. , (B.9) 
For the staples we write 



p(as) _ p( — ) _ p( — ) _ p( — ) _ v. r 



o(",+) _ C _j_ c 

"fi ~~ ~ ^iM ' 

S^ = S v%ll -S- v ^, (B.10) 

C , ( I/ ' + ) C _l_ Q 

u — J u,~fM ~r >->—v,—fi i 
q( u >~) c o 

u — J U,—fM >~> — U,—fM • 

The subscript ±/x denotes the direction of the staple, the superscript specifies 
the plane \iv and the parity in v. 
For a line of length 2 we have 



(B.ll) 



v fMU 


= V,,u 


+ v»,- 


-u + V- 


-fM,U 


+ v_ 


■fi,—v 




= Vw 


- v, r 


-u + V. 


-fM,U 


- V- 


H,—u 




= Vw 


+ v*- 


- V -V. 


-fM,U 


- V- 


■fi,—v 






- y,,- 


-v-V- 


-fM,U 


+ v_ 


■H -V 



B.2. The Offset r = (0, 0, 0, 0) 
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and analogously for longer lines. In analogy to eq. ( |B.S| ) we introduce the 
completely (anti)symmetric combinations V^ ym \ V^l for the line products. 
We also introduce the "4d plaquette" 

Qh,i2,h,U = U([h,h, h, h, —h, ~h, ~h, ~h]) , (B.12) 

and the odd combination Qi~i~ 2 ~i~i\- 

To simplify the expressions below, we require that all the directions entering 
Ph h< &h h< Vh h < etc - are different, i.e. they are taken to be zero if e.g. 

\h\ = \hl 

Below we list terms for the different choices of Fq and Iq used in current 
implementation of the parametrized FP Dirac operator code. The notation Yl 
indicates that all indices in the corresponding sum are taken to be different. 



B.2 The Offset r = (0, 0, 0, 0) 

r = i, k = [} 

1. (B.13) 

r = i, i = [1,2, -1,-2] 

r = 7l, i = [1,2,-1,-2] 

^E'>( P ^ +) - h - c -) • (B.15) 

flU 

F = i<Ti2, lo = [1,2,-1,-2] 

i^m^f). (B.16) 

pL<U 

r = 75, h = [1,2,-1,-2,3,4,-3,-4] 

^75 E'lW^^- (B.17) 

flispcr 

r = 75, lo = [1,2,3,4,-1,-2,-3,-4] 

^775 E e ^P^Q^po r ~ ) ■ ( B - 18 ) 

fiupa 
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r = 7175, = [1, 2, -1, -2, 3, 4, -3, -4] 

7^ E {Pfr-lPP + h-c) . (B.19) 

fiupa 

r = 7175, lo = [2, 1, -2, -1, 3, 4, -3, -4] 

77^ E (^; +) ^ s) + h-c.) . (B.20) 

flfpCT 

B.3 The Offset r = (1, 0, 0, 0) 

r = i,io = [i] 

E^ + ^-m)- (B.21) 

r = l, i = [2, i, -2] 

^E'(* +)+5 -t ) )- (B.22) 
r = l, Z = [2,3,1,-3,-2] 

2i E' + V -P s i' +)V P + h - c -) • (B.23) 

iivp 

r = 7i, «o = [i] 

E>(^"^)- (B.24) 
/* 

r = 7 i- = [2,1,-2] 

*E^(* +) -^ +) )- ( b - 25 ) 

r = 72, / = [1,2,3,-2,-3] 



^E'>(^ +) -^ } ^- h - c -) 



(B.26) 



B.4- The Offset r = (1, 1, 0, 0) 
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r = mi2, l Q = [2,1,-2] 

^'^i(5^-)-5^-)-h.c.) . (B.27) 

hp 

r = ia 23 , k = [1,2,3,-2,-3] 

Yq { V p P ^ + P ^ )V p ~ hx ') • (B.28) 

Hup 

r = 75, io = [2,1,-2,3,4,-3,-4] 

^75 E ' ^ W (4^ } ^ S) + + h ' c -) • (B.29) 

/If per 

r = 7i75, k = [2,1,-2,3,4,-3,-4] 

^ E ' w (4"'" )p iS' ) - + h - c -) • ( B - 3 °) 

flUpcr 

r = 7275, / = [1,3,4,-3,-4] 

^ E ( f ^ s) + p £°^ + h - c ) • (B - 31) 

/If per 

B.4 The Offset r = (1,1,0,0) 
r = i, z = [1,2] 

^ E ^ ym) • (B.32) 

\x<v 

r = 7i, io = [1, 2] 

^E'>(^ +) - h - c -) • (B.33) 

/If 

r = 73, /o = [1,3,2,-3] 

£ E'> (^' _) + ^^ P '" ) " *' _) ^ " " h - c -) • (B.34) 

iivp 
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r = *<7i2, i = [1,2] 

5E*Vt#°- (B.35) 

/^<^ 

r = i<7is, = [1,3,2,-3] 
J E'^p (^ ,_) " + - ^^- P r } - h-c.) . (B.36) 

pup 

r = ia 34 , l = [1,2,3,4,-3, -4] 

Y 2 E , ^(^ ++)j $' ) - h - c -) • (B.37) 
r = 75, = [1,2,3,4,-3,-4] 

^75 E ' ^ W (Vj"^ + h-c) • (B.38) 
r = 7175, h = [1,2,3,4, -3,-4] 

7^ E ' 7^ W + h-c) . (B.39) 

Fo = 7375, k = [1,4,2,-4] 

I E '7p75W " 

+ V-^Sf > + h-c) . (B.40) 

B.5 The Offset r = (1,1,1,0) 
r = i, = [i,2,3] 

I E v iZ m) - ( B - 41 ) 

p<V<p 



B.5. The Offset r = (1, 1, 1, 0) 
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r = 71. to = [1,2,3] 

lE^(^& ++) + ^" ) )- (B.42) 
r = 74, fo = [1,2,4,3,-4] 

YA E ' ^ (^ ym) * _) + ^ ym) ^T } " h - c -) • (B.43) 

\XVpo 

r = *<7i2, = [1,2,3] 

J E'^^ + V ti~ ] - h - c ) • (B.44) 

p,vp 

r = ia 14 , l = [1,4,2,-4,3] 

pupa 

+V ll Sp~W- + V- a Sj?'-\ - V-pSp-lV-a - V-rSp-lV-p - h.c.) . 

(B.45) 

r = 75, lo = [1,4,2,-4,3] 

pupcr 

-VpS^-lV-,, + + V-pSp-lV-r - V-aSp-lV-p + h.c.) . 

(B.46) 

r = 7i75, *o = [1,4,2,-4,3] 

pi^pcr 

- F_ CT 5^-)V^ - V-nS^V-v - V-aS^V-p + h.c.) . 

(B.47) 

ro = 7475, = [1,2,3] 

I E '-VK^W (^T _) + h-c) . (B.48) 

/IJ/pcr 
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B.6 The Offset r = (1,1,1,1) 

r = l, l = [1,2,3,4] 

^^iS? • (B.49) 

r = 7i. = [1,2,3,4] 

^ E'>(^ ++) - h - c -) • ( B - 5 °) 

r =«<7i2, = [1,2,3,4] 

I E'^(^ ++) -h.c) • (B.51) 
r = 75, = [1,2,3,4] 

^75 £ ' ^ + h-c.) . (B.52) 
r = 7175, = [1,2,3,4] 

£ ' WW + h -^-) • (B.53) 

fiupa 



Appendix C 

Parametrizations 



In this appendix we give the parameters of £> par and i? par at different stages 
of the parametrization procedure^]. We give only the parameters of the oper- 
ators used in the most important steps of the parametrization. These are the 
parametrizations of the free operators, which are the starting point of the pa- 
rametrization procedures described in Chapter |3[ the parametrizations on the 
minimized (3 ~ 3.0 configurations and finally we specify the parametrizations 
for £>fp and Rfp obtained from the last step of the parametrization procedure. 
These are the operators which used in all the production runs. 



1 Note that the parameters of the R operator do not exactly correspond to the param- 
eters in the parameter files for the production run code, because the definitions of the 
operators in this code differs slightly from the ones given in this work. The parameters 
given in this appendix, however, are rewritten such that they conform to the conventions 
used in this work. 
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Appendix C. Parametrizations 



C.l Parametrizations for R^p 



offset 


reference path l 


CO 


(0000) 


D 


0.36647189 


(1000) 


[1] 


0.01491600 


(1100) 


[1,2] 


0.00299651 


(1110) 


[1,2,3] 


0.00062104 


(1111) 


[1,2,3,4] 


0.00013253 



Table C.l: The reference paths and the corresponding couplings for all offsets 
of the operator i? par = Rpp in the free case. Note that the determination of 
the parameters of the free Rfp can be done analytically. Together with the fact 
that i?FP is strictly a hypercubic operator this makes that this parametrization 
is exact, except for the truncation of the coefficients. 



offset 


reference path Iq 


CO 


Cl 


(0000) 


D 

[1,2,-1,-2] 


0.363762 
0.002710 


0.000240 


(1000) 


[1] 

[2,1,-2] 


0.011249 
0.003667 


0.001769 


(1100) 


[1,2] 


0.002996 


0.000710 


(1110) 


[1,2,3] 


0.000621 


0.000247 


(1111) 


[1,2,3,4] 


0.000132 


0.000062 



Table C.2: The reference paths and the corresponding couplings for all offsets 
of the operator i? par used on minimized f3 3.0 configurations. It is used 
in connection with a 1-level APE smearing with coefficients c\ = 0.356 and 
c 2 = -0.079 (see Section |J). 



C.2. Parametrizations for SDfp 
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offset 


reference path l 


CO 


Cl 


(0000) 


D 


0.35793322 


0.00105816 




[1,2,-1,-2] 


0.00853867 




(1000) 


[1] 


0.01948398 


-0.00971484 




[2,1,-2] 


-0.00456798 




(1100) 


[1,2] 


0.00299651 


-0.00069294 


(1110) 


[1,2,3] 


0.00062104 


0.00009383 


(1111) 


[1,2,3,4] 


0.00013253 


0.00004382 



Table C.3: The reference paths and the corresponding couplings for all off- 
sets of the operator i? par used in production runs. 



C.2 Parametrizations for ©fp 



offset 


r 


reference path Iq 


CO 


(0000) 


i 


D 


2.68037621 


(1000) 


i 


[1] 


-0.15397946 




7i 


[1] 


0.10830529 


(1100) 


1 


[1,2] 


-0.03798073 




7i 


[1,2] 


0.03545779 


(1110) 


1 


[1,2,3] 


-0.01355349 




7i 


[1,2,3] 


0.01206385 


(1111) 


1 


[1,2,3,4] 


-0.00645569 




7i 


[1,2,3,4] 


0.00427270 



Table C.4: The reference paths and the corresponding couplings for the different 
Tq's and offsets of the hypercubic parametrization of the free massless £>fp- 
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otrset 


To 


reference path 1$ 


CO 


Cl 


(0000) 


1 


D 


1.94363679 


11.19824151 






[1,2,-1,-2] 


0.73666167 


-11.29811978 




7i 


[1,2,-1,-2] 


-0.11480487 


0.10127720 






[1,2,-1,-2] 


-0.04800436 


0.06833909 




75 


[1,2,-1,-2,3,4,-3,-4] 


-0.04851795 


4.32121617 




7i75 


[1,2,-1,-2,3,4,-3,-4] 


0.15675378 


3.07442815 






[2,1,-2,-1,3,4,-3,-4] 


0.44943367 


0.51046076 


(1000) 


1 


[1] 


-0.09051958 


-0.25661340 






[2,1,-2] 


-0.06354147 


0.23665490 




7i 


[1] 


0.01888183 


-0.01941104 






[2,1,-2] 


0.08936156 


0.05080964 




72 


[1,2,3,-2,-3] 


-0.08485598 


-0.17419947 




^12 


[2,1,-2] 


0.02286974 


-0.01987646 




^23 


[1,2,3,-2,-3] 


-0.03402500 


-0.02913557 




75 


[2,1,-2,3,4,-3,-4] 


/-\ h -1 r~ A\ A\ A\ *-v 

0.01152283 


-0.92662206 




7175 


[2,1,-2,3,4,-3,-4] 


-0.06931361 


0.88455911 




7275 


[1,3,4,-3,-4] 


U.Ulo9zUoD 


0. 0b2466o/ 


(1100) 


1 


[1,2] 


-0.03810518 


-0.01157131 




7i 


[1,2] 


0.03481370 


-0.01010100 




73 


[1,3,2,-3] 


-0.00216685 


-0.10400256 






[1,2] 


-0.01946887 


0.15803834 




icii3 


[1,3,2,-3] 


0.00035914 


0.05117233 




icr 34 


[1,2,3,4,-3,-4] 


-0.02584512 


-0.05780439 




75 


[1,2,3,4,-3,-4] 


0.00833421 


*-\ /~ A A\ — 7 A" S~\ 

-0.23649769 




7175 


[~t O O A O A~\ 

[1,2,3,4,-3,-4] 


-0.01619093 


-0.01034868 




7375 


[1,4, 2, -4J 


n mi nv/ioo 


U.U/ 4/ ODDU 


(1110) 


1 


[1,2,3] 


A\ A\ H /—l |— 1— S~\ 

-0.01355233 


-0.00847175 




7i 


[1,2,3] 


0.01221388 


0.01717281 




74 


[1,2,4,3,-4] 


-0.00696166 


A\ A\ A\ /^v r - A' A"\ Jl 

0.08095624 




icri2 


[1,2,3] 


-0.01690089 


A\ A\ Jl - JAlA/ 7J1 

0.04723674 




icri4 


[1,4,2,-4,3] 


0.00098380 


0.04858656 




75 


[1,4,2,-4,3] 


-0.00265353 


0.03418916 




7175 


[1,4,2,-4,3] 


-0.00878439 


-0.01941618 




7475 


[1,2,3] 


0.00614657 


0.19247807 


(1111) 


1 


[1,2,3,4] 


-0.00622566 


-0.00847456 




7i 


[1,2,3,4] 


0.00453820 


0.01531472 




io"12 


[1,2,3,4] 


-0.00525155 


-0.15768691 




75 


[1,2,3,4] 


-0.20323184 


3.11220088 




7175 


[1,2,3,4] 


0.00347359 


-0.23403519 



Table C.5: The reference paths and the corresponding couplings for the different 
To's and offsets of the hypercubic parametrization S) par used on minimized 
« 3.0 configurations. It is obtained from the overlap reparametrization step 



discussed in Section [3.2.4 It is used in connection with a 2-level APE smearing 
with coefficients c\ 



-0.10 and c 2 = -0.70 (see Section gj) 



C.2. Parametrizations for SDfp 
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offset 


r 


reference path l 


CO 


Cl 


(0000) 


1 


D 


2.83651159 


-1.40064002 






[1,2,-1,-2] 


-0.22445666 


1.98366915 




7i 


[1,2,-1,-2] 


-0.12076995 


-0.42014243 






[1,2,-1,-2] 


-0.03056282 


-0.07057469 




75 


[1,2,-1,-2,3,4,-3,-4] 


-0.03056051 


0.31628438 






[1,2,3,4,-1,-2,-3,-4] 


0.64742088 


0.40015556 




7i75 


[1,2,-1,-2,3,4,-3,-4] 


0.15555202 


-0.60960766 






[2,1,-2,-1,3,4,-3,-4] 


0.03179317 


1.99545374 


(1000) 


1 


[1] 


-0.17409523 


-0.02912516 






[2,1,-2] 


0.02792864 


-0.04200253 






[2,3,1,-3,-2] 


-0.00981706 


0.13878744 




7i 


[1] 


0.06080162 


-0.10246365 






[2,1,-2] 


0.04196345 


0.16430626 




72 


[1,2,3,-2,-3] 


-0.07123397 


-0.09855541 




2C"12 


[2,1,-2] 


0.00687342 


0.05725091 




2023 


[1,2,3,-2,-3] 


-0.04266141 


-0.06541919 




75 


[2,1,-2,3,4,-3,-4] 


0.04175240 


-0.66992860 




7175 


[2,1,-2,3,4,-3,-4] 


-0.04050605 


0.02898812 




7275 


[1,3,4,-3,-4] 


0.01942004 


0.06487760 


(1100) 


1 


[1,2] 


-0.04134473 


-0.01726090 




7i 


[1,2] 


0.03393463 


0.04718603 




73 


[1,3,2,-3] 


0.00194291 


-0.00255123 






[1,2] 


-0.01086183 


-0.00693550 






[1,3,2,-3] 


-0.00116166 


0.05587536 




icr 34 


[1,2,3,4,-3, -4] 


-0.04012812 


-0.00457675 




75 


[1,2,3,4,-3,-4] 


-0.01627782 


0.34101473 




7175 


[1,2,3,4,-3,-4] 


-0.01478707 


-0.11230535 




7375 


[1,4,2,-4] 


0.01667895 


0.14258977 


(1110) 


1 


[1,2,3] 


-0.01219653 


-0.03212035 




7i 


[1,2,3] 


0.01204799 


0.02121544 




74 


[1,2,4,3,-4] 


0.00081979 


0.04339593 




«CT12 


[1,2,3] 


-0.00948029 


-0.06624570 




iau 


[1,4,2,-4,3] 


-0.00300787 


0.05319220 




75 


[1,4,2,-4,3] 


-0.00900855 


0.14650243 




7175 


[1,4,2,-4,3] 


-0.00253705 


-0.08197206 




7475 


[1,2,3] 


0.00112681 


0.11754621 


(1111) 


1 


[1,2,3,4] 


0.00119452 


-0.04279686 




7i 


[1,2,3,4] 


0.00613421 


0.00906434 






[1,2,3,4] 


0.00474526 


-0.07095205 




75 


[1,2,3,4] 


-0.23018075 


0.71623195 




7175 


[1,2,3,4] 


0.00053879 


-0.07524030 



Table C.6: The reference paths and the corresponding couplings for the 
different To's and offsets of the hypercubic parametrization © par used in the 
production runs. 



Appendix D 



Implementation of the 
Overlap Dirac Operator with 

i?^ 1/2 

In this appendix we give an overview of the details of the implementation of 



Neuberger's construction with a general R, as defined in eqs. (|5.l| ) and (|5.3| ). 
We assume that the operator R from the GW relation eq. fll.lD commutes with 
75, is trivial in Dirac space, local, invertible and hermitian. The main problem in 
numerical simulations is the fact that the general overlap construction, starting 
from a Dirac operator D and the operator R, involves the square root of the 
operator R or the inverse of its square root. This becomes numerically an 
expensive part of the whole construction and therefore one would definitely like 
to circumvent the calculation of this square root as often as possible. In the 
following we give the general ideas, how one can perform calculations with only 
very few applications of the square root of R, if any at all. These ideas that are 
based on a simple change of basis have been successfully been implemented in 
all our codes where the chiral properties of © par or P par had to be improved. 



D.l Change of Basis 

Let us first of all recall the facts about the change of basis in linear algebra, 
before we make use of it in the overlap construction. We denote quantities in 
the changed basis by primes. In matrix form a linear map of the space C n onto 
itself in a certain basis can be written as 

y = Ax, (D.l) 

where A is a n x n-matrix. A change of basis can be expressed also in terms 
of a matrix S G GL(C n ). This matrix describes the old basis © in terms of the 
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new basis 2d'. This makes clear that a vector x under a change of basis becomes 

x' = Sx. (D.2) 

The linear map itself does not depend on the choice of basis and therefore we 
have 

y' = A'x' 

A' = SAS- 1 . (D.3) 

Next we have a look at the effect a change of basis has on hermitian bilinear 
forms. In matrix form a hermitian form in a given basis can be written as: 

f(x,y)={x,y) B = x^By. (D.4) 

Using again the fact that the choice of basis does not have an influence on the 
hermitian form, i.e. 

(x,y) B = (x',y') B ,, (D.5) 

one sees that the transformation behaviour of B is: 

x^By = x'\S- 1 )^B'S- 1 y' 

B' = (S-^BS- 1 . (D.6) 

This means that the standard inner product, where B = 1, under a change of 
basis goes over into 

(x',y') B , =z'W)-y, (D.7) 

which means that the form of the standard inner product is only invariant under 
unitary transformations, i.e. under a change to a different orthonormal basis. 

With the form of the inner product also the form of the norm changes under 
a change of basis. The norm of a vector is defined through the standard inner 
product in an orthonormal basis: 

\\ x \\=VxJx~. (D.8) 
Under a change of basis this amounts to 

||x'|| = ^x^(SS^)- l x>, (D.9) 
which leaves the norm of the vector invariant under the transformation. 



D.2 Application to the Overlap Expansion 

In this section we discuss the use of a specific change of the basis, which will 
be used to get rid of the square roots of R in the overlap construction. First we 
explain the benefits of this change of basis for the operator V and later for $D. 
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D.2.1 Tricks for V 

The overlap formula in terms of D and R is given by: 

V = 1 - A(A^A)- 1 / 2 

A = l + s- {2R) l / 2 D(2R) 1 ' 2 , (D.10) 

where the square root of the operator ^^4[]has to be approximated in a certain 
way, as discussed in Chapter|5|. The approximation schemes used in the literature 
[98, 120, 131-135| have in common that the square root term is rewritten in a 



more or less complicated way as a function f{A^ A) that under a change of basis 
transforms like Sf{A^A)S~ 1 = /(SA^S^SAS' 1 ). This fact can now be used 
to rewrite the overlap Dirac operator in a new basis, with the change of basis 
given by0 

S = (2R)- 1 / 2 . (D.ll) 
In this basis the operator A' takes the form: 

A' = 1 + s — D2R (D.12) 
and the overlap Dirac operator 

V = 1 - (1 + s - D2R)f{{\ + s- DhR){\ + s - D2R)) . (D.13) 

This means that the multiplication of V times a vector x can be done as a 
sequence of the following three steps: 

1. Change the vector x to the new basis: x' = Sx. 



2. Multiply x' with V from eq. ( pl3|) : y' = V'x'. 

3. Change the resulting vector y back to the old basis: y = S^ 1 y' . 

At the first sight, it doesn't seem as we would have gained much, because we 
are still left with a multiplication of (2i?) 1 / 2 and an other one with (2R)^ 1 ^ 2 . 
The crucial point is now that one typically wants to calculate quantities like 
propagators, eigenvectors or eigenvalues. The calculation of eigenvalues anyway 
becomes very simple as the eigenvalues are invariant under a change of basis and 
therefore we can simply calculate the eigenvalues of V without ever having to 
change the basis. For the calculation of eigenvectors and propagators the story 
is a little more involved, but here the crucial point is that Krylov space methods 



1 Even though the discussion is done with the square root of A' A, the statements are 
also valid for approximations to the matrix sign function e. 

2 S = (2-R) 1 ' 2 is also possible as change of basis, but it has the drawback that the inner 
product in this basis is given by \2R)~ 1 y, while for the choice in the text it is x^2Ry, 
which is clearly simpler to handle numerically. 
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involved in the inversion of a matrix or the calculation of eigenvectors do not 
depend on the choice of the basis and therefore the algorithms can simply be 
rewritten in a arbitrarily chosen basis by adjusting the dot products, the norms 
and the matrix vector multiplications. This finally makes that the sequence of 
steps for the calculation of e.g. a propagator is still similar to the one shown 
before, however, step 2 is now replaced by typically hundreds of multiplications 
with V . That one still has to perform steps 1 and 3 becomes negligible in terms 
of computer time. 



D.2.2 Tricks for © 

In the overlap construction £> is defined as 

© = (2 J R)- 1 / 2 (1 - A{A^A)- 1 / 2 ){2R)- 1 / 2 

A = l + s-{2R) l / 2 D(2R) 1 / 2 . (D.14) 

In terms of the basis change matrix S = (2i?)~ 1 / 2 the relation between £> and 
V is: 

© = SVS (D.15) 
and using the representation of V in the changed basis we have 

© = V'S 2 = V'(2RY l . (D.16) 

In this form the multiplication of 35 with a vector involves no calculation of the 
square root of R. The only inconvenience is the inverse of the operator R. 
However, one of the most important and also most expensive numerical tasks is 
the calculation of propagators; and for this particular calculation the calculation 
of the inverse of R is not even needed, because we have 

8T 1 = 2R{V')- 1 = 2R[1 - (1 - D2R)f((l - DhR){l - D2R))] ^ , 

(D.17) 

where the function / is again an approximation to the square root of A' A 
defined in the previous section and therefore we simply have to invert the same 
operator as for the calculation of the inverse of V . It is very important to note 
that here one does actually not perform a change of basis, one just makes use of 
the simpler structure of V. Hence, contrarily to the case of V, the algorithms 
must not be adapted to another choice of basis. 



One of the biggest advantages of writing £> in the form (p.lq) is that its 



massive version (see Section 5.1 for more details) 



1 - —J ©(0) +m(2R)- 1 (D.18) 
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can be written as 

»(m)= (l-^V + m (2R)- 1 =V'(m)(2R)- 1 (D.19) 

and in this form it can be used in the multi-mass Krylov space inverters [186, 



187]. The reason for this is clearly that V'{m) and T>(m), respectively have 



the shift structure which is needed in multi-mass Krylov space inverters [17£|, 
because 

V'{m) = (l - tt) P '(°) + m ( D - 2 °) 
can be rewritten in the form 

V'(m) = a(m)(P' (0)+mM), (D.21) 

with 

vn Tn 

a(m) = l-- and /x(m) = j—^ (D.22) 

and therefore the inverse is 

V\m)- X = a(m)- l (V'(0) + ^(m))" 1 , (D.23) 

which makes the linear shift structure obvious. The corresponding expression 
for is 

©(m)- 1 = a(m)- 1 2R(V'(0) + /x(m))" 1 , (D.24) 
which allows again the use of a multi-mass solver. 



Appendix E 

Details of the Determination 
of the Low-Energy Constant 

E 



In this appendix we give a detailed description of our determination of the low- 
energy constant S in quenched QCD from the technical point of view. We focus 
especially on the calculation of the trace of the subtracted quark propagator, 
which is used to extract E. 



E.l Calculating the Trace 

The finite-volume and finite-mass behaviour of (S) m y : Q (see Chapter ^ for 
definitions) can be extracted from the trace 



(S) m ,v,Q = ^Tr' 



(®(m)2i?)- 1 - \ 



(E.l) 



where the prime denotes the trace without the contribution of the zero modesQ, 
because these cause a divergence in the quantity (S) m> v > Q/m which is a topo- 
logical finite volume effect. The contribution of the zero modes is exactly known 
in the analytical formula for (S) m ,V,Q ( see ec l- ( |7.4| ), where the corresponding 
divergence is given for (iptp) m,v,q) and therefore we can compare the subtracted 
values for {S) m y ) Q/m. The calculation of the full trace in eq. ( |E.1| ) is com- 
putationally too expensive to be done for lattice sizes used in real simulations. 
Hence, one has to resort to other techniques. Stochastic estimators with ran- 
dom sources are the method of choice in cases, where the operator of which one 
wants to calculate the trace is essentially dominated by its diagonal elements. 
This is actually the case for the operator (S)(m)2i?) _1 — 1/2. 

Actually, with zero modes we mean the modes that in the limit m — > become zero 
modes. 
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Let us be a little bit more specific about this method: We take a set of N 



random vectors \rj^),a = 1, . . . , N such that their elements 77] ■ , i = 1 
have the properties 

1 N 

& ]vEw =fo> = ° 



, n 



N 

\'-' x .Y 



lim 4 E ^ a) \- a -^;> = . ( E - 2 ) 



a=l 

where the brackets denote the expectation value. These properties are fulfilled 



by, for example, Gaussian or Z(2) (see [177] for details) random number distri 



butions. Using vectors whose elements are taken out of such a random number 
distribution one can show that trace of a n x n matrix M can be written as 

N , TV n 



1 N n n 

J im 4 E ( E Mu(vt ] rvt ] + E M^Tvf) 



a=l i=l ijLj 
n 1 N n 1 N n 

E ^ ^ w E#N a) +E ^ ^ E^ (B) N 0) = E M - 

j=l a=l i^j a=l i 

(E.3) 

where we have used the properties from eq. QE.2| ) to make the reduction to the 
sum over the diagonal elements of M. 

In order to calculate the trace without the zero mode contributions there are 
two different ways one can proceed. The most obvious way is to calculate the 
zero modes and to project them out of the random vectors, using a projection 
operator that projects orthogonal to the subspace of the zero modes. This 
direct method, however, shows to be rather expensive because the precision to 
which the Dirac operator has satisfy to the GW relation has to be very high 
in order that this method works properly [ 12C ] and therefore, we apply the 



second method, which makes use of the definite chirality of the zero modes of 
a GW Dirac operator (see e.g. Section 1.2.2| ). The crucial point is that one can 



calculate the trace in the chiral sector that has opposite chirality with respect to 
the zero modes. This is achieved by using only random vectors whose chirality 
are opposite to the zero modes. The chiral random vectors from one single 
sector obviously no longer form a complete basis. But the contributions of the 
two chiral sectors are equal, apart from the contributions of the real modes of 
150. In order to show this we use a complete set of eigenstates \tpi),i = 1, . . . ,N 



2 For reasons of simplicity we use the fact that Tr' [(©(m)2i?) 1 — 1/2] is equal to 



E.l. Calculating the Trace 



137 



of "Df^ and define the eigenstates projected to the two chiral sectors 

1 



(1+75)1^), 
(1-75)1^)- 



(E.4) 
(E.5) 



Having defined these projected eigenstates we can evaluate the contributions of 
the two chiral sectors to the trace of the operator M = V(m)~ l — 1/2 

N . N 

^(^+|M|^ + ) = - ^(Vi|M + 75 M + M 75 + M%) = (E.6) 
i=l i=l 

i(TrM + TY 75 M), (E.7) 

JV N 

J2(^\M\i>l) = - J2(i>i\M - 75 M - M 75 + Mt|V») = (E.8) 



: (TVM - Tr 75 M) . 



(E.9) 



Hence the difference between the trace in the two chiral sectors is given by 
Tr 75 M which for M = V(m)~ l — 1/2 is exactly the contribution of the zero 
modes, since the contributions of all the non-real eigenvalues cancel as they 
come in pairs of opposite chirality and the contribution of the modes at A = 2 
are cancelled by the subtraction of 1/2. Thus for GW fermions we can indeed 
take the trace in the chiral sector opposite to the zero modes and then double 
the result at the end to get the full trace in eq. (|E. If) . 

Technically, there remains still a problem, which is that the zero modes 
are only exactly chiral when the GW relation is fulfilled exactly and this is not 
possible for an ultralocal Dirac operator [|82-84]. In order to make use of the 
method for the calculation of the trace discussed above, we rewrite the trace 
from eq. (E.l) as follows 



TV 



(^(m)2R)- 1 



1 



TV 



©(m) t (©(m)2 J R©(m) t ) 1 -- 



Tr' 



Tr' 



®(m) t 2 J R((S(m) + ®(m) t )2i?) 1 -- 



, (E.10) 



where the last line of the rewriting is only exact, when $D(m) satisfies the GW 
relation. We use the form of the trace from eq. ( |E.10[) in our calculation, 

Tr'fX^m) -1 — 1/2] to do the derivation in terms of T>. The equality of the two traces 
can be shown easily using the change of basis used in eq. (D.ll) and the invariance of the 
trace under a change of basis. 

3 Note that the operator T> is normal and therefore a complete set of eigenstates exists. 
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because the zero modes of £>(m) +S)(m)' have definite chirality and therefore 
the divergent contribution of the zero modes^] is completely eliminated from 
the calculation. Furthermore, the insertion of a 1 in form of 2i?(2i?) 1 in 
eq. ( |E.10|) has a purely technical reason. It makes that the operator, which has 
to be inverted, can be brought in such a form that a multi-mass Krylow space 
solver can be used for the inversion (see Appendix 

Forcing the GW relation to hold for a Dirac operator which does not fulfill 
it exactly, however, introduces a systematic error in the calculation of the trace. 
But it is possible to control this systematic error by increasing the accuracy 
of the approximation of the overlap operator to machine precision and then 
compare it to the result of a less accurate approximation. In our calculations we 
choose the accuracy of the appoximation to the overlap operator such that the 
systematic error is roughly 2 orders of magnitude smaller than the statistical 
error of our calculation and therefore can safely be neglected. Figure 7.1 in 
Chapter ^ shows the dependence of the relative difference between a high order 
Legendre expansion and lower order approximations and clearly shows that the 
systematic error introduced by the rewriting of the trace in eq. flE.lCp is under 
control. Let us finally remark that after rewriting the trace it is possible to 
reduce the precision to which the Dirac operator used in the calculations has to 
satisfy the GW relation. This is surely supports the hope that for the practical 
implementation of chiral fermions on the lattice the GW relation has not always 
to be satisfied to machine precision in order perform chiral measurements with 
small systematic errors. 



4 In numerical simulations one always has to use approximations to GW fermions and 
therefore the zero modes are not exactly zero modes, but real modes with an ei genva lue 
very close to 0. The statement that the divergent contribution to the trace eq. (E.1C) is 
eliminated from the calculation of the trace, however, also holds for these approximate 
zero modes. 
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